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ABSTRACT 

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical 
problems. This paper is a study of the science questions that can be addressed using an instrument 
with flexible scheduling that delivers parallaxes at about 4 microarcsec (/zas) on targets as faint as 
V = 20, and differential accuracy of 0.6 /xas on bright targets. The science topics are drawn primarily 
from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM 
PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of 
astrometric precision in modern astrophysics. 

SIM PlanetQuest is currently in the detailed design phase, having completed in 2005 all of the 
enabling technologies needed for the flight instrument. It will be the first space-based long baseline 
Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many 
astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, 
and the study of quasar and AGN nuclei. Using differential astrometry SIM will search for planets with 
masses as small as an Earth orbiting in the 'habitable zone' around the nearest stars, and could discover 
many dozen if Earth-like planets are common. It will characterize the multiple-planet systems that 
are now known to exist, and it will be able to search for terrestrial planets around all of the candidate 
target stars in the Terrestrial Planet Finder and Darwin mission lists. It will be capable of detecting 
planets around young stars, thereby providing insights into how planetary systems are born and how 
they evolve with time. Precision astrometry allows the measurement of accurate dynamical masses 
for stars in binary systems. SIM will observe significant numbers of very high- and low-mass stars, 
providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision 
proper motion measurements, SIM will probe the Galactic mass distribution, and through studies of 
tidal tails, the formation and evolution of the Galactic halo. SIM will contribute to cosmology through 
improved accuracy of the Hubble Constant. With repeated astrometric measurements of the nuclei of 
active galaxies, SIM will probe the dynamics of accretion disks around supermassive black holes, and 
the relativistic jets that emerge from them. 

Subject headings: Extrasolar Planets, Stars, Galaxies, Quasars and Active Galactic Nuclei, Astronom- 
ical Instrumentation 



^ Jet Propulsion Laboratory, California Institute of Tech- 
nology, 4800 Oak Grove Drive, Pasadena, CA 91109; 
Stephen. unwinOj pi. nasa. gov 

^ Jet Propulsion Laboratory, California Institute of Technology, 
4800 Oak Grove Drive, Pasadena, CA 91109 

^ Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218 

* Michelson Science Center, California Institute of Technology, 
770 S. Wilson Ave., Pasadena, CA 91125 

^ United States Naval Observatory, 3450 Massachusetts Avenue 
NW, Washington DC 20392 

® Department of Physics and Astronomy, 6127 Wilder Labora- 
tory, Dartmouth College, Hanover, NH 03755 

Department of Physics and Astronomy, San Francisco State 
University San Francisco, CA 95064 

* Spitzer Science Center, 1200 E. California Blvd., Pasadena, 
CA 91125 

^ Department of Astronomy, The Ohio State University, 140 W. 
18th Avenue, Columbus, OH 43210 

Department of Physics and Astronomy, Georgia State Uni- 
versity, Atlanta, GA 30303 

Van Vleck Observatory, Wesleyan University, Middletown, 



CT 06459 

Columbia University, Pupin Physics Laboratories, 550 West 
120th Street, New York, NY 10027 

University of Virginia, Department of Astronomy, P.O. Box 
400325, Charlottesville, VA 22904-4325 

University of California, Berkeley, 417 Campbell Hall, Berke- 
ley CA 94720 

University of Maryland, Astronomy Department, 0227 CSS 
College Park, MD 20742 

Landessternwarte Koenigstuhl 12, 69117 Heidelberg, Ger- 
many 

University of California at Irvine, Department of Physics and 
Astronomy, Irvine, CA 92697 

University of California at San Diego, Center for Astrophysics 
and Space Sciences, 9500 Oilman Drive, La JoUa, CA 92093 

Space Sciences Laboratory, 7 Gauss Way, University of 
California, Berkeley CA 94720-7450 

Space Science Institute, 4750 Walnut Street, Suite 205, 
Boulder, CO 80301 

Department of Physics and Astronomy, Washington State 
University, Department of Physics and Astronomy, Webster Hall 
1245, Pullman, WA 99164-2814 



2 



Unwin et al. 



1. INTRODUCTION 

Astrometry is perhaps the most fundamental, and old- 
est of all areas in astronomy, and it remains a corner- 
stone of the field for the twenty-first century. Accurate 
distances to astronomical objects are essential for de- 
riving fundamental quantities like mass and luminosity. 
Photographic astrometry in the nineteenth and twen- 
tieth centuries laid the foundation for our understand- 
ing of local stellar populations by i dentifying the inhab- 
itant s of the solar neighborhood (|Gliesel 119691 : iLuvtenI 
[19791) . The Second US Naval Observatory CCD Astro- 
graph Catalog (UCAC2) is a CCD-based survey covering 
most of the sky, with accuracies of 15—70 milliarcse c 
(mas), depending on brightness (jZacharias et al.ll2004D . 
and utilizing the Hipparcos and Tycho-2 reference frame. 
Recent CCD based astrometry over narrow fields has 
achieved a n accuracy of less th an 1 mas in a single mea- 
surement fPravdo ct al. 2005). On still smaller scales, 
iLane &: Mu terspaugh (2004) have demonstrated ~ 16 
microarcsec measurements between the components of a 
0.25-arcsecond binary, using the Palomar Testbed Inter- 
ferometer at 2 /Ltm. CCD parallaxes ii ow achieve typical 
errors of 0.5 mas (jHarris et al.ll2 00^ . Wide-angle as- 
trometry using ground-based o ptical and near-IR in ter- 
ferometers now reaches 20 mas (jHummel et al.lll994[ ). In 
the radio range, very long baseline interferometry (VLBI) 
astrometry of quasars has allowed the creation of a quasi- 
inertial reference fram e, the ICRF (In ternational Celes- 
tial Reference Frame; iMa et afl Il998l ) with wide-angle 
accuracy 0.25 mas. 

Space-based astrometry has brought about a renais- 
sance in the field. The ESA Hipparcos mission, which 
operated from 1989-1993, yielded an astrometric cata- 
log of 118,000 stars down to 12.5 magnitude, with posi- 
tional accuracy of 1 mas for stars brighter than V ^ 11. 
The European Space Agency (ESA) is now developing 
the Gaia mission as a next generation astromet ric survey 
mission (jPerrvman et al.ll2001l ; lPerrvmaiill2002D , which is 
expected to develop a catalog of ~ 10^ stars, with accu- 
racy ~ 20 — 25 microarcsec (/xas) for stars brighter than 
V = 15. 

In this paper, we present an overview of the impact 
of precision astrometry in many fields of astrophysics. 
We use NASA's SIM PlanetQuest mission, hereinafter 
SIM, as a specific example of a space-based facility in- 
strument for astrometry. This mission has been under 
active development since 1 996, based o n concept studies 
made several years earlier (|Shaol[l993l ). A major objec- 
tive of this paper is to show how microarcsecond-level 
astrometry is a powerful tool for 21-st century astron- 
omy. There are observing opportunities for new experi- 
ments with SIM, and this paper is intended as a resource 
for astronomers using precision astrometry in their re- 
search. Although presented in the context of the specific 
capabilities inherent to the SIM design, the topics repre- 
sent very clearly the impact across many areas of astron- 
omy in which precision astrometry plays a fundamental 
role. Most of the science investigations described here are 
drawn from the Key Projects of the SIM Science Team, 
which was selecte d via a NASA Announcement of Op- 
portunity in 2000 (jUnwinI [20051 and references therein). 
The SIM Science Team members are co-authors on this 
paper. 



Recornmend ed by the 1990 NRC Decadal Survey 
(|Bahcalll990( l. SIM PlanetQuest entered its Formulation 
Phase (Phase A) in October 1997 and was approved to 
enter Phase B in August 2003. SIM w as again endorsed 
by t he 2000 NRC Decadal Survey (i McKee fc Tavlod 
l200Ct ) wherein it was assumed that SIM would be com- 
pleted, making it unnecessary to rank it against new mis- 
sion recommendations in that report. Technology devel- 
opment was completed in July 2005 and formally signed 
off by NASA Headquarters in March 2006 after extensive 
external independent review. Having completed nearly 
all of the Formulation Phase (Phase A/B), SIM is ready 
to enter the Implementation Phase, with mature designs, 
well understood schedule and cost, and low technical and 
cost risk. Unfortunately, there is no official launch date, 
since budget pressures on NASA's Science Mission Direc- 
torate have resulted in NASA delaying the Implementa- 
tion Phase. 

This paper covers the expected science contributions 
of SIM but does not describe any of the technical details 
of the instrument or mission. Brief descriptions of the 
instrument itself and the supportin g technologies may 
be found in seve ral technical papers (iLaskinl 120061 : IMml 
l200llShaoll200l . A companion paper (Shao & Nemati, 
in preparation) explains the SIM instrument design, op- 
eration, performance and calibration in more detail. The 
astrometric performance of SIM is based on an hierar- 
chical error budget with more than 1000 terms, and with 
key sets of parameters verified in a series of testbeds de- 
veloped during Formulation Phase. Quoted performance 
numbers are current best estimates from the error budget 
and detailed instrument design. 

The acronym SIM stands for Space Interferometry Mis- 
sion. SIM will be the first space-based Michelson inter- 
ferometer for astrometry. The instrument will operate 
in the optical waveband using a 9-m baseline between 
the apertures. With a global astrometry accuracy of 
3 /ias for stars brighter than V = 20, it will measure 
parallaxes and proper motions of stars throughout the 
Galaxy with unprecedented accuracy. Operating in a 
narrow-angle mode, it will achieve a positional accuracy 
of 0.6 /ias for a single measurement, equivalent to a dif- 
ferential positional accuracy at the end of the nominal 
5-year mission of < 0.1 ^as. This performance is about 
1000 times better than existing capabilities on the ground 
or in space, and about 100 times better than the upcom- 
ing Gaia mission, for differential measurements. Such 
high accuracy will allow SIM to detect and measure un- 
ambiguous masses of terrestrial planets around stars in 
our Galactic neighborhood. 

SIM is a targeted mission which measures the astro- 
metric positions of stars, referencing the measurements 
to a grid of 1302 stars covering the entire sky. Its 
scheduling is highly flexible, in both the order of ob- 
servations, their cadence, and the accuracy of each indi- 
vidual measurement. This contrasts with the Hipparcos 
and Gaia missions, which scan the entire sky according 
to a pre-determined scanning pattern. Many astrome- 
try experiments can make effective use, or in some cases 
require, this pointing capability ~ for instance, searches 
for terrestrial planets (especially in multiple planet sys- 
tems), stellar microlensing events, orbits of eccentric bi- 
nary systems, and variable targets such as X-ray bina- 
ries and active galactic nuclei. Currently, the ICRF, de- 
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fined by the locations of 212 extragal actic radio sources 
(jJohnston at al.lll995l : fMa et al.lll998l ) with most having 
errors less than 1 mas, is the standard frame for astrom- 
etry. SIM is expected to yield an optical reference frame 
at a level of about 3 /las; it will be 'tied' to the ICRF by 
observing a number of radio-loud quasars in common. 

This paper is divided into a number of sections, each 
covering a major area of SIM astrophysics. In §[2] we 
show how SIM can be used to search a large sample of 
nearby solar-type stars for Earth-like planets orbiting in 
the "habitable zone" , and to take a planetary census of 
an even larger population of stars with a variety of spec- 
tral types, ages and multiplicities. Section [3] describes 
a search for planetary systems around young stars with 
ages of 1-100 Myr, which will provide knowledge of the 
evolutionary history of planetary systems. In §3] we ex- 
plore how combining SIM with existing datasets extends 
our sensitivity to long-period planets and allows a very 
complete picture of planetary systems to be made. Sec- 
tion [5] shows how microarcsecond astrometry allows us to 
make fundamental advances in understanding the stellar 
mass-luminosity relation. In §[6] we show that SIM con- 
tributes to a range of problems dealing with the physics 
of 'exotic' stellar objects such as neutron stars and black 
holes, and stars with circumstellar maser emission. 

Sections [71 [3 and [9] cover stellar evolution, Cepheids, 
and the luminosity-age relation in globular clusters and 
constraints on the ages of clusters and the Galaxy. Sec- 
tion [To] explores the dynamics and evolution of our 
Galaxy using tidal streams, and §[TT] explains how astro- 
metric microlensing provides insight into the mass spec- 
trum of dark bodies in the Galaxy. §[12] and [13] cover 
astrometric studies of the dynamical properties of our 
Galaxy, galaxies out to 5 Mpc, and the structure and 
properties of active galactic nuclei. In §[T3] we present 
an application of SIM astrometry to cosmology. Sec- 
tion [15] describes how SIM will be able to make high 
dynamic range and high angular resolution images using 
the technique of aperture synthesis. In addition to its 
scientific significance, such data represent a demonstra- 
tion of the future of high-resolution optical/IR imaging 
in space, similar to the way that the NRAO VLA revo- 
lutionized ground-based radio imaging in the late 1970s. 
In §[T6]we discuss a fundamental physics experiment with 
SIM. 

Astrometric measurement techniques which support 
the science objectives are covered in the Appendices. 
Appendix [Al discusses SIM's 'narrow-angle' precision as- 
trometry mode, which is used for discovering planets 
around other stars and measuring their masses, as well as 
other science. In Appendix[B]we show how SIM measure- 
ments of 'grid stars' are used to construct an astrometric 
reference frame for wide-angle astrometry, and how these 
measurements are tied to an inertial reference frame de- 
fined by ICRF quasars (Appendix [C|). Defining a non- 
rotating frame to high precision is essential for some of 
the science described in this paper. Appendix iP] explains 
the astrophysical criteria used to select suitable stars to 
serve as reference objects. 

We conclude the paper (§[T2l) with a recap of the ma- 
jor science areas, and an indication of new areas where 
the implications of microarcsecond-accuracy astrometry 
have yet to be explored in detail. This paper may serve as 
a guide to those interested in using precision astrometry 



to help further their research interests. SIM is a facil- 
ity instrument, and there will exist opportunities for the 
science community to propose new experiments. About 
half of the total observing time is assigned to the Sci- 
ence Team, but the remaining time is open and is not 
yet allocated. To assist researchers, §[17] includes a table 
showing the overall assignment of mission time, including 
open time for new programs. 

2. THE SEARCH FOR POTENTIALLY HABITABLE 
PLANETS 

Greek philosophers Epicurus, Metrodorus of Chios, 
and Lucretius pondered the possibility that many worlds 
like Earth existed. Aristotle and Plato argued that our 
world was unique. We are now in a position to resolve 
this 2400-year debate with NASA space missions such as 
Kepler, SIM, and the Terrestrial Planet Finder. SIM 
measures three key characteristics of a planet (1) its 
mass, (2) the size and shape of its orbit (semimajor axis 
and eccentricity), (3) the inclination of its orbit - if there 
are multiple planets, this will tell if their orbits are co- 
planar. The planet's mass and orbit determine whether 
it can retain an atmosphere, develop a molten core and 
protective dynamo-generated magnetic field, and harbor 
oceans of liquid water. Such planet characteristics are 
believed to play vital roles in the formation and evolu- 
tion of organic life. 

Since 1995, over 200 exoplanets have been discovered, 
most by using the Doppler technique to monitor the 
gravitationally-indu ced radial velocity (RV) 'wobble' in- 
duced by a planet (jMarcv et al.l 120051 ) . Precise knowl- 
edge of the orbits and minimum plan et masses are given 
in the Catalog of Nearby Exoplanets ijButler et al.ll200^ 
which provides the physical and observational properties 
of known exoplanets orbiting within 3 AU. Twenty multi- 
planet systems have been discovered, spawning theoret- 
ical studies of the interactions between planets, their 
nascent protoplanetary disks, and other planets |Ford| 
20051 : lidafcLin] 120051: [Kiev et all l2005t lAlibert et al.l 



2005l:lChiang fc Murravir2002l : iTanaka &: Wardll200ir 

Most exoplanets found so far are gas giants or ice 
giants, with minimum masses greater than Mnop and 
Jupiter- like radii gleaned from transit observations. The 
planet with lowest minimum mass found thus far by 
the Doppler technique has Ms'mi — 5.9 Mm and a pe- 
riod, P = 1.94 d, orbiting the star, GJ876 (iRiver a et all 
l2005f ). This discovery demonstrated that planet forma- 
tion yields masses below 10 and motivated ques- 
tions regarding the occurrence and properties of rocky 
planets. Indeed, the distribution of masses of the well- 
characterized exoplanets around nearby stars (within 
200 pc) rises steeply toward lower masses, at least as 
fast as dN/dM cx M"! "^. While planet detectabil- 
ity in current RV surveys becomes poor for masses be- 
low the mass of Saturn, the rise toward lower masses 
and the correlation of exoplanets with metal abun- 
dance suggest that planets grow from rocky/icy em- 
bryos toward larger masses. Such growth suggests that 
rocky planets should be at least as common as the gi- 
ant planets, forming f rom the leftover planet esimals in 
a protoplanetarv disk (jColdreich. Lithwick fc' Sari 200^; 
iKenyon &: Bromlevl i20061. The observed semimajor axes 
span a range of 0.02-6.0 AU with a rise observed in 
diV/dloga. This suggests that planets of terrestrial 
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masses may also be found at a wide range of orbital dis- 
tances. 

Re cent simulations of the giant planet formation pro- 
cess ()Benz et al.ll2006[) produce large populations of low- 
mass planets whose growth was halted before they could 
become giant planets. These planets orbit beyond 1 AU 
of the parent star. Their results indicate that for every 
currently known exoplanet, there should be many 'failed' 
giant planet cores with masses smaller than 5 M^. Every 
solar-type star may have one or two low-mass planets. 

Extrapolation of the RV-discovered distribution to 10 
AU and integration over the entire range of semimajor 
axes indicates that at least 10% of all nearby FGK stars 
harbor gas giants in the inner few AU. Given the wide 
range of masses and orbital sizes of known planets, we ex- 
pect that many rocky planets will orbit between 0.1 and 
2 AU, having masses above IMg, just the domain in 
which SIM is uniquely sensitive. Note that the RV tech- 
nique, with a precision of 1 m s~^, cannot detect planets 
of 1 Mgg orbiting near 1 AU, as the RV semi-amplitude 
will be ~0.1 m s~^. Moreover, stellar surface jitter of 1 
m s~^ makes improvement in the Doppler technique un- 
likely. Thus SIM offers a unique opportunity to detect 
Earth analogs, planets of one Earth mass in the habitable 
zones of nearby Solar- type stars. 

Gaia is expected to detect many exoplanets, but for 
individual observations its astrometric precision will be 
many times lower than that of SIM, so its main discovery 
space is that of gas-giant planets (see Fig. [T]). Searching 
the new domain of terrestrial planets will require the 
precision and flexibility offered by SIM, which can select 
the number and timing of observations, along with the 
number of reference stars, allowing a tailored study of 
each target. 

2.1. Astrometric Detection of Terrestrial Planets 

The angular wobble induced in a star by an orbiting 
planet is given by: 



Mo Mp a Ipc 



(1) 



Me lAU D 

where a is the angular semi-amplitude of the wobble in 
/ias, Mp is the mass of the planet, M* is the mass of 
the star, a is the orbital semimajor axis, and D is the 
distance to the system. 

SIM's narrow angle observing mode will allow for a 
single astrometric measurement precision of 0.6 /zas for 
stars brighter than V — 7. Narrow-angle astrometry of 
each target star will be made relative to at least 3 ref- 
erence stars selected to evenly surround the target star 
within 1?5. The reference stars are K giants brighter 
than V = 10, within roughly 600 pc, so that the as- 
trometric 'noise' due to orbiting planets is minimized. 
Radial velocity (RV) observations prior to launch will de- 
tect b rown dwa rf and stellar companions of the reference 
stars ([Frink e t al. 2001). A ten-chop sequence between 
a target and a reference star, with 30 sec integrations 
per chop, will achieve 0.85 /Ltas differential measurement 
precision for V ^ 7 stars, including instrumental and 
photon- limited errors; this is more conservative than a 
scenario individually optimized for each targets, which 
delivers 0.6 /las. 

It is important to clearly define what constitutes astro- 
metric detection of a planet. In this paper, for a star of 



a given mass and at a given distance, we define the effec- 
tive mass sensitivity as the mass of a planet that SIM can 
detect with false-alarm probability (FAP) of 1%, and a 
detection probability of 50%. Effective mass sensitivities 
can be determined for lists of actual planet-search tar- 
get stars using Monte Carlo simulations of detection of 
stellar reflex motion due to Keplerian planet orbits. Ef- 
fective mass sensitivity provides a good metric of SIM's 
planet-finding capability for a given target star, since it 
depends only on assumptions about the SIM instrument 
and the known characteristics of target stars, without 
assumptions about the poorly-known properties of the 
planets under study, e.g., their mass distribution, semi- 
major axis distribution, and frequency of occurrence. 

SIM's capability of detecting planets orbiting in the 
habitable zones of nearby stars for several hypotheti- 
cal planet survey s has been investigated in detail by 
ICatanzarite et al] (|2006f ). In their simulations, each tar- 
get was allocated the same amount of observing time. 
In this Section we present a different approach, in which 
each target star is searched to a specific mass sensitiv- 
ity. Thus the simulated observing program computes 
observing time based on the mass and distance of each 
star individually. This takes advantage of recent Micro- 
Arcsecond Metrology (MAM) testbed results at JPL, in- 
dicating that SIM's systematic noise floor is below 0.1 
/ias after many repeated measurements, opening up the 
possibility of detecting sub-Earth mass planets around 
the closest stars. The current best estimate of SIM's 
single-measurement accuracy is 0.6 fxas. Both the num- 
ber of 'visits' during the 5-year mission (nominally ~ 100 
1-D measurements in each of two orthogonal baseline ori- 
entations) and the spacing of those visits, are flexibly 
scheduled, allowing foUowup of the most interesting tar- 
gets during the mission. Approximately 40% of SIM's 
five-year mission time is available for planet searching in 
narrow-angle mode. 

2.2. The Lowest Mass Planets Detectable by SIM 

The threshold planet mass detectable by SIM, for a 
given orbital radius around a star of given mass and dis- 
tance, can be estimated in several ways. ISozzetti et al.l 
(|2002| ) used a chi-squared-based test of the null hypoth- 
esis for detection. They derived a detection threshold of 
S = 2.2, where S is the "scaled signal", the ratio of the 
angular radius of the astrometric wobble and the astro- 
metric measurement accuracy. This criterion has been 
widely quoted in the planet search community. But it is 
overly simplistic to deem a planet detectable only if the 
amplitude of the angular wobble is greater than the as- 
trometric measurement accuracy. Such an estimate fails 
to account for both the advantage of large numbers of 
observations, iVobs and for the temporal coherence of the 
orbital position. 

In our own Monte Carlo study (jCatanzarite et ahl 
l2006f ). planet detection is accomplished by using a joint 
periodogram, the sum of the periodogram power in the 
astrometric measurements along two orthogonal baseline 
directions, after fitting out a model of position offset, 
proper motion, and parallax. From the simulations we 
derive and validate a more appropriate planet detectabil- 
ity criterion. 

We define SNR as the ratio of the angular wobble to 
the standard error of the observation: 
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Fig. 1. — (a) Detailed view of the discovery space for rocky Earth-like (~ 1 — lOMjg) planets in the habitable zone {~ 0. 7-1.5 AU for a G 
star), for the 'Earth Analog Survey' of 129 stars described in §[2] The small dots represent a theoretical planet distribution lllda fc Linll2005l ) 
for planets of 0.1 — SOOOAfg. In this distribution gas giants have an envelope mass at least 10 times the mass of the central core; terrestrial 
planets initially formed within the ice line (2.7 AU for a solar luminosity star); icy planets formed outside the ice line; and hot Jupiters 
have periods < 0.1 AU. Exoplanets discovered as of early 2007 and with semimajor axes > 0.03 AU are shown as filled circles. Planets in 
our Solar System are labeled with single letters. Labeled curves represent the estimated sensitivity limits of indirect detection methods: 
for radial velocity method (RV at 1 m s~^), and astrometry with SIM and Gaia. The SIM sensitivity in this space is a broad band, defined 
by the three lowest curves (labeled with specific Hipparcos star numbers). The lowest curve shows the 'best' star (as computed from star 
mass and distance); the middle curve represents the median star; and the upper curve shows the least favorable star in the sample. Also 
shown is the eff'ective sensitivity of Gaia for stars at 50 pc, a typical distance for Gaia targets, (b) Detailed view of the discovery space for 
the SIM 'Broad Survey' of 2100 stars, in which a much larger sample of stars is surveyed with less sensitivity than in the 'Earth Analog 
Survey'. Symbols and curves have the same meaning as in (a). 
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SNR 



(2) 



2D 



TABLE 1 

Planet Mass-Limited Surveys with SIM 



where N2D is the number of two-dimensional measure- 
ments, fj is the single-measurement accuracy and the 
factor of V2 is inserted because the measurement is dif- 
ferential. SIM measurements are one-dimensional; each 
target is measured and then re-observed within a day or 
so with a baseline orientation that is quasi-orthogonal to 
that of the first measurement. So the number of two- 
dimension measurements is half the number N of 1-D 
measurements. Replacing N2D with N/2, we have 



SNR^ 



N 



2a 



or SNR = 0.5SVN 



(3) 



in terms of the scaled signal S. It is the SNR rather than 
"scaled signal S" which properly determines detectabil- 
ity; we find that a planet with a SNR of 5.4 is detectable 
half the time at the 99% confidence level. With 200 ob- 
servations, our detectability criterion is eq uivalent to S = 
0.76, a factor of three below that quoted in lSozzetti et al.l 
With 200 observations at 0.85 /las differential accuracy, 
SIM achieves mass sensitivities of w 0.2 M0 at the mid- 
habitable zones of the nearest few targets. 

To illustrate the detection process we show in Figure [2] 
how the joint periodogram, by simultaneously using data 
from the two orthogonal baseline directions, is able to 
reliably detect a planet in the low signal-to-noise regime 
where a simple test would reject it. 

2.3. SIM Planet Surveys 

Discovery and characterization of many Earthlike plan- 
ets is one of SIM's most important scientific objectives. 
Up to half of the SIM observing time will be devoted 
to three planet surveys, each with distinct science objec- 
tives: 

• A "Deep Survey" of up to several hundred stars 
located within 30 pc, for the lowest-mass planets 
detectable by SIM. The Deep Survey is expected 
to yield a significant number of Earth- like planets. 

• A "Broad Survey" of ^ 2100 stars over a variety of 
spectral types, ages, and multiplicities, for planets 
with masses of a few Mq and greater. It will ex- 
plore the diversity of planetary systems, providing 
a more complete picture of planetary systems than 
is possible with, say, RV or direct imaging surveys 
alone. 

• A "Young Planet Survey" of ~ 200 stars with ages 
in the range 1 — 100 Myr. This survey, when com- 
bined with the results of planetary searches of ma- 
ture stars, will allow us to test theories of planetary 
formation and early Solar System evolution. 

In this Section we discuss the objectives of the Deep 
Survey and Broad Survey. See §[3]for details of the Young 
Planet Survey. 

The most effective observing strategy for detecting low- 
mass planets will depend on the fraction of stars expected 
to have terrestrial planets (?7©). We expect results from 
the Kepler mission (planet detection via transits) to in- 
form that strategy. The basic argument is simple: if 





Mass sensitivity 


Number of stars surveyed 


Survey 1 


1.0 Me 


129 


Survey 2 


2.0 M® 


297 


Survey 3 


3.0 M© 


465 



Note. — In a mass-limited survey, observing time on 
each star is calculated to yield the given mass sensitivity 
for a planet at the center of the 'habitable zone', com- 
puted from the star 's distance, estimated mass, and spec- 
tral type (see ii l2.3ll . Stars are rank-ordered in observing 
time, and the resulting number of stars which can be sur- 
veyed to a given mass sensitivity is shown as a function 
of planet mass, for three different sensitivity levels. 



Earths are rare, then SIM should concentrate on a larger 
sample of hundreds of stars to get as much information 
on as many systems as possible. If on the other hand, 
terrestrial planets are common, we would like to probe as 
many systems as possible for potentially habitable plan- 
ets. With this goal in mind, it makes sense to search each 
target to the same mass sensitivity, instead of measuring 
each target to the same accuracy, as w as adopted in our 
previous work (jCatanzarite et al1l2006f) . 

To emphasize this 'mass-limited' approach to the 
search, we term this survey the "Earth Analog Survey" . 
We allocate to each target enough integration time to 
allow a planet of 1 Mq to be detected at the radius of 
its mid-habitable zone (as determined from its spectral 
type). We find that with an assignment of 40% of a 
5-year mission, and a single-measurement accuracy of 
0.6 fia.s, SIM can probe the mid-habitable zone of 129 
stars for 1.0 planets (Table [ij. Essentially, one can 
regard this as the survey yield for the idealized case of 
delta- function distributions for planet mass and planet 
orbit radius (IM® at 1 AU for a G2V star). Although 
unrealistic, this measure of performance avoids having 
to make assumptions about the distributions. (In S I2.4I 
below, we show a second simulation of the planet yield, 
this t ime basing it on distributions from Gumming et al.l 
120071 ) 

The discovery space (planet mass vs. orbit radius) for 
the "Earth Analog Survey" is shown in Figure [1^, with 
the 129 stars filling a band in the lower portion of the 
plot. Distributing the observing time over a larger target 
list allows one to detect more terrestrial planets, albeit 
at higher masses. Table [T] shows the expected SIM yield 
for three different values of the search depth. A survey 
to a sensitivity of would encompass more nearby 

stars than would likely be observed by the Terrestrial 
Planet Finder (TPF) Mission. In each survey, the mass 
sensitivity improves with orbit radius, out as far as orbits 
with periods < 5 yr (see Fig.[T^). Note that these surveys 
are intended to be illustrative; at the time of SIM launch, 
the best available data from all sources will be used to 
design a survey which might represent a combination of 
the approaches in Table [TJ 

In the "Broad Survey", SIM will probe 2100 stars for 
planets. As its name implies, this planetary census in- 
cludes stars of all spectral types (including 0,B,A and 
early F, which are not accessible to RV measurements), 
binary stars, stars with a broad range of age and metal- 
licity, stars with dust disks, evolved stars, white dwarfs. 
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Fig. 2. — Simulation of astrometric detection of a planet with 100 SIM measurements in RA and 100 in Dec, over a 5-year time baseline. 
The planet has a mass of 1.5 M^g and orbits at 1.16 AU from a 1 Mq star at a distance of 10 pc; it was chosen to illustrate an orbit 
close to the limit of detectability. (a) Sky plot showing the astrometric orbit (solid curve) and the SIM measurements with error bars, 
for the observing scenario described in §[2T] (b) and (c), the same data as in (a) but shown as time series along with the astrometric 
signal projected onto RA and Dec. (d) Periodograms of the data plotted in (b) and (c). (e) Joint periodogram of data from RA and Dec 
simultaneously. Note that the planet is not reliably detected in RA or Dec, but is detected with a false-alarm probability (FAP) of well 
below 1% in the joint periodogram. This illustrates the power of the joint periodogram relative to the method which does not use any 
time information. In this example, the data shown in (a) have (reduced) x = 1-22, 
null hypothesis with > 99% confidence. 



slightly less than the x = 1-25 required to reject the 



and stars with planets discovered with RV surveys. Each 
class addresses specific features of the planet-formation 
process: Are metals necessary for giant planet forma- 
tion? Does the number of planets decline slowly with 
time due to dynamical evolution? What is the relation 
between dust disks and planets? 



Using about 4% of a 5-year mission, each star will be 
measured 100 times at 4 fias per measurement. Figure [TJd 
shows the discovery space for the Broad Survey, which is 
expected to yield a large sample of hot, cold, rocky, ice 
giant and gas giant planets, as well as multiple-planet 
systems for tests of planet-formation theories. Orbit so- 
lutions will determine masses and inclinations, and elu- 
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cidate planetary system architecture for multiple-planet 
systems. 

SIM's discoveries will complement future exoplanet 
missions. SIM will complete the planetary system ar- 
chitecture for stars with planets identified by Kepler and 
COROT. Where Kepler and COROT find the rocky plan- 
ets SIM will find the gas giant planets. Furthermore, it 
will provide high quality parallaxes (and thus accurate 
angular diameters) of stars around which planets have 
been detected by transits. SIM resolves the uncertainty 
in determining planetary orbit radii from transits. SIM 
determines the orbit, so that it can provide the time- 
dependent location of the planet in the sky, which is 
critical for any follow-up program, such as TPF. 

2.4. Expected Planet Yields for SIM Surveys 

We estimate the likely yield of planets from SIM ob- 
servations under plausible assumptions regarding their 
frequency of occurrence and distributions as a func- 
tion of mass and orbit radius. Tentative target lists 
have been se lected for the survey of nearb y main se- 
quence stars ()Marcv et al.ll2005t IShaoll2006[ ). Our sim- 
ulations use actual star lists, since catalogs of nearby 
stars are almost coniplete. except for late-type stars 
([Duauennov fc Mavoil 119911) . Target hsts for the sim- 
ulations are derived from an initial list of 2350 stars 
taken from the Hipparcos cata-log, w ith distances of less 
than 30 pc (|Turnbull fc TarteHl2003l ) . We excluded stars 
with luminosity greater than 25 L©, thereby eliminat- 
ing giants from our sample. To eliminate the possibility 
of fringe contamination from a binary companion, we 
applied the following selection rules: stars with a com- 
panion closer than 0'.'4 were excluded; for stars with a 
companion that was separated by 0'.'4 to 1'.'5, both were 
included as target-star candidates if the magnitude dif- 
ference was greater than 1; otherwise, both companions 
were excluded. If the target-star candidate had a wide 
binary companion that was separated by more than 1'.'5, 
the companion was added to the list of target-star can- 
didates. Surviving candidates were rank-ordered by ef- 
fective mass sensitivity. 

Although the sensitivity for planet detection at each 
target is of primary significance in assessing the capabil- 
ity of any proposed planet survey, it is also important 
to understand the yield: how many planets we do we 
expect to find and what is their expected mass distribu- 
tion? As discussed in ji l2.3i planet detection sensitivity 
may be derived from knowledge of instrument perfor- 
mance, the target list, the observing scenario, and the 
available observation time for each star. To predict the 
expected survey yield requires knowledge (or plausible 
assumptions) of mass and orbit distributions of planets 
and their occurrence frequency, for solar-type stars. 

Discov eries from the gold en age of radial velocity (RV) 
surveys (|Butler et al.ll2006( l have given us robust knowl- 
edge of planetary statistics for orbits out to 3 AU and 
masses down to a few Saturns. But the surveys arc in- 
complete for planets on more distant orbits; and though 
a handful of planets with masses in the terrestrial range 
have been discovered, these are in very close-in orbits. 
Though RV is advancing toward detection of Earth mass 
planets orbiting M-stars, terrestrial planets in the hab- 
itable zones of solar-type stars will remain beyond its 
capability, except possibly for a handful of nearby stars 



with extremely low variability such as Alpha CenB. On 
the other hand, information on orbital and mass distri- 
butions and occurrence frequency of terrestrial planets 
around Sun-like stars will be forthcoming from the Ke- 
pler mission in a few years; and COROT will very soon 
yield statistics of Neptune-class planets. 

At the present time we can only estimate planetary 
statistics in the terrestrial mass regime by extrapolation 
from observational results and the expectations of planet- 
formation theorists. To this end we created a simple 
hybrid model based on the power-law mass and period 
distributions derived from the RV observ ational data, ob- 
taine d from surveys of solar-type stars (jCumming et al.l 
|2007() . For simplicity, in our model we assume that 
each star has a maximum of one planet. We extrap- 
olated these power laws to orbits out to 10 AU, and 
to masses down to the terrestrial mass regime. To ac- 
count for the preval ence of 'failed c ores' expected by 
many theorists (e.g.. Ilda fc LhJl2005D . we increased the 
occurrence frequency of terrestrial planets by a factor of 
five. The model distributions are depicted graphically in 
Figure [T] Some recent studies suggest tha .t the planet 
occurrence rate is lowe r in low-mass stars (iButle r et al.l 
2004 120061: lEndl et al.l [200l iJohnson et al. 2007), and 



Gould et al.l (|2006D deduce that about a third of low- 



mass star s may have cold Neptu nes, whereas extrapola- 
tion from lCumming et al.l ()2007f ) indicates that only 5% 
of solar-type stars have Neptunes at all separations. 

Our hybrid power-law model has the following prop- 
erties for solar-type stars: 73% of stars have terres- 
trial mass planets (0.3 — lOM^); ~ 10% of stars have 
terrestrial-mass planets in the habitable zone (0.7 — 1.5 
AU); 5% of stars have Neptune-class planets (lOMg to 
O.lMjup); and 16% of stars have Jupiter-class planets 
(0.1 — lOMjup). According to this model, the overall oc- 
currence frequency of planets is 95%. Our predictions of 
Neptunes and Jupiters are probably close to reality, since 
they involve little extrapolation from observational data, 
but the terrestrial mass planet prediction is sensitive to 
our extrapolation. 

Using this hybrid power-law model, we estimate planet 
yields for the 'Earth Analog' and 'Broad surveys' via 
Monte Carlo simulation. For each survey target star, we 
generate 1000 planets, with masses and periods drawn 
randomly from the model described above. For each 
planet we generate a circular orbit, with parameters 
other than mass and period randomized. We calculate 
the reflex motion trajectory of each target star due to 
its planet, and sample it 100 times uniformly over a 
time baseline of 5 years. This results in a time series 
of 100 pairs of simulated RA and Dec true star posi- 
tions. This database of planets and orbits is then stored 
away. Next we create 1000 'sky realizations'; each real- 
ization results from assigning to 95% of the target stars 
a planet randomly drawn from the database; according 
to the statistics of our model, 5% of the targets have no 
planet. Finally, we generate a simulated survey for each 
sky realization by perturbing each stellar reflex motion 
trajectory with parallax, and single-measurement error of 
0.6^as. We pre-process the simulated observational data 
by fltting out a model of position, parallax and proper 
motion, running the flt residuals through the joint peri- 
odogram (see §[2T2]) with the detection threshold set to 
allow only a 1% chance of false detections. Each simu- 
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lated survey therefore has a set of 'input' planets, and 
for each, a subset of those are detected with SIM. The 
most useful representation of the results are histograms 
of the ensemble averages of input and detected masses 
for the simulated surveys and planets. 

Figure [3^ shows the expected histogram for input vs 
detected terrestrial planet masses in the 'Earth Analog' 
survey. The histogram shows fractional counts because 
it is a mean over 1000 simulated surveys. Results for the 
complete range of planet masses are shown in Table O 
In the habitable zone, SIM would detect 61% of all the 
terrestrial planets, including almost half all planets with 
masses in the range 1 — 1.5M0, and nearly every planet 
of higher mass. 

We repeated the methodology described above, with 
the same hybrid power-law distributions for the input 
planets, on the SIM 'Broad Survey' of 2100 stars. In 
Figure[3}D we show the mean of the mass histograms (log- 
arithmic mass bins extending over entire mass range) for 
1000 simulated surveys. Table [3] shows that we expect 
SIM to find 7% of the terrestrial planets, 2% of all terres- 
trial planets in the habitable zone, 47% of the Neptune- 
class planets and 87% of the Jupiter-class planets. 

It is important to realize that the planet yields pre- 
dicted by these simulations depend on many parameters, 
e.g., the SIM single measurement accuracy, the observing 
scenario and time devoted to each target, the mass and 
orbit radius distributions of the planets, and of course, 
the frequency of occurrence of those planets. In partic- 
ular, we note that the fractions of habitable-zone ter- 
restrial planets which are input to the simulations are 
different in Tables [2] and [3l due to different characteris- 
tics of the survey stars. The 'Broad Survey' target list 
includes a larger number of low-mass stars; about half 
have masses < 0.5Mq. Though our hybrid power-law 
model is derived from observations of solar-type stars, 
we have assumed that it also applies to low-mass stars. 
One feature of the model is that there is a decrease in 
the number of planets per dex as the orbit radius be- 
comes smaller than about 0.7 AU. Since the habitable 
zones of low-mass stars are entirely within 0.7 AU, these 
stars will accordingly have fewer habitable zone planets 
than solar-type stars, and this is reflected in the tables. 

To summarize, SIM will be capable of detecting a sig- 
nificant fraction of the expected population of planets 
for a large sample of stars within 30 pc. As the first 
planned instrument capable of detecting terrestrial plan- 
ets around nearby stars, the planet yield from SIM will 
in fact test the degree to which the above model assump- 
tions are valid. SIM's scientific discoveries will likely re- 
veal the erstwhile hidden regime of rocky planets, and 
make possible the first thorough checks of the predictions 
of current theories of planet formation. 

2.5. Physical Parameters of Habitable Planets 

SIM provides a wealth of planetary astrophysics, in- 
cluding the masses, orbital radii, and orbital eccentric- 
ities of rocky planets around the nearest stars. It will 
also find correlations between rocky planets and stellar 
properties such as metallicity and rotation. 

SIM and TPF/Darwin together, along with Kepler, 
provide a valuable combination of information about 
rocky planets. Each mission brings results that illumi- 
nate a different portion of the multidimensional space 



■ Input terrestrial planets 

- Detected terrestrial planets 

■ Input habitable zone terrestrial pianets 

- Detected tiabitable zone terrestrial planets 
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. Input planets 
- Detected planets 




Fig. 3. — (a) Histogram of the expected yield of terrestrial planets 
for the SIM 'Earth Analog Survey' observing program and assumed 
planetary orbit and mass distributions, and a normalized planet oc- 
currence rate, described in fi \2A\ This histogram is a mean of 1000 
simulated surveys in which geometric parameters of the model or- 
bits were randomized. The mean input distribution of terrestrial 
{M < lOAfg) planets is shown as the dotted curve, and the mean 
number of terrestrial planets detected in the survey by SIM as the 
solid curve. For terrestrial planets in the habitable zone only, the 
corresponding curves are shown as dash-dotted and dashed respec- 
tively. The integral planet counts are summarized in Table [21 Note 
that the yield of planets, especially at the low-mass end, depends 
sensitively on both the assumed planetary model parameters and 
the observing strategy, (b) The same as in (a) but for the 2100 
star 'Broad Survey', which includes stars spanning the entire main 
sequence. This survey would discover planets over a wide range 
of masses and orbit radii which are largely unexplored by other 
detection methods. The integral planet counts are summarized in 
Table [3] 

that represents the field of exoplanet research. Kepler 
offers the occurrence rate of small planets. SIM provides 
the masses and orbits of planets around nearby stars, 
identifying the candidate Earths. TPF/Darwin measures 
radii, chemical composition, and atmospheres. In some 
cases, images from TPF/Darwin may provide feedback 
which allows re-analysis of old SIM data, helping orbit 
determination, especially for multiple planet systems. 

Imaging surveys (such as TPF and Darwin) require 
lists of target stars for observation, ideally those for 
which rocky planets have been detected. Assuming that 
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TABLE 2 

Expected yield of SIM Earth-Analog Survey of 129 stars 



Planet type Number Total Number Fraction 

Detected in Sample Detected 

Terrestrial, 0.3 to lOM® 43.0 ± 5.0 94.0 ±5.1 0.46 

Terrestrial (habitable zone) 9.6 ± 2.9 15.7 ± 3.7 0.61 

Ice giant, lOM® to O.lMj^p 5.3 ± 2.2 6.1 ± 2.4 0.87 

Gas giant, 0.1 to lOMjup 21.1 ± 4.2 21.3 ± 4.3 0.99 



Note. — Based on 1000 Monte Carlo survey realizations, as sum ing 
the distribution of planets from the hybrid model discussed in ii l2.4l 



TABLE 3 

Expected yield of SIM Broad Survey of 2100 stars 



Planet type Number Total Number Fraction 

Detected in Sample Detected 



Terrestrial, 0.3 to lOM® 98.3 ± 9.5 1511.6 ± 20.4 0.07 

Terrestrial (habitable zone) 2.4 ± 1.6 154.2 ± 12.0 0.02 

Ice giant, lOM© to O.lMjup 47.1 ± 6.6 99.8 ± 9.6 0.47 

Gas giant, 0.1 to lOMjup 303.5 ± 16.6 347.4 ± 17.4 0.87 



Note. — Based on 1000 Monte Carlo survey realizations, assuming the 
distribution of planets from the hybrid model discussed in i) l2.4l 



the fraction of stars with Earths in the habitable zone, 
?70, is 0.1, SIM will produce a list of target stars for 
TPF enriched by a factor of at least 2 in rocky planets 
between 0.5 and 2.0 AU relative to a TPF-only sample. 
For many of these stars, SIM's orbital solution will be 
precise enough to predict the best timing for a direct 
observation. This information is crucial for direct imag- 
ing, since a planet in the habitable zone can spend much 
of its time hidden in the glare of the parent star. In- 
deed, habitable rocky planets detected by SIM will likely 
reside at angular separations of at least 100 mas from 
the host star. Such tantalizing rocky planets will be- 
come high priority targets for those instruments, both 
on the ground and in space, that can perform high con- 
trast imaging. With sufficiently long integration times 
and on-band, off-band filters, early imaging of Earthlike 
planets around the very nearest stars may be achieved in 
advance of TPF and Darwin. SIM also identifies those 
stars that TPF and Darwin should avoid, notably those 
with large planets near the habitable zone that render 
any Earths dynamically unstable. Of course, SIM also 
detects those Saturn or Neptune-mass planets located at 
2 AU, valuable in themselves for planetary astrophysics. 

As a benchmark, one may assume that at least 10% of 
stars have a rocky planet between 0.5 and 2.0 AU. If so, 
Kepler will find them in its transit survey of stars 400- 
1000 pc away; and SIM is likely to find the first rocky 
planets orbiting in the habitable zones of Sun-like stars 
closer than 30 pc. Although no rocky planets will be 
detected in common between the two missions, SIM could 
detect gas giants orbiting Kepler target stars for which 
rocky planets have been detected via transits. (Multiple 
planet detections by Kepler will likely be rare due to the 
very stringent coplanarity requirement). 

Detections of rocky planets will spawn theoretical work 
about geophysically plausible interior structures for such 
planets. SIM measures planet masses, which is the basic 
physical parameter for any planet. Imaging of Earth- 
mass planets around all stars within 30 pc remains be- 



yond current technical ground-based capabilities as such 
planets are 10^" times fainter than the host star and 
will be 28th magnitude, comparable to the background 
patchwork of high redshift galaxies. Adding to the chal- 
lenge, planets in somewhat edge-on orbital planes will 
spend a significant fraction of the time located within the 
diffraction-limited angle of the host star. The next gener- 
ation of space-based telescopes, represented by TPF and 
Darwin, will have a rich discovery space to explore. SIM 
will pave the way by conducting an inventory of rocky 
planets around nearby stars. 

2.6. The Impact of Starspots on Astrometric Planet 
Detection 

Stellar variability manifests itself in different ways in 
photometric, astrometric, and radial velocity data. In 
this subsection we estimate the expected astrometric cen- 
troid jitter due to variability of the planet-search target 
stars, and assess the impact on astrometric planet detec- 
tion. 

The 30 year record of satellite observations of the 
Sun's photomet ric variability shows an RMS of 0.042% 
(jFrohlichI |2006[ ). Variations on timescales of days to 
decades can be attributed to the evolution and rota- 
tional modulation of magnetic surface phenomena, e.g., 
sunspots and faculae (Wcnzler, Solanki & Krivova 2005). 
In general, photometric variability in a star due to 
starspots introduces noise in measurements of both its 
photocenter and radial velocity. This noise, in turn, im- 
poses limits on the mass of a planet detectable by these 
types of measurements. 

To investigate the size of the effect, we developed a 
simple dynamic sunspot model that accurately captures 
the known behavior of Sun's photometric variations in 
both time and frequency domains (Catanzarite, Shao & 
Law, in preparation). Starspot noise has a 'red' power 
spectral density (PSD), showing strong variation with 
frequency, and our model takes account of this. The 
important frequencies are those associated with the du- 
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ration of a measurement (about an hour), an observing 
campaign (up to a few years) and with the orbital period 
of the planet one is trying to detect. 

We used our dynamic sunspot model to characterize 
the jitter in the radial velocity and in the astrometric 
centroid. For the Sun, we find typical RMS jitter of 
7x 10^^ AU in the astrometric centroid, and 0.3 m s~^ in 
the radial velocity. Because of the shape of the PSD, 
a simple RMS does not adequately represent the noise 
contribution to planet detection. To gauge the impact 
on planet detection, we sampled the centroid and the RV 
signal from our sunspot model once every 11 days (100 
epochs) over three years. From the PSD of the resulting 
time series, we found that the noise level in the centroid 
due to starspots is 4 x 10~^ AU for orbit periods longer 
than 0.6 years, equivalent to the astrometric signal (at 
10 pc) of a O.lMg planet in a 1 AU orbit, and well below 
the sensitivity of SIM at this distance. 

This level of centroid jitter translates (at 10 pc) to an 
astrometric noise of 0.04 /ias, substantially below SIM's 
noise floor of 0.085 /ias achieved with 100 observations 
with a differential accuracy of 0.85 //as (see ^ 12. ip . We 
therefore conclude that if the Sun were at 10 pc, starspot 
noise would not impact the astrometric detection of ter- 
restrial planets with orbit periods longer than 0.6 years. 

Radial velocity measurements of solar-type stars are 
subject to variability due to starspots. The PSD in ra- 
dial velocity is flat in the same region of frequencies, with 
a noise level of 0.2 m s~^, comparable to the signal of a 
IM® planet in a 1 AU orbit. In addition, RV measure- 
ments may also be subject to astrophysical noise from 
other processes involving velocity field fluctuations, such 
as p-modes. For this reason, the estimated RV jitter 
due to starspots is only a lower bound to the noise in RV 
measurements. Astrometry is not affected by these other 
processes, so our dynamic starspot model is a good rep- 
resentation of the dominant astrophysical noise source 
for astrometric planet detection. A detailed discussion 
of our starspot simulations is forthcoming (Catanzarite, 
Shao & Law, in preparation). 

2.7. Detecting Multiple Planets 

Of the planet-bearing stars identified by the RV tech- 
nique, 20 are revealed to have multiple-plan et systems, 
comprising 13% of the sample. ISozzetti et al.l (2003) and 
iFordI (|2006 ) have investigated astrometric orbit fitting for 
multiple planet systems. Our own simulations show that 
with 200 observations, SIM can detect and characterize 
systems with two or three short-period planets as long 
as their periods are well-separated, which should be the 
case if they are in dynamically stable orbits. Gas giant 
companions with periods longer than the mission length 
are hard to detect, because SIM would detect an accelera- 
tion, but not obtain data for a closed orbit (Gould 2001). 
However, SIM can generate valuable statistical data on 
long-period planets, even if the periods are very uncer- 
tain, because these planets are undetectable with RV 
measurements and too faint for direct imaging. In inter- 
mediate cases, combined RV and astrometric data should 
constrain the orbits and make orbit solutions tractable 
(|Eisner fc Kulkarnil[200l ). 

2.8. Early-Mission Detection of Planets 



Precision astrometry requires knowledge of the SIM 
baseline length and orientation. A set of baseline vec- 
tors for each tile is derived as part of the astrometric 
grid solution. Early in the mission, the grid accuracy, 
and the reconstruction of baseline vectors, is relatively 
poor, but it improves rapidly after about 9 months of 
data have been taken. An observing and analysis tech- 
nique termed grid-based differential astrometry (GBDA) 
has been developed to make effective use of early obser- 
vations of planet-search targets. Details of the method 
are given in Appendix [X] To demonstrate the GBDA 
approach, we modeled the detection of the planet orbit- 
ing Tau Boo, previously de tected by the radial velocity 
method (jButler et al.lll997l ). It has a Jupiter-hke planet 
in a 3-day orbit, and an expected astrometric signature 
of 9.0 /ias. The model shows that this planet would be 
readily detected, even with limited baseline knowledge 
from the grid. Thus SIM will be able to make useful de- 
tections of planetary systems very early in the mission. 

3. JUPITER MASS PLANETS AROUND YOUNG STARS 

A SIM "Young Planets" survey, targeted toward 150- 
200 stars with ages from 1 Myr to 100 Myr, will help 
us understand the formation and dynamical evolution of 
gas giant planets. The host stars of the majority of the 
more than 200 exoplanets found to date are mature main 
sequence stars which were chosen based on their having 
quiescent photospheres for the successful measurement of 
small Doppler velocities (<10 m s^^). Similarly, stellar 
photospheres must be quiescent at the milli-niagnitude 
level for transit detections since a Jupiter mass planet 
transiting a solar type star reduces the photometric sig- 
nal by about 1.4%. Since young stars have RV fluctua- 
tions or rotationally broadened line widths of at least 500 
m s~^ and brightness fluctuations of many percent, opti- 
cal RV measurements accurate to < 100 m s~^ or transit 
observations cannot be used to detect planets around 
young stars. The near-IR is more promising, and a num- 
ber of groups are attempting RV observations to improve 
on these limits and find a few "hot Jupiters" within 0.1 
AU. 

A few potentially planetary mass objects have been 
detected at 20-100 AU from young (< 10 Myr) host 
stars by direct, co ronagraphic iniaging e.g., 2MASSW 
J1207334-393254 (iChauvin et al.l WM) and GQ Lup 
(|Neuhauser et al.l 120051 ). However, these companions 
are only inferred to be of planetary mass by compari- 
son to uncertain evolutionary models that predict the 
brig htness of young Jupiters as a function of mass and 
age (IWuchterl fc Tscha rnuteil l2003t iBaraffe et all 120021 : 
iBurrows et air i997). Since dynamical determinations of 
mass are impossible for objects on such distant orbits, 
it is difficult to be sure that these are planets and not 
brown dwarfs. Nor is it even clear that the origin of these 
distant young Jupiters is due to the same formation pro- 
cesses as_planets found closer-in. Multiple fragmentation 
events (Boss 2001), rather than core accretion in a dense 
disk (Ida fc Lin.2005i) , may be responsible for the forma- 
tion of these objects orbiting so far from their star. 

As a result of the selection biases of the radial veloc- 
ity, transit and direct imaging techniques, we know little 
about the incidence of planets around young stars in or- 
bits close to their stars, leaving us with many questions 
about the formation and evolution of gas giant planets. 
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Fig. 4. — Planet mass detection sensitivity for the SIM-YSO survey (solid curve) in Mj versus orbital semi-major axis. Estimated 
capabilities for a large ground- based coronagraph (taken to be a diameter d = 30 m telescope at A = 1.6 ^m operating an angular distance 
of 4X/d) and a near-IR interferometer (85-m baseline at A = 1.6 /^m and an intermediate stellar distance of 50 pc) are shown as shaded 
boxes. Also plotted are the properties of the known radial velocity (RV) detected planets (dots). RV accuracy for YSOs (dashed curve) is 
limited to about 100 m s~^ (§[3ll. Except for RV, the sensitivity limits assume a distance of 140 pc. 



Using Equation 1, a Jupiter orbiting 5.2 AU away from 
a 0.8 Mq star at the distance of the youngest stellar 
associations (1-10 Myr) such as Taurus (140 pc) and 
Chamaeleon would produce an astrometric amplitude of 
44 /las. At the 25-50 pc distance of the nearest young 
stars (10-50 Myr) such as members of the f3 Pic and 
TWHya groups, the same system would have an astro- 
metric amplitude in excess of 100 /las. Moving a Jupiter 
into a 1 AU orbit would reduce the signal by a factor 
of 5.2, or 50 //as for a star at 25 pc and 8 /las for one 
in Taurus. In its narrow-angle mode, SIM will have a 
Single Measurement Accuracy (SMA) of 0.6 /xas (Itr); 
observations made during wide angle observations (Ap- 
pendix [B| will have SMA ~ 4/ias. Thus a search for gas 
giants falls well within SIM's capabilities and forms the 
core of the SIM-YSO program. Figure H] shows SIM's ex- 
pected astrometric accuracy for the SIM-YSO survey as 
a function of planet mass and semi-major axis. Also plot- 
ted is the expected RV accuracy achievable with present 
day infrared echelle spectrometers. Unlike RV surveys, 
SIM will be able to detect Jupiter mass planets at radii 
out to several AU around young stars. 

3.1. Science Goals 

In a SIM survey of 200 young stars, we expect to find 
anywhere from 10-20 (assuming that only the presently 
known fraction of stars, 5-10%, have planets) to 200 (all 
young stars have planets) planetary systems. We have 
set our sensitivity threshold to ensure the detection of 
Jupiter-mass planets in the critical orbital range of 1 to 5 
AU. These observations, when combined with the results 
of planetary searches of mature stars, will allow us to test 
theories of planetary formation and early Solar System 



evolution. By searching for planets around pre-main se- 
quence stars carefully selected to span an age range from 
1 to 100 Myr, we will learn at what epoch and with 
what frequency giant planets are found at the water-ice 
"snowline' where they are expected to form ([Pollack et al.l 
I19961 ). This will provide insight into the physical mech- 
anisms by which planets form and migrate from their 
place of birth, and about their survival rate. 

With these SIM observations in hand, we will have 
data, for the first time, on a series of important ques- 
tions: What processes affect the formation and dynam- 
ical evolution of planets? When and where do planets 
form? What is the initial mass distribution of plane- 
tary systems around young stars? How might planets 
be destroyed? What is the origin of the eccentricity of 
planetary orbits? What is the origin of the apparent 
dearth of companion objects between planets and brown 
dwarfs seen in mature stars? How might the formation 
and migration of gas giant planets affect the formation 
of terrestrial planets? 

Our observational strategy is a compromise between 
the desire to extend the planetary mass function as low 
as possible and the essential need to build up sufhcient 
statistics on planetary occurrence. About half of the 
sample will be used to address the "where" and "when" 
of planet formation. We will study classical T Tauri stars 
(cTTs) which have massive accretion disks as well as 
post-accretion, weak-lined T Tauri stars (wTTs). Pre- 
liminary studies suggest the sample will consist of ~30% 
cTTs and ^--^70% wTTs, driven in part by the difficulty of 
making accurate astrometric measurements toward ob- 
jects with strong variability or prominent disks. The 
second half of the sample will be drawn from the closest. 
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young clusters with ages starting around 5 Myr, to the 
10 Myr thought to mark the end of prominent disks, and 
ending around the 100 Myr age at which theory suggests 
that the properties of young planetary systems should be- 
come indistinguishable from those of mature stars. The 
properties of the planetary systems found around stars 
in these later age bins will be used to address the effects 
of dy namical evolution and planet destruction (|Lin et al.l 
Since we will also measure accurate parallaxes, we 
will have good luminosities for the host stars, and will 
use these to help estimate ages. 

3.2. Astrophysical Noise 

The photospheric activity that affects radial veloc- 
ity and transit measurements affects astrometric mea- 
surements, but, as we will now show, at a level con- 
sistent with the secure detection of gas giant plan- 
ets. From measurements of photometric var iabil- 
ity (|Bouvier fc Bertouil 119891: iBouvier et al.l I1995D plus 



Doppler imaging I Strassmeier fc Ricd 19981 ). T Tauri 
stars are known to have active photospheres with large 
starspots coveri i ig sig nificant portions of their surfaces 
(jSchussler et al.l Il996l) as well as hot spots due to in- 
falling, accreting material (jMekkadenl |1998[ ). These ef- 
fects can produce large photometric variations which can 
significantly shift the photocenter of a star. Using a 
simple mode l for the effect of s tarspots on the stellar 
photocenter (jTanner et al.ll200"7[ ). for a typical T Tauri 
star radius of 3 J?©, we see that the astrometric jitter 
is less than 3 /^as for R-band variability less than 0.05 
mag. Thus, the search for Jovian planets is plausible 
for young stars less variable than about 0.05 mag in the 
visible even without a correction for jitter that may be 
possible using astrometric information at multiple wave- 
lengths. Note that since both the astrometric signal and 
the astrometric jitter scale inversely with distance, there 
is no advantage (from the jitter standpoint) to examin- 
ing nearby stars even despite their larger absolute astro- 
metric signal. Other astrophysical noise sources, such as 
offsets induced by the presence of nebulosity and stellar 
motions due to non-axisymmetric forces arising in the 
disk itself are negligible for appropriately selected stars. 
Finally, it is worth noting that searching for terrestrial 
planets will be difficult until stars reach an age such that 
their photometric variability falls below 0.001 mag and 
the corresponding astrometric jitter below 0.5 /zas. 

3.3. The Sample 

The youngest stars in the sample will be located in well 
known star-forming regions such as Taurus, the Pleiades, 
Sco Cen, and TW Hydra (Tanner et al. 2007) and will 
be observed in Narrow Angle mode, which is capable 
of achieving a single measurement accuracy of 0.6 fias. 
Somewhat older stars, such as those in the /3 Pictoris and 
TW Hydrae Associations, are only 25-50 pc away and can 
be observed with less mission time in Wide Angle mode, 
capable of 4 /ias single measurement accuracy. We have 
adopted the following criteria in developing our initial list 
of candidates: a) stellar mass between 0.2 and 2.0 M©; b) 
ii < 12 mag for reasonable integration times; c) distance 
less than 140 pc to maximize the astrometric signal to 
be greater than 6 /ias; d) no companions within 2" or 
100 AU for instrumental and scientific considerations, 
respectively; e) no nebulosity to confuse the astrometric 



measurements; f) variability AR <0.1 mag; and g) a 
spread of ages between 1 Myr and 100 Myr to encompass 
the expected time period of planet-disk and early planet- 
planet interactions. 

A literature sear ch and precu rsor observing program 
(described in Tan ner et al.ll2007f l was carried out to iden- 
tify and validate stars according to these criteria. The 
observing program included adaptive optics imaging with 
the Palomar 5m, VLT 8m, and Keck 10m telescopes to 
look for M star or brown dwarf companions; RV measure- 
ments with the McDonald 2.7 m and HET telescopes, as 
well as the Magellan telescope to look for M star or brown 
dwarf companions; and photometric observations with 
smaller telescopes to look for variability. The variabil- 
ity program proved to be the most stringent filter with 
roughly 50% of the sample showing photometric disper- 
sion in excess of 0.1 mag. We now have a validated list 
of 75 stars meeting all of the above criteria. More stars 
will be added to the precursor program to bring the total 
up to the desired number of 150 — 200 stars. With the 
available observing time allocated to this program (1,600 
hours), we will be able to make 75-100 visits to each star 
(up to 50 2-D visits) which, spread over 5 years, will 
be enough to identify and characterize up to 3 planets 
per star having periods ranging from less than a year up 
to 2.5 years. For narrow angle targets we will take ad- 
vantage of the natural clustering of young stars to share 
requisite observations of ~ 5 reference stars, typically 
i? = 10 — 12 mag K giants, with multiple (2-8) science 
stars within a 1° radius. With additional observations 
during a 10 year extended mission, it will possible to 
find planets out to 5 AU. 

A secondary goal of the program is put our knowledge 
of stellar evolution on a firmer footing by measuring the 
distances and orbital properties of '--^lOO stars precisely 
enough to determine the masses of single and binary stars 
to an accuracy of 1%. This information is r equired to cal- 
ibrat e the pre-main sequence tracks fe.g.. iBaraffe et al] 
[200l) that serve as a chronometer ordering the events 
that occur during the evolution of young stars and plan- 
etary systems. To accomplish the goals of this program, 
we will observe a few dozen binary T Tauri stars as well 
as stars with gas disks observed (in millimeter lines of 
CO) to exhibit Keplerian rotation. With accurate orbits 
and distances for these systems, it will be possible to de- 
termine accurate stellar properties for comparison with 
stellar evolution models. 

4. HOW UNIQUE IS THE SOLAR SYSTEM? 

SIM is most sensitive to orbits with periods in the 
range of ~ 1 — 5 years (Fig.[T]), with the most sensitivity 
to periods close to the mission length. SIM is well-suited 
to detect Earth-like planet s in the habitable zone a round 
nearby stars (see § [2] and ICatanzarite et ahl I2006D . For 
periods up to 10 years, estimates of orbital parameters, 
including period, can be made, but as the period length- 
ens, the uncertainties grow quickly. In the limiting case, 
SIM can only make a dete ction of the acceleration due to 
a companion ()Gouldll2001l ) . Even though the parameters 
of any one such target may not be well-determined, im- 
portant statistical conclusions can be drawn from an en- 
semble of long-period systems. Independent data, espe- 
cially over a long time baseline, can greatly improve our 
knowledge of long-period companions. For many targets. 
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there wiU be a 10-15 year basehne of RV measurements 
to draw on. 

Combining astrometry from a "SIM quick look" (SQL) 
survey with data from the Hipparcos astrometric mission 
(|ESAlll997 ) would constrain orbits of 100 years or more, 
and this could be done for several thousand Hipparcos 
stars. Below, we show results from a simulation of the 
extraction of planets from the combined dataset. Orbits 
of a 1 Mjup planet can be reliably characterized up to 
periods of about 10 years, 10 Afjup planets up to 80 years, 
and stellar companions up to 320 years. 

4.1. Masses and Periods of Long-Period Planets 

Long-period extrasolar giant planets appear to be rare: 
only 25 have periods above 5 years and just one has 
a period slightly longer than t hat of Jupiter. Taking 
the se lection effects into account . iTabachnik fc Tremaind 
()2002t l estimate that 3% of Sun-like stars should have a 
planet with a period between 2 days and 10 years and a 
mass of 1 — 10 Mjup. For our simulation, we adopt the 
normalization oflSozzettil (I2005D, which is 1. 62 times that 
computed by Tab achnik fc Tremaind (|2002l ). 

We define a Solar System Giant Analog (SSGA) as a 
planet (or planets) whose mass and period fall within the 
range of the giants of our Solar System (i.e., with a pe- 
riod between 12 and 165 yr and mass of 0.05 — lAfjup). 
Such systems may or may not contain lower-mass plan- 
ets, in closer orbits, but the astrometric signatures of 
SSGAs would normally dominate, and would remain de- 
tectable in distant systems for which terrestrial planets 
are below the detection limit. Because the extrasolar 
giant planet period distribution function increases with 
period, systems dominated by giant planets should be 
rather common, and we estimate that 12.6% of Sun-like 
stars could harbor SSGAs. 

We also define a more massive version of the Solar Sys- 
tem Giant Analog, with mass between 1 and 13 Mjup, as 
a Massive Solar System Giant Analog or MSSGA. These 
are predicted to be quite abundant, and of course are eas- 
ier to detect: 20% of the total number of planetary sys- 
tems with periods up to 165 years, and occurring around 
7.9% of single stars. 

4.2. A Survey for Solar-System Analogs 

To identify likely long-period planetary systems we 
combine data from a "SIM Quick Look" (SQL) survey 
with Hipparcos data (|ESA1 Il997f ). This method uses 
the astrometric parameters as determined from a fit to 
the SIM data to predict the position at the Hipparcos 
epoch. Differences between the observed and predicted 
positio ns indicate the presence of a companion (OlliriJ 
l2007bf l. SIM data allow for the determination of the 
seven astrometric parameters of an "acceleration" solu- 
tion (in addition to the two positions, two proper mo- 
tions, and parallax). Additional SIM observations would 
help with reliable extraction of accelerations. In any case, 
the aim of the survey would primarily be to reject the 
main-sequence (MS) binaries that have huge signals. Ei- 
ther way, a SIM survey would produce a sample rich in 
planetary and/or brown-dwarf (BD) companions. 

For truly single stars, the SIM data will be an excellent 
predictor of positions recorded in the twentieth century. 
However, if the star has a companion, the short-term 
proper motion determined from a SQL survey can be very 



different from the center-of-mass motion. For a face-on, 
circular system, the semi-major axis of the orbit, orbital 
speed, acceleration, and the derivative of the acceleration 
are all substantial for nearby MSSGA systems. For a 1 
M© star at a distance of 20 pc, and a lAfjup planet, we 
can show that MSSGAs with periods in the range of 5 to 
160 years are readily detectable. This issue has been well- 
st udied, in the context of FK5 and Hipparcos astrometry, 
bv lWielen et al.l (|2001h and references therein. 

Due to the short observing span with respect to the 
orbital period, SQL data effectively determine the in- 
stantaneous proper motion and acceleration due to or- 
bital motion. The long time baseline t between the 
SQL and Hipparcos epochs allows us to compute a met- 
ric, Axyir) = y^A2(r) -I- A^(t) which is independent of 

phase for circular, face-on orbits (jOUingl [2007br ). The 
Aa;y(r) diagnostic is useful when it exceeds the astro- 
metric error. 

The xt figure of merit is useful in revealing MSSGA, 
BD and MS companions. For orbital periods between 
5 and 320 years, the (reduced) xt values uncover 11%, 
39%, and 73% of the companions in the MSSGA, BD, 
and MS mass range respectively, if we use only the SIM 
data to compute xt- Here we ignore the effects of incli- 
nation and eccentricity, which complicate the character- 
ization of the companion, although they will not lower 
the \l and Aat; v alues very much (jMakarov fc Kaplan! 
l2005HOllinj|2007bf ). Including the available non-SIM as- 
trometry significantly increases the yield to 46%, 90%, 
and 99.8% for the three mass ranges respectively. These 
results indicate that low-mass companions can be effi- 
ciently detected by combining SQL and Hipparcos data. 
We find that the orbits of a 1 or lOMjup MSSGA can 
be characterized up to periods of 10 or 80 years, respec- 
tively, and for stellar companions up to 320 years (jOlhnd 
I2007bl) . The A^^y,;. values are significant up to 1,200 or 
5,000 years for companions with mass 0.08 or 1.0 Mq, 
respectively. 

4.3. Very Long-period Companions 

There is strong evidence that the intrinsic multiplic- 
ity rate due to either stars or planets is close to 100% 
among Hipparcos MS stars (jOllinel l2005l ). The lack of 
cataloged companions is most likely due to selection ef- 
fects. Thus, those systems without signs of binarity in 
a SQL-|-Hipparcos survey are likely to have either sub- 
stellar companions with an unknown period, or very 
long-period stellar companions. An extended SIM sur- 
vey would further explore these poorly-characterized sys- 
tems. 

An extended SIM astrometric survey would be signif- 
icantly more sensitive than the initial Quick Look Sur- 
vey. Applying the A^y analysis presented above to the 
extended SIM survey indicates that the MSSGAs can 
be detected with masses as low as O.lAfjup in 10-yr or- 
bits. The maximum period for which a lAfjup planet can 
be reliably detected is extended by a factor four (to 40 
years) and for lOAfjup by a factor two (to 160 years). 
Thus SQL-|-Hipparcos plus an extended SIM survey can 
uncover a very significant part of the MSSGA parameter 
space. 

Given the importance of accurate pre-SIM astrometry, 
we note that large-scale ground-based photometric sur- 
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veys such as Pan-STARRs will also provide astr ometry 
at the required (sub-mas) lev el (IChambersll2005f) . Also, 
data from the Gaia mission (jPerrvmanl l2002r )"~will help 
explore and characterize SSGAs more fully. 

5. PRECISION M-L RELATION FOR EXTREME STELLAR 

TYPES 

Mass is the most fundamental characteristic of a star. 
It governs a star's entire evolution — determining which 
fuels it will burn, what color it will be, and how long it 
will live. It is crucial to our understanding of stellar as- 
trophysics that we determine stellar masses to high accu- 
racy. Knowing the masses of main sequence stars answers 
basic stellar astrophysics questions such as. What is the 
mass-luminosity relation for the highest mass and also 
the lowest mass stars? What is the initial mass function? 
What is the mass content of the Galaxy and how does it 
evolve? In fact, the dependence of luminosity upon mass 
— the mass-luminosity relation (MLR) — is one of the 
few stellar relations sufficiently fundamental to be appli- 
cable to many areas of astronomy. With the exception of 
the H-R Diagram, it is the single most important "map" 
of stellar astronomy. To answer truly fundamental as- 
trophysical questions about stars, the ultimate goal is to 
determine masses to 1% accuracy, which will allow us to 
challenge stellar models more severely than ever before. 
Because of SIM's exquisitely accurate astrometric capa- 
bilities, coupled with its faint magnitude limit, we can 
develop a well-stocked "toolbox" of MLRs that can be- 
come the standards against which all stars are measured. 

Here, we consider the extreme ends of the main se- 
quence, where SIM will be crucial in making real progress 
in defining the MLR. In the case of the most massive 
stars, SIM's extreme accuracy will allow us to reach fur- 
ther across the Galaxy to pick up the rare O and B type 
binaries needed for mass determinations. For their much 
less massive cousins, the red M dwarfs, SIM's faint limit 
provides the ability to measure the orbital motions of ob- 
jects all the way to the end of the stellar main sequence, 
and into the regime of the substellar brown dwarfs. 

There are two tactics that can be used to pin down the 
most massive and least massive stars — measurements 
of individual systems in the field, and the calibration 
of the so-called third and fourth parameters, metallicity 
and age, by targeting stars in clusters for which those 
quantities are known. 

5.1. Massive O and B Stars 

Massive stars are key contributors to the energy bud- 
get and chemical enrichment of the Galaxy, but little 
is presently known about their masses (see Figure [5|). 
There are only five known eclipsing binaries among the 
O stars that have reasonably well established masses 
(jHarries et al.1119981 ). and this lack of data has seriously 
hindered our understanding of the evolution of massive 
stars. One unknown, for example, is the maximum mass 
possible for a star. Interior models for massive stars pre- 
dict that stable stars can exist with initial masses of 120 
M0, but the most massive object among the five eclipsing 
binaries is only 33Mq. Furthermore, indirect methods 
of estimating mass by comparison with model evolution- 
ary tracks and through spectroscopic diag nostics lead to 
discr epancies as large as a factor of two ()Herrero et al.l 
I2OOOD . SIM will record the photocentric and/or absolute 
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Fig. 5. — The mass-luminosity relati on in 2007, using eclips- 
ing binary data (open circles) from A ndersenI 111 99 II ) and others, 
supplemented with visual, speckle and interferometric binary data 
(filled circles). Model curves for the mass-luminosity relatio n at the 
indicated ages and solar metallicity are shown, from Schall e? et al.l 
III992 I at the higher masses and Baraffe et al. (1998) at the lower 
masses. The empirical fit of iHenrv fc McCarthvi (il993i ) for stars 
with masses 0.08 to 2.0 Mq is indicated with a dotted line. 

orbits of many binaries and by combining this informa- 
tion with spectroscopic data it will be possible to deter- 
mine accurate distances, inclinations, and masses. An 
excellent example is the massive binary HD 93205, which 
consists of an 03V -I- 08V pair in a 6.08 day orbit. SIM 
observations will show a 45 /^as photocentric variation 
that will yield the first accurate mass for a star at the 
top of the mai n sequence (only known to be in the range 
of 32-154 M0; lAntokhina et al.ll2000D . 

SIM will also provide the first accurate masses of 
the evolutionary descendants of massive stars. The 
most massive stars develop strong outflows later in 
life and appear as Wolf-Rayet (WR) stars. SIM mea- 
surements oOh^WRHnary, WR 22 (WN7 + 09III; 
Schwcickh ajdt et al.lfl99i ). wiU show a 250 /ias astromet- 
ric variation through the 80.3 day orbit. These measure- 
ments will provide the mass of this extraordinary object, 
currently estimated to be 55 ± 7 Mq, the most massive 
star known. Intermediate-mass B stars in close binaries 
are believed to suffer extensive mass transfer and mass 
loss during the Roche lobe overflow phase. The best ex- 
ample of this evolutionary stage is the enigmatic binary, 
77 Lyr (jBisikalo et al.l[2b00l ). which consists of a bright, 3 
Mq star losing mass to a 13 Mq star hidden in an ex- 
tensive accretion disk. The astrometric orbital motion of 
the bright component will amount to 820 /las, and SIM 
will provide accurate mass estimates at this key evolu- 
tionary stage. Finally, two other examples of massive 
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Fig. 6. — A zoom in of the mass-luminosity relation in 2007, 
focusing on red dwarfs. Eclipsing binary data are represented by 
open blue points, and visual binary data from the MASSIF Team 
(s olid red poin t s) and others (open red points). The empirical fit 
of'Henrv et al.l l|1999l ) is indicated with a dotted line, as well as the 
range of minimum masses for main sequence stars, depending on 
metallicities. Note the disconnect between some of the eclipsing 
and visual binary points, as well as the need for a revision to the 
fit (toward higher masses) for the lowest mass stars. 

stars with longer periods include HD 15558 (05e) and 
HD 193793 (WR), with periods of 1.2 and 7.9 years, re- 
spectively. At distances of 1.3 and 2.6 kpc, each system 
has a semimajor axis of 5-10 mas, easily within reach of 
SIM. 

5.2. Low Mass M Stars and Brown Dwarfs 

Red dwarfs dominate the solar neighborhood, account- 
ing for at least 70% of all stars, an d represent nearly 
half of the Galaxy's total stellar mass (jHenrv et al.lll999[ ) 
and Figure [6l These stars have sp ectral type M, V = 
9 - 20, and masses 0.0 8-0.50 Mf;, (jHenrv fc McCarthy! 
119931: iHenrv et al.]|1994D . The MLR remains ill-defined 
for M dwarfs, so their contribution to the mass of the 
Galaxy is a guess at best, and the conversion of a lu- 
minosity function to a mass function is problematic. At 
masses less than ^ 0.20Mq an accurate MLR can pro- 
vide a strict test of stellar evolutionary models that sug- 
gest the luminosity of such a low-mass star is highly de- 
pendent upon age and metallicity. Finally, the MLR be- 
low O.IO Mq is critical for brown-dwarf studies because 
accurately known masses can convincingly turn a candi- 
date brown dwarf into a bona fide brown dwarf. 

In recent decades, the masses of red dwarfs have been 
determined using a combination of infrared speckle in- 
terferometry and HST-FGS, and occasionally via radial 
velocity efforts. The number of red dwarfs with accurate 
mass measurements less than 0.20 Mfn has increased from 
four in 1980 (Popper 1980) to 22 (jHenrv et al.lll999l and 
unpublished). The sample of more massive red dwarfs 
in the range O.5OM0 > AI > O.2OM0 has also increased, 
with particular improvement in the quality of the avail- 
able masses. 

SIM is critical for M dwarf systems because they are 
typically faint and do not allow high-precision RV mea- 
surements due to their slow orbital motion and poorly 
separated spectral lines. In addition to accurate orbital 



monitoring, SIM will provide two crucial pieces of infor- 
mation required to reduce mass errors to the 1% level, 
where they become astrophysically interesting: paral- 
laxes and mass fractions. As an example, we examine 
the nearby binary Gl 748, which represents the current 
state-of-the art accuracies for red dwarf masses (2.4%). 
SIM can improve the mass by reducing the error in the 
semi-major axis of the absolute orbit (147.0 ± 0.7 mas) 
by a factor of 18 (to 0.04 mas, or 10 times the nominal 
astrometric accuracy of SIM for Global Astrometry) and 
the error in the parallax (98.06 ± 0.39 mas) by a factor 
of 10. The result would be mass errors of only 0.1%. 

Mass is the best discriminator between stars and brown 
dwarfs. An object's mass determines whether or not tem- 
peratures in the object's core are sufficiently high to sus- 
tain hydrogen fusion — the defining attribute of a star. L 
dwarfs are objects with smaller masses and cooler tem- 
peratures ('-1500-2000K) than those of M dwarfs, but 
no accurate masses of L dwarfs have yet been measured, 
so the models of L dwarfs are completely untested by 
data. Several hundred L dwarfs have been discovered to 
date, and appropriate systems observable at SIM's faint 
limit are being found. One example is GJ 1001 BC at a 
distance of 13 pc. With an orbital period of ~4 years, 
this system is ideally suited to the nominal SIM mission 
lifetime. 

6. THE LATE STAGES OF STELLAR EVOLUTION 

The first discovery of an X-ray binary occurred over 
four decades ago when Scorpius X-1 was detecte d dur- 
ing a New Mexico rocket flight (jGiacconi et al.l [l962). 
Although the nature of Sc o X-1 was not im mediately 
clear, it was not long before IShklovskiil (jl967j ) suggested 
the correct explanation that Sco X-1 is a neutron star 
accreting from a stellar companion. Over the years, in 
excess of 25 X-ray satellites have found hundreds of X- 
ray binaries with neutron star or black hole accretors ex- 
hibiting a rich variety of physical phenomena. Studies of 
these systems allow us to probe the most extreme phys- 
ical conditions in the universe, including magnetic fields 
at the surface s of neutron stars t hat can be in excess 
of 10^2 Gauss (jCoburn et al.ll2002f) . densities in neutron 
star cores that may be a s much as an order of magnitude 
above nuclear densities (jLattimer fc Prakashl l2004) . and 
gravitational fields near black holes and neutron stars 
can provide tests of strong gravity (Psaltis 2004). 

Studies of accreting stellar mass black holes also 
improve our understanding of Active Galactic Nuclei 
(AGN) and quasars. X-ray bina ries with relativistic 
jets are often called microquasars (jMirabel et al.|[r992f) 
because of the similarities between these systems and 
quasars. However, detailed comparisons between the two 
populations are hampered by uncertain distances to the 
microquasars, making parameters like total luminosity 
and jet velocity uncertain. 

While much has been learned about the physics of 
X-ray binaries, it is evident that more precise measure- 
ments of physical properties are required to make further 
progress in testing theory. Some of the properties that 
are the least accessible using current instrumentation, 
such as the source distance (d), proper motion, and bi- 
nary inclination (i), will be readily measured using SIM. 
Here, we discuss some of the issues related to the physics 
of X-ray binaries that SIM will help to address. 
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For these sources as well as for other radio-emitting 
stars, it will be possible to combine SIM's observations 
with Very Long Baseline Interferometry (VLBI) observa- 
tions to place the radio components within an absolute 
reference frame that is accurate to 3 /las. 

6.1. Masses of Neutron Stars 

Neutron stars (NSs) provide a unique opportunity to 
understand what happens to matter as densities are in- 
creased beyond the density of nuclei. Thus, measuring 
the NS equation of state (EOS) has important impli- 
cations for nuclear physics, particle physics, and astro- 
physics, and measuring NS masses, radii, or both pro- 
vide constraints on the EOS. The NSs for which accu- 
rate mass measurements have b een made lie very close to 
the canonical value of 1.4 Mq (|Thorsett fc Chakrabartvl 
[1999). and EOSs with normal matter (neutrons and 
protons) as well as exotic matter (e.g., hyperons, kaon 
condensates, and quark matter) can reproduce this 
mass for a large range of radii (|Lattimer fc PrakashI 
l2004f) . However, more recently, there are indications that 
some systems may harbor higher mas s . 1. 8-2. 5 Mp,, NSs 
(jBarziv et all [20011 Icla^k et all 12001 iNice et all (20n5h 
and confirming these high NS masses by reducing the un- 
certainties would lead immediately to ruling out a large 
fraction of the proposed EOSs. 

SIM will be capable of making precise orbital mea- 
surements for a large number of High-Mass X-ray Bi- 
nary (HMXB) systems. These systems typically have 
O- or B-type companions with ~25 HMXBs being 
brighter than V ~ 15. They also have orbital peri- 
ods (Porb) of days to a couple of years, and their wide 
orbits give large astrometric signatures. Taking esti- 
mates of HMXB p arameters (Pprb, d, and the compo- 
nent masses) from iLiu. van Paradiis fc van den Heuvel 
(|2000l ) as well as more recent literature, we find that 
21 likely NS HMXBs have orbital signatures (the semi- 
major axis of the optical companion's orbit) of Ugig > 
5 lias and 8 HMXBs have Ogig > 40 /las. Detailed sim- 
ulations that account for the optical so urce brightnesses 
(jTomsick. Quirrenbach fc Reffertl[2005[ ) show that SIM is 
expected to be capable of detecting orbital motion for 16 
NS HMXBs (see Figure [7|). 

The most interesting among these 16 systems are those 
for which the projected size of the NS's orlpit (o xSini) 
has already been measured (jBildsten et al.|[l997[ ). The 
five sources for which this is the case are Vela X-1, 
X Per, 3A 0535-^262, GX 301-2, and PSR B1259-63. 
SIM measurements, along with sini, will immediately 
yield a NS mass measurement. Perhaps Vela X-1 is the 
most tantalizing as it is suspected of having an over- 
massive NS. The current NS mass measure ment for Vela 
X-1 i s Mx = 1.86 ± 0.16Mq (1-cr errors) ([Barziv et al.l 
[200l . In 40 hours of narrow-angle SIM observations 
of Vela X-1, it will be possible to mea sure to 3.9% 
(jTomsick. Quirrenbach fc ReffertI [20051 ) . This IS a ma- 
jor improvement over the current mass measurement and 
will be sufficient to determine if Vela X-1 harbors an over- 
massive NS. As our estimate of the astrometric signature 
for Vela X-1 is 9.5 fias, the microarcsecond measurement 
accuracy provided by SIM is critical. 

6.2. Masses of Stellar Black Holes 
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Fig. 7. — The expected astrometric signature from orbital mo- 
tion vs. V-band magnitude for the 16 neutron star and neutron 
star candidate High Mass X-ray Binaries which we expect to have 
large enough signatures to be detected by SIM. The solid line shows 
the lO-g narrow-angle SIM sensitivity limit found from simulations 
IjTomsick. Quirrenbach fc Reffert^i2005 ^ where 40 hour-long obser- 
vations are made, and we have a priori knowledge of the binary 
orbital period. The dashed line shows the 10-cr wide-angle SIM 
sensitivity limit for 20 half-hour-long observations. 

Black holes (BHs) are among the most fascinating ce- 
lestial objects. Stellar BHs in our galaxy accreting from 
a normal star facilitate investigation of disk accretion 
and relativistic jets. The mass is a fundamental prop- 
erty of a stellar BH and has critical implications for 
the evolution of BH binaries. Although BHs and NSs 
are often difficult to distinguish, it is believed that a 
compact object that is more than 3 Mq is probably a 
BH while objects less than 3 M© could be either a BH 
or a NS. Present measured masses of twenty confirmed 
BHs have possible masses ran ging from 3 to 18 Mq 
(iRemillard fc McClintockl [2006h . The masses of stellar 
BHs have large uncertainties due to the unknown orbital 
inclination, parallax, and other systematic errors. 

As an example, our understanding of the nature of 
the galactic BH Cyg X-1 could change substantially de- 
pending on its true distance and the companion's mass. 
The mass of the BH in Cyg X-1 is estimated as 10 Mq 
([Herrero et al.|[r995f ). However, the companion star of 
Cyg X-1 might be undermassive for its early spectral 
type. Additionally, there is a huge range of distance esti- 
mates for the system: Hipparcos measurements place it 
at 1724 ± 1000 pc; VLBI estimates 1400 ± 900 pc; spec- 
tral analysis places it at 2000 — 2500 pc. If the lower 
companion mass and nearest distance are adopted, the 
mass of Cyg X-1 could be as low as 3 Mq. SIM can re- 
fine the mass measurements for X-ray binaries that are 
thought to harbor BHs. Currently, the main uncertainty 
in the component masses arises from uncertainty in the 
binary inclination and distance to the system, quantities 
that will be measured accurately by SIM. For the case 
of Cyg X-1, the ~20 Mq supergiant companion has an 
orbital astrometric signature of 27 fias at a distance of 
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2.5 kpc, or 34 /ias at 2.0 kpc. Thus, the orbit of Cyg X-1 
can be easily resolved by SIM, allowing the binary incli- 
nation to be determined accurately and the semi-major 
axis t o be determined to better than 2% ()Pan &: ShaklanI 
12005^. SIM's measurements, combined with X-ray spec- 
troscopic and photometric observations and VLBl obser- 
vations, will give us a complete physical picture of stellar 
BHs for the first time. 

SS 433 is an HMXB and microquasar with unique rel- 
ativistic b aryonic jets that precess with a 162 day period 
(|Margon|[T984i ). Although there is evidence to suppor t 
the presence of a BH in the system (jHillwig et al.lf200^ . 
the nature of the compact object is still debated. Due to 
the uncertainty in the compact object's mass, the orbital 
astrometric signature is also uncertain, but could easily 
be 10 — 30/^as. Because the object is bright {V — 14), 
SIM narrow-angle measurements are feasible, and a de- 
tection of its orbital motion would allow for a definitive 
answer regarding the nature of the compact object. 

6.3. Formation and Evolution of Black Holes and 

Neutron Stars 

The evolutionary endpoints of massive stars result in 
compact objects, such as white dwarfs, NSs, or BHs. It 
is of great importance to investigate their birth place, 
asymmetric birth kicks, and the path of formation and 
evolution of these compact objects. In theory, a BH can 
be formed in two different ways: a supernova (SN) ex- 
plosion or collapse of a massive star without an energetic 
explosion. These formation mechanisms can be distin- 
guished by unique information contained in the kine- 
matic history. If a galactic BH or a NS is formed in a su- 
pernova event, very often a supernova remnant is nearby. 
It is crucial to use 3-dimensional space velocity measure- 
ments to determine the object's runaway kinematics and 
the galactocentric orbits (see Mirabel & Rodrigucs 200j, 
and references therein) . Unfortunately, the current mea- 
surement precision of transverse velocities is limited. The 
image superposition technique with HST has a precision 
of ~ 1 mas f or proper motion s and an error of 16% for 
velocities ( Mirabel et al.|[T99^ . SIM can provide at least 
two orders of magnitude improvement for kick velocity 
measurements, and can identify associations between a 
compact object and supernova remnant for many X-ray 
binaries. 

For the scenario where a BH is formed without a super- 
nova event, the most important issues include determina- 
tion of its birth place, measurements of its galactocentric 
orbit, and a thorough invest igation of its space environ- 
ment (iPan fc ShaklanI [2005h . So far, only Cyg X-1 pro- 
vides observational evidence of this outcome, because it 
appears to have proper motion in common with the asso- 
ciation OB-3, and there is no nearby supernova remnant. 
For this type of BH formation, key parameters, such as 
inclination, kick velocity, distance and masses, are either 
indirectly known or lack sufScient precision from current 
observations. 

6.4. LMXB Distances and Constraints on Physical 

Parameters 

As most X-ray binaries are too far away for parallax 
measurements, distance measurements for these systems 
are, for the most part, highly uncertain. It is not unusual 
for an X-ray binary's only distance estimate to be based 



on a companion's spectral type and a system brightness, 
and these estimates can be uncertain by a factor of 2 
or more. This leads to uncertainty about many basic 
system parameters such as luminosities, mass accretion 
rates, radii of NSs (Rutledgc et al. 2002), sizes of accre- 
tion disks, and jet velocities (Fender 2006) . 

SIM will be able to measure the distances of Low- 
Mass X-ray Binaries (LMXBs), and we have used 
van Paradijs fc van den Heuvell ()2001f ) as well as 
more recent literature to compile a list of LMXBs. As 
LMXBs tend to be optically faint, the main selection cri- 
terion is V-band magnitude. Although the optical bright- 
ness can be strongly variable for transients, these sys- 
tems spend most of their time in quiescent (i.e., low flux) 
states. There are 27 LMXBs with ^ < 20 for which SIM 
parallax measurements will be feasible. For the bright- 
est few sources {V — 12-13), it will be possible to mea- 
sure distances to accuracies as high as 2% using 1 hour 
of SIM time. While more time will be required for the 
fainter sources, accurate distance measurements will still 
be feasible. For the 19 sources on our list with V > 17, 
we typically expect to obtain distance measurements to 
5% accuracy using, on average, 5 hours of SIM time per 
source. Thus, these systems make excellent use of SIM's 
ability to observe fainter targets. 

The LMXBs on our target list include microquasars 
such as V4641Sgr and GROJ1655-40 for which accu- 
rate distance measurements will provide a test of whether 
their jet velocities actually exceed 0.9 c. Observations of 
Cen X-4 will allow improved constraints on NS radius 
measurements, and we will be able to determine if the 
brightest persistent NS systems (such as Sco X-1 and Cyg 
X-2) and the brightest X-ray bursters (such as 4U 1636- 
536) reach the Eddington limit. 

6.5. Active Stars and Micro-quasars 

There are various types of stars that produce contin- 
uum emission at radio wavelengths including: RS CVn 
binaries, eclipsing Algol-type binaries, X-ray binaries, 
novae, pre-main sequence stars, and micro-quasars. Two 
areas in which a SIM astrometric mission would have a 
significant impact in the study of radio stars are: (1) 
establishing a link between the ICRF and the optical 
reference frame and (2) in the study of the stars them- 
selves, specifically, the location of the radio emission and 
the mechanism by which it is generated. 

Traditionally, links between the radio and optical 
frames have been determined through observations of ra- 
dio stars. At optical wavelengths, the Hipparcos Cat- 
alogue currently serves as the primary realization of 
the celestial reference system. The link between the 
Hipparcos Catalogue and the ICRF was accomplished 
through a variety of ground and sp ace-based efforts 
(jKovalevskv et al.lll997l : lLestrade et"al]|1999[) . The stan- 
dard error of the alignment was estimated to be 0.6 mas 
at epoch 1991.25, with an estimated error in the sys- 
tem r otation of 0.25 mas yr~^ per axis ([Kovalevskv et al"] 
Il997| ). For future astrometric missions such as SIM, the 
link between the ICRF and the optical frame will be 
established through direct observations of the quasars. 
However, observations of a number of radio stars should 
provide a useful check on this important frame tie. 

In a series of radio obser vations made with connected 
element interferometers pohnston. de Vegt fc Gaumd 
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20031 iBoboltz et 011120031 : iFev et al.|[200l iBoboltz et all 
2007( ). positions and proper motions of ~50 radio stars 
were determined in the ICRF. One goal of this program 
was to investigate the current ac curacy of the ICRF - 
Hipparcos frame tie. Most recently. Bobo ltz et alJ (|2007l ) 
compared radio star positions and proper motions with 
the Hipparcos Catalogue data, and obtained results con- 
sistent with a non-rotating Hipparcos frame with respect 
to the ICRF. These studies demonstrate the methods by 
which the optical and radio frames can be linked on levels 
of a few milliarcsec using radio stars. Such a connection 
between a future SIM optical frame and the ICRF will 
require much more accurate VLBI observations in the 
radio, and will take into account phenomena related to 
the orbits of the close binary companions. 

In addition to establishing a link between frames, ob- 
servations of active radio stars performed with ground- 
based VLBI and SIM will greatly enhance our under- 
standing of these objects. Many of the stars emitting 
in the radio are close RS CVn and Algol-type binaries 
with separations <20 mas and orbital periods <20 days. 
SIM will provide unprecedented insight into the process 
of radio emission and mass transfer for such stars. For 
example, the prototype radio star, Algol, is a triple sys- 
tem. From VLBI measurements, it was concluded that 
the two orbital planes of the close and far pairs are per- 
pendicular to each other, ra ther than being co-planar 
(|Pan. Shao &: Colavital [l993l) . The cause of such per- 
endicular orbits in stellar evolution theory is on-going 
Lestrade et~allll993f ). 

In both RS CVn and Algol-type binaries it is unclear 
where exactly the radio emission originates relative to 
the two stars in the system. Competing mechanisms for 
generating radio emission are reviewed in iRansom et al.l 
2002^) and include phenomena such as: gyrosynchrotron 
radiation from polar regions of the active K-giant star 
(|Mutel et al.i ll998'). emi ssion from coronal loops originat- 
ing on the K-giant (.Franciosini et al1ll999f ). and emission 
from active regions near the surface of both stars with 
possible channeling of energetic electrons along intercon- 
necting magnetic field lines (Ransom ct al. 2002). An as- 
trometric mission such as SIM should provide stellar po- 
sitions on the 10 /xas level, and the full three-dimensional 
orbits required to distinguish between the various emis- 
sion mechanisms. 

SIM will also have the flexibility to coordinate observa- 
tions with ground-based instruments such as the VLBA 
to allow the location of the radio emission relative to the 
stars as a function of time, even for the shortest period 
(~1 day) binaries. In addition, studies of the dynamics 
of radio jets from micro-quasars will take advantage of 
SIM's flexible 'Target of Opportunity' scheduling. Most 
micro-quasars are X-ray transients, and when they un- 
dergo their month- to year-long outbursts they become 
millions of times brighter in X-rays, thousands of times 
brighter in the optical, and they often produce observable 
radio jets. SIM observations of an outburst will provide 
the precise absolute location of the compact object and 
accretion disk, which is critical for interpretation of the 
locations and velocities of the jets. 

A final issue regarding radio stars and micro-quasars 
is the establishment of the linear scale sizes of the sys- 
tems through accurate parallax measurements. Existing 
parallax measurements are sometimes in confiict. For ex- 



ample, with Hipparcos, the distance to the micro-quasar 
LS I -|-61 303 was found to be 190 pc; however, VLBI ob- 
serva tions place it at a distance of 1150 pc (jLestradd 
l2000f ). Through accurate parallax measurements, SIM 
will provide the linear scale sizes necessary to relate the 
radio emission to stellar positions and to constrain theo- 
retical models of radio star and micro-quasar emission. 

6.6. Late-type Stars with Maser Emission 

The evolution of stars along the asymptotic giant 
branch (AGB), including Mira variables, semi-regular 
variables, and supergiants, is accompanied by signifi- 
cant mass loss to the circumstellar envelope (CSE). The 
nature of this mass-loss process and the mechanism by 
which spherically symmetric AGB stars evolve to form 
axisymmetric planetary nebulae (PNe) is not well un- 
derstood. 

The circumstellar maser emission (OH, H2O and SiO) 
associated with many AGB stars provides a useful probe 
of the structure and kinematics of the nearby circumstel- 
lar environment. Figure [5] shows a schematic view of the 
inner CSE of a typical AGB star with masers. The var- 
ious maser regions can be studied at radio wavelengths 
with VLBI, while the star itself, the molecular atmo- 
sphere, and the circumstellar dust can be studied using 
long baseline interferometry in the optical and infrared. 
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Fig. 8. — A schematic view of the radial structure of the en- 
velope of a typical AGB star with circumstellar maser emission. 
Interferometry in the optical, near-infrared, and mid-infrared can 
be used to study the photosphere, the molecular atmosphere and 
the circumstellar dust. VLBI at radio wavelengths can be used to 
study the circumstellar SiO and H2O masers. 

While ground-based techniques provide a powerful tool 
to study AGB stars, there are still unanswered questions 
for which SIM could provide crucial information. For 
example: (1) What is the underlying cause of the tran- 
sition of symmetrical AGB stars to asymmetrical PNe 
(e.g., unseen binary companions, non-radial pulsations)? 
(2) What are the positions of AGB stars relative to the 
circumstellar masers within the CSE? (3) What is the 
linear scale size of the CSE? 
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A review of the c urrent research regardin g PNe shap- 
nig IS presented in iBaUck fc Fra^ (|2002f ). Theoreti- 
cal models involve interacting stellar winds, magnetic 
field shaping, astrophysical jets, and unseen companions. 
Observational radio/IR/optical interferometric studies 
probe the inner regions of the prog enitor AGB stars and 
provide evidence for asym metry (|Monnier eTaDHool 

f boltz & Diamond' '2005f). significant magnetic fields 
lemmings ct al. 200(|) , and highly collimated astro - 
physical jets (Imai et al] 120021 : iBoboltz fc Marvelll2005l ). 
Whatever the mechanism for shaping PNe from AGB 
stars, it must be operating at the innermost scales of the 
CSE. This is just the regime that SIM will be able to 
probe. 

Figure [9] illustrates the problem of referencing op- 
tical/IR interferometry to radio interferometry results. 
Shown are the results of a joint Very Large Tele- 
scope Interferometer (VLTI) and Very Long Base- 
line Array (VLBA) stu dy of the Mira variable S Ori 
(IWittkowski et al.]l2007f l. The stellar diameter, repre- 
sented by the dark circle in the center, was measured 
with the VLTI while the circumstellar SiO masers were 
imaged concurrently with the VLBA. The referencing of 
the star to the masers is purely conjectural, however, 
with additional astrometric information from SIM, this 
assumption would become unnecessary. 
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Fig. 9. — The v = l,J=l — SiO maser emission toward the 
Mira variable S Ori at a stellar phase 'I'vis = 0.56, as measured 
by the VLBA. The top panel shows the spectrum formed by plot- 
ting maser intensity versus velocity. The bottom panel plots the 
spatial and velocity distribution of the masers with point shading 
representing the corresponding velocity bin in the spectrum and 
point size proportional to the logarithm of the flux density. The 
dashed circle is based on the mean angular distance of the SiO 
masers from the center of the distribution. The dark circle in the 
center illustrates the angular size of th e continuum photosphere as 
determined from VLTI measurements llWittkowski et al.ll2007fl . 

A similar astrometric problem is demonstrated by 



recent II2O maser observations of disks toward sili- 
cate carbon stars. I n the case of the star V778 Cyg, 
ISzczerba et al.l (|2Q06f ) were able to use Tycho astromet- 
ric data to associate the H2O masers with an unseen 
companion orbiting the carbon-rich AGB star. A simi- 
lar disk has been observed for the carbon star EU And, 
also traced by H2O masers; however, it is impossible to 
determine whether the disk is associated with the AGB 
star with the available astrometric data. With SIM as- 
trometry, such a determination would be routine. 

Finally, SIM will greatly improve the study of AGB 
stars by providing the precise parallax distances essen- 
tial to establishing a linear scale size for the star and 
the various regions of the CSE. Knowledge of these scale 
lengths is important for theories relating to the chemistry 
of CSEs, the formation points of circumstellar masers 
and dust, and the strength of the stellar magnetic field. 
Furthermore, the linear velocity of circumstellar gas as 
traced by maser proper motions has yet to be accurately 
determined for many stars without precise distances. 

7. STELLAR EVOLUTION, EXTRAGALACTIC DISTANCES, 
AND GALAXY FORMATION 

The study of normal galaxy evolution is greatly en- 
hanced not only by SIM projects that target dynam- 
ics and dark matter, but also by those that target the 
distance scale. We show in this Section that better dis- 
tances translate to much more precise information on the 
chemical and age structure of the various stellar popu- 
lations that make up the Galaxy and external galaxies. 
SIM can be used to obtain parallax distances to Galac- 
tic (disk; Population I) clusters. This will complement 
SIM's Population II (halo and thick disk) distance scale 
investigations described in § [51 

A critical step in studying ages and chemical composi- 
tions of stellar populations is to establish a collection of 
standard clusters, mostly Galactic clusters, for which dis- 
tances, reddenings, abundances, and ages will be derived 
with unprecedented accuracy. SIM is critical to this task 
by providing accurate parallax distances. These standard 
clusters can then be used to tightly constrain theoretical 
isochrone sets more stringently than ever before. The 
isochrones, in turn, give ages for clusters and also for 
galaxies via integrated- light models, and precision stud- 
ies of galaxy evolution are the ultimate aims enabled by 
this SI M study . Morphological lookback studies (e.g., 
GEMS: iRix et a l. 2004) and s pectroscopic surveys (e.g., 
COMBO-17: IWolf et al.i[200l dovetail nicely with stel- 
lar populations studies that earlier predicted in a broad 
way what the direct observations are finding. That is, 
spiral galaxies look as if they have had quasi-continuous 
star formation for long epochs, as expected, but ellip- 
tical galaxies, while mostly dead today, have also had 
much more complex star formation histories than one 
would suppose, not too drastically different from spirals 
(jWort hcv 1998). 

A plethora of questions regarding field versus cluster 
environment, chemical evolution, morphological evolu- 
tion remain, and are likely to remain for many years. 
Present age error s intrinsic to isochr one-based models 
are of order 30% (iCharlot et al.|[l996l ). but the clusters 
studied by SIM should allow for increased precision for 
better understanding of galaxy evolution at all redshifts. 
A reasonable goal in this regard is 5% age precision for 
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favorable, well-observed extragalactic stellar populations 
from the next generation of large, ground-based spectro- 
scopic telescopes. 

But an absolute 5% age precision requires a much bet- 
ter grip on the systematics of stellar populations than 
presently exists. With precise distances from SIM, dis- 
tance will no longer contribute significantly to the uncer- 
tainties, and instead other effects will dominate. For ex- 
ample, the uncertainty in heavy element abundance (Z) 
propagates approximately as Slog age = —3/2 Slog Z 
(|Worthev..l994 ) using stellar temperatures as age indica- 
tors, as one is forced to do for integrated-light applica- 
tions. This implies that overall heavy element abundance 
uncertainty be less than 0.02 dex, a goal reached only 
rarely at present, but which should be very common in 
the near future. The de tailed, element- by-element com- 
position also matters. I Wort he vl (|1998D estimates that 
abundance ratio effects need to be tracked and calibrated 
if they induce more than a 7 K shift in stellar tempera- 
ture. Progress on such detailed effects is underway and 
should be available in a few years. Progress on bolomet- 
ric corrections, absolute flux scale, and stellar color-Tcg 
relations can also be expected shortly. What a standard 
cluster set does is to provide tie-down points for stellar 
modelers, which relates intimately to the interpretation 
of high-redshift stellar population studies. 

SIM will measure parallax distances to the Galactic 
clusters (supplemented by the globular cluster distances 
described in § [3]) . The luminosity of the main sequence 
turnoff in the color-magnitude diagram is the best age 
indicator: the one with the smallest errors (Chabovcr 
I1995D and the one that ties most directly to the "fu- 
sion clock" of the hydrogen-burning star. Distance un- 
certainty is currently the dominant uncertainty, and that 
will be removed by SIM (to less than 1% for most indi- 
vidual Galactic clusters, and perhaps 1% for the Globu- 
lar clusters in aggregate). After abundance effects, the 
remaining uncertainty is that of interstellar extinction, 
which may prove to be the dominant uncertainty in the 
end, although progress is being made in that area as well. 

The target clusters were chosen to fulfill the following 
science goals. (1) We would like to see extragalactic stel- 
lar population age estimates with 5% absolute precision, 
at least for "red envelope" galaxies. This requires that 
isochrone sets be calibrated to the standard cluster set, 
and that the clusters themselves have well-determined 
ages (TableS]). (2) Of the many distance-scale issues that 
are benefitted by SIM, another impo rtant one is the sur- 
face brightness fluctu ation method (|Tonrv fc Schneideil 
119881 iMei et al.l l2005f) that can be used to chart local 
galaxy flows and matter distributions. This method de- 
pends directly on the isochrone sets and is thus tied to 
the standard cluster set. (3) In our Galaxy, the stan- 
dard clusters can be tied in to the photometry of the 
rest of the globular and open cluster system in order to 
investigate the chemical and dynamical history of the 
Galaxy. In external galaxies this is filtered through the 
isochrone sets. (4) The cluster ages themselves are im- 
portant. (5) Finally, the clusters were chosen to be as 
massive as possible in order to attempt to populate the 
rarer, post-hydrogen-burning phases of evolution. These 
phases are the current frontier of stellar evolution, and 
the more constraints we can place, the better off we are. 

The primary selection criterion is to cover as much 



TABLE 4 

Clusters Selected for Population Studies 



Cluster 


Distance (kpc) 


E(B-V) 


[Fo/H] 


Age (Gyr) 


NGC 6528 


9.1 


0.6 


-0.2 


12 


Palomar 6 


7.3 


1.5 


~0.0 


12 


NGC 6440 


8.4 


1.1 


-0.3 


12 


CoUinder 261 


2.2 


0.27 


-0.2 


9 


NGC 6791 


4.2 


0.1 


0.4 


8 


Melotte 66 


2.9 


0.2 


-0.4 


7 


NGC 6253 


1.5 


0.2 


0.4 


5 


Messier 67 


0.8 


0.02 


-0.1 


4 


NGC 2420 


2.2 


0.02 


-0.4 


4 


Berkeley 18 


5.8 


0.46 


0.0 


4 


NGC 6819 


2.4 


0.05 


0.1 


2 


NGC 7789 


1.9 


0.22 


-0.2 


1.7 


IC 4651 


0.9 


0.15 


0.1 


1.5 


NGC 2243 


4.5 


0.05 


-0.5 


1.1 


NGC 2477 


1.2 


0.3 


0.0 


1.0 


NGC 6134 


0.9 


0.4 


0.3 


0.9 


Messier 44 


0.2 


0.0 


0.2 


0.7 


NGC 1817 


2.0 


0.33 


-0.3 


0.4 


NGC 2324 


3.8 


0.11 


-0.8 


0.4 


NGC 2099 


1.4 


0.3 


0.1 


0.4 


Note. — This list is given in 


order of decreasin; 


5 age. All 



parameters given are approximate. 

age versus metallicity parameter space as the Galaxy al- 
lows. This makes "oddball" clusters with atypical abun- 
dances very important. For instance, young and metal- 
poor clusters are rare in the Galaxy, so NGC 2243 be- 
comes a very important cluster. The globular clusters 
do not reach to supersolar abundance, so the old, metal- 
rich cluster NGC 6791 becomes a valuable tie-point. The 
clusters to be observed with SIM cover about 1.5 dcx in 
age, and are listed in Table HI 

8. CEPHEIDS IN THE MILKY WAY 

SIM's contributions to Cepheid science are at least 
fourfold: 1) approximate Cepheid distances are fairly 
easily estimated, so that once a variable star is identi- 
fied as a Cepheid it will be useful for Galactic rotation- 
curve studies if a SIM-based parallax and proper motion 
are available, 2) an accurate distance calibration allows 
for an accurate determination of extinction and metal- 
licity effects on the inferred absolute luminosity, 3) the 
physics of the pulsation mechanism (including the myste- 
rious amplitude decline of Polaris) can be studied in great 
detail for those nearby Cepheids where the extinction is 
small and/or well measured (such as in clusters), and 
for Cepheids that are members of binary systems where 
accurate mass measurements can be made, and 4) the 
changing of the color of the Cepheids during its pulsation 
phase (A(F-7) ~ 1.1 AF; [piling (2007b), private com- 
munication) could help calibrate SIM's color-dependent 
astrometric terms, while on a verage, A incre ases with 
period: AF ~ logPdays (e.g.. iBono et al.ll2000D . 

The second and third points are essential for our under- 
standing and usage of Cepheids as extra-galactic distance 
indicators. Currently, our lack of detailed understand- 
ing of the physics of the pulsation mechanism (i.e., the 
calibration of the period-luminosity-color relation) yields 
galaxy dista nces which carry systernatic uncertain ties of 
order ±5% (|Pietrzvnski et al.l ' i2006t iMacri et all [20061. 
A better understanding of the physics would likely re- 
sult in smaller systematic errors, and hence, in an eas- 
ier method for determining accurate distances to a large 
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number of galaxies. For example, it has been claimed 
that 'bump Cepheids' can be used to determine distances 
below the 2% level. This method is based on a detailed 
analysis of the light-profile an d nonlinear pulsation mod- 
els fe.g.. ' Keller fc Wood '2006V 

The Milky Way is the only galaxy for which we can per- 
form detailed three-dimensional dynamical studies be- 
cause all six phase-space parameters can be determined 
for a number of tracers of the gravitational potential. 
Such studies are essential for the interpretation of ve- 
locity fields of external galaxies, especially those at high 
redshift which are used to infer galaxy-formation sce- 
narios. Young stars such as Cepheids (age ^^50 Myr) 
are very sensitive to sma l l- and large-scale perturbations 
of the potential (jMavoil ri974f) . and are thus very use- 
ful to study the dynamical effects of, for example: 1) 
the bar, 2) spiral structure, 3) the Gould Belt, and 4) 
the warp. For such studies, the apparent magnitude 
is not important, just the distance and space velocity. 
Cepheids are useful for these kind of studies because 
they can be identified based o n their periodic signal 
()Metzger, Caldwell fc Schechteii ll998). A t otal of abou t 
900 Galactic Cepheids are currently known (|Welchlll998f ) , 
while only the 200-odd nearest of these sta rs are typ- 
ically used in s tudies of Galactic dynamics (Zhu' 2000"; 
Metzger. Caldwell &: Schcchte r 19_98,; ,Feast_&: W hitelock. 
1997: Po nteTaLl |l997l: iPont. Mavor fc Burkil il99l 
Caldwell fc CoulsonI |1987[ ). The Cepheid sample pro- 
vides a unique opportunity to perform very detailed stud- 
ies of the dynamics of disk galaxies. Many of these 
Cepheids are too distant for Gaia, but are easy targets 
for SIM. Bec ause the binarit y rate amongst Cepheids is 
large (> 80%: ISzabadosl[2003h . it is crucial to monitor the 
Cepheids astrometrically throughout the SIM mission. 

Cepheids in the Milky Way suffer a significant 
amount of extinction (Ay). For example, the near- 
est 180 stars in the sam ple o f Pont and collabora- 
tors (jPont. Mavor fc Burkil [1991 iPont et all 179971 ) have 
Ay = 1.7 ± 1 mag, where the extinction correction is 
uncertain by about 0.1 magnitude. In general, it is hard 
to determine extinction better than to ±0.05 mag for 
stars with Cepheid colors, e ven with the Gaia instru- 
ment suite (jJordi et al.ll2006f ). Currently, the extinction 
is estimated frorn an intrinsic period-color relatio n (e.g., 
iLanev fc Stobi3 [1991 ICaldwell fc Coulsmll fl986h which 
is calibrated on Cepheids in open clusters. If more accu- 
rate methods become available to determine extinction, 
the Galactic relation between period, luminosity, color, 
metallicity, etc., would be very- well calibrated. Cross- 
validation of such a calibration would be available via the 
Cepheids in galaxies with rotational-parallax distances 
(see gll3|): M31 and M33 (employing SIM) and the 
LMC (from Gaia data). 

The 'expanding photosphere' or 'Baade-Wesselink' or 
'Barnes-Evans' or 'infrared surface brightness' method 
has been used for many years to yield 'geometric' dis- 
tances for Cepheids. In this method, one can equate the 
integral of the changing radial velocity of the stellar en- 
velope during the pulsation cycle to observed changes in 
radius. In principle, this method is very accurate be- 
cause the radial velocities can be measured very accu- 
rately, while the radii of nearby Cepheids can be mea- 
sured employing ground-based interferometry or via a 
surface-brightness color relation. However, this method 



also suffers from zero-point issues ([Gieren et al.l l2005l) 
that may depend on, for example, metallicity, period and 
pulsation mode. Thus, a large sample of Cepheids with a 
range of physical properties is required to establish this 
relation firmly. SIM could provide a much better cali- 
bration of this method than Gaia because SIM can reach 
the required distance accuracies at both faint and bright 
magnitudes. 

A final, perhaps philosophical, point is that Cepheids 
are variable stars, and it is through this variation that 
we can learn much more about the internal structure and 
atmospheric physics than for normal stars. For example, 
the confirmation by helioseismology of the standard solar 
model firmly established neutrino oscillations, and hence 
proved that neutrinos are massive. 

9. ACCURATE AGES AND DISTANCES FOR POPULATION 
II OBJECTS 

The metal-poor stars in the halo of the Milky Way 
galaxy were among the first objects formed in our 
Galaxy. These Population II stars are the oldest ob- 
jects in the universe whose ages can be accurately de- 
termined. Age determinations for these stars allow us 
to set a firm lower limit to the age of the universe and 
to probe the early formation history of the Milky Way. 
The age of the universe determined from studies of Pop- 
ulation II stars may be compared to the expansion age 
of the universe and used to constrain cosmological mod- 
els. Globular clusters (GCs) provide the best opportu- 
nity to determine ages of Population II (hereafter Pop 
II) stars, as it is easy to identify the various evolutionary 
sequences in a GC color-magnitude diagram. The main 
sequence turnoff (MSTO) luminosity is the best stellar 
'clock' which can be used to d e termine the absolute ages 
of GCs (e.g. 'Demarquc IgSO"; 'Rood'l990|; IVandenBerd 
1990; Renzini 1991; Chaboyer et al. 1996). 

The theoretical isochrones in Figure \W\ demonstrate 
how age affects the color-magnitude diagram for a cluster 
of stars with uniform age and metallicity. It is immedi- 
ately apparent that the MSTO and sub-giant regions are 
most sensitive to age differences. The MSTO becomes 
redder and fainter as a cluster of stars gets older. Thus, 
in principle one could determine the age from the color 
of the turn-off, independent of distance. 

The predicted colors of MSTO stars are subject to a 
great deal of uncertainty. To quantify this, we ran a 
Monte Carlo simulation to determine the uncertainty in 
the calculation of isochrones. First, distribution func- 
tions for the various input parameters used in a stel- 
lar evolution code (such as the opacities, nuclear reac- 
tion rates, treatment of convection, oxygen abundance, 
etc) are determined based upon the known uncertain- 
ties in the determination of each of the various quan- 
tities. Isochrones are then calculated for a given set 
of input parameters which are drawn at random from 
the specified distribution function. This procedure was 
then repeated 1119 times in order to determine how the 
known uncertainties in the input parameters required 
for stellar evolution calculations affect the theoretical 
isochrones. The results may be represented as confidence 
contours in the color-magnitude diag ram. Full details 
are given m iBiork fc Chaboved (|2006l ). Figure [lOl illus- 
trates the large uncertainty in the predicted colors of 
the MSTO for a given age and metallicity, by comparing 
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the Monte Carlo isochrones to standard isochrones with 
different ages. The uncertainty in the theoretical calcu- 
lation of the MSTO color leads to an error of approxi- 
mately ±2 Gyr (1 a). In contrast, theoretical uncertainty 
in the MSTO luminosity is considerably lower, of order 
±lGyr. 

The largest uncertainty in the determination of glob- 
ular cluster ages based upon the MSTO luminosity is 
the distance scale for Pop II objects. A 1% error in 
the distance leads to a :^ 2% error in the derived age 
(e.g.. lChabover et al.|[l99l . SIM wiU be able to deter- 
mine distances to globular clusters and other stars in 
the halo with unprecedented accuracy, thereby signifi- 
cantly reducing the uncertainty in the derived ages of 
metal-poor stars. Table [5] provides basic data on the 21 
globular clusters that will be observed with SIM. These 
clusters were chosen based upon the following properties: 
(a) distance from the Sun, (b) metallicity, (c) reddening, 
and (d) whether the cluster is thought to belong to the 
Old Halo (OH), Young Halo (YH) or Thick Disk (TD) 
(column 3). This grouping of globular clusters is based 
upon their ki nematics, metallicity and horiz ontal branch 
morphology (|Da Costa fc Armandrofilfl995f ). 

For each globular cluster, we plan to observe approx- 
imately six red giant branch stars with an accuracy of 
7/xas. We will average together the parallaxes of all the 
confirmed members to obtain a final parallax to the clus- 
ter with an accuracy of 4 /j,as, which supported by the 
expected accuracy of the SIM grid of about 3 /zas. This 
will determine the distances to the individual globular 




B-V 

Fig. 10. — Theoretical isochrones with [Fe/H ] = —1. 6 calculated 
using the stellar evolution code described in IBioriTfc Chaboye^ 
II2006I) . showing the change in the position of the MSTO (solid 
circles) when the age of the standard 13 Gyr isochrone is changed 
by ±3 Gyr (dashed curves) . The effects of changes in other pa- 
rameters (e.g., nuclear reaction rates, opacities, treatment of con- 
vection, and oxygen abundance) were evaluated using the Monte 
Carlo simulation described in §(9] and are shown as the ±68% and 
±95% confidence levels in a standard isochrone (solid curves), with 
corresponding MSTO (triangles). 



clusters with an accuracy of 1% to 5% (column 7 in Ta- 
ble [5]), which is a factor of 2 to 10 times better than cur- 
rently achieved. The current uncertainties in the globular 
cluster distance scale is dominated by systematic errors, 
while our distance scale will be dominated by random er- 
rors. This will allow us to average together the age of the 
most metal-poor globular clusters in our sample, thereby 
significantly reducing the uncertainty in the determina- 
tion of the mean age of the oldest globular clusters. 

In order to determine the expected accuracy in our ab- 
solute age estimate for the oldest, most metal-poor glob- 
ular clusters, a Monte Ca rlo simulation was perform ed, 
similar to that outlined in lBjork fc Chaboveil (|2006l ). In 
this simulation, we varied all of the sources of error in our 
age determinations within their expected uncertainties, 
including: the reddening determinations, photometric 
zero-points, parallax uncertainties, uncertainties in the 
exact composition of the stars (helium abundance, oxy- 
gen abundance and iron abundances), and uncertainties 
in the stellar models (including nuclear reaction rates, 
opacities, treatment of convection, model atmospheres, 
diffusion, etc). The distribution function for each of the 
individual input parameters was determined by a careful 
consideration of the expected uncertainties in the vari- 
ous quantities when SIM will deliver its final parallaxes. 
The simulation showed that we will determine the abso- 
lute age of the oldest globular clusters to an accuracy of 
±3%, or ±0.4 Gyr. 

To study the relative age distribution of stars in the 
halo, SIM will observe 60 metal-poor turn-off/subgiant 
branch stars in the field. To illustrate the expected ac- 
curacy of our relative age determinations for the field 
halo stars and the globular clusters in our SIM program, 
we ran a Monte Carlo simulation which allowed for the 
true distance to the object to vary within its current es- 
timated uncertainties, and which took into account the 
uncertainties in the SIM distance determination, redden- 
ing determinations and in the chemical composition of 
the stars. We find that the field stars will have a typical 
uncertainty of ±0.6 Gyr, while the globular cluster ages 
will have an error of ±0.9 Gyr (column 8 in Table [5]). 
From these simulations we conclude that we will be able 
to detect age differences of the order of 1 Gyr between 
various stellar populations. 

10. EXPLORING GALACTIC STELLAR POPULATIONS 
AND DARK MATTER ON GALACTIC SCALES 

N-body simulations of the formation of structure in 
the Universe in the presence of dark matter (and dark 
energy) show galaxies (and all large structures) build- 
ing up hierarchically. The active merging history on 
all scales demonstrated by high resolution. Cold Dark 
Matter (CDM) numerical simulations has had remark- 
able success in matching the observed properties of the 
largest structures in the Universe, like galaxy clusters, 
but are a challenge to reconcile with the observed prop- 
erties of structures on galactic scales. The Milky Way 
and its satellite system is a particularly important labo- 
ratory for testing specific predictions of the CDM mod- 
els, most especially because high accuracy astrometric 
observations enabled by SIM allow definitive tests of dy- 
namical effects specifically predicted by CDM. 

SIM will make possible unprecedented opportunities 
to explore stellar dynamics with a precision that will 
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TABLE 5 
Target Globular Clusters 



NGC 


[Fe/H] 


Group 


Nrr 


E(B - V) 


Dq (kpc) 


ao (%) 


0"agc {%) 


6341 


—2 29 


OH 


25 


0.02 


8.1 


4 


g 


7099 


—2.12 


OH 


10 


0.03 


7.9 


4 


g 


^uyu 


9 nfi 
— z.uu 


YH 


41 




1 n 1 

lU. -L 


o 


1 n 

-LU 


6397 


— 1 95 


OH 




0.18 


2.2 




5 


6541 


— 1.83 


OH 




0.12 


7.4 


3 


g 


6809 


— 1.81 


OH 


10 


0.07 


5.3 


2 


g 


5139 


— 1 62 


OH 


152 


0.12 


5.1 


2 


7 


5272 


— 1 57 


YH 


260 


0.01 


10.0 


5 


10 


6752 


— 1.55 


OH 




0.04 


3.9 


2 


5 


6205 


-1.54 


OH 


3 


0.02 


7.0 


3 


7 


6218 


— 1.48 


OH 




0.19 


4.7 


2 


g 


3201 


-1.48 


YH 


85 


0.21 


5.1 


2 


8 


5904 


-1.29 


OH 


123 


0.03 


7.3 


3 


8 


288 


-1.24 


OH 




0.03 


8.1 


4 


8 


362 


-1.16 


YH 


13 


0.05 


8.3 


4 


8 


6723 


-1.12 


OH 


29 


0.05 


8.6 


4 


9 


6362 


-1.06 


OH 


33 


0.09 


7.5 


3 


8 


6652 


-0.85 


YH 




0.09 


9.4 


4 


10 


104 


-0.76 


TD 


1 


0.05 


4.3 


2 


5 


6838 


-0.73 


TD 




0.25 


3.8 


2 


7 


6352 


-0.70 


TD 




0.21 


5.6 


3 


9 



Note. — Values taken from the lHarri^ l|1996l ) compilation, unless otherwise 
noted. 



Notes on individual columns: 

Group: OH = Old Halo; YH = Young Halo; D = Thick Disk 
IIDa Costa fc Armand"ro3IT995l ') 

Nrr — Number of RR Lyrae stars in the cluster (from compilation by Carney, 
private communication) 

Dq — distance from the Sun; these values are uncertain by ±10%. 

(T^ — percent uncertainty in the parallax, assuming 4.6 ^as accuracy 

Cage - percent uncertainty in the age estimate including contributions due to 

uncertainties in the distance determination, reddening, photometric zero-point 

and metallicity determination. 



allow critical measurements of gravitational potentials 
from Local Group size (see Section [TT|) to Milky Way, 
dwarf galaxy and star cluster scales. We outline in more 
detail below several specific important SIM contributions 
that will bear directly on tests of dark matter (DM) and 
the evolution of galaxies like the Milky Way and its stel- 
lar populations. 

10.1. Probing the Outer Halo with Tidal Tails 

At large distances from the Galactic center (> 30 kpc), 
the stellar distribution is far from homogeneous. Stan- 
dard methods of estimating the depth of the Milky Way's 
gravitational potential using a tracer population whose 
orbits are assumed to be random and well-mixed would 
be systematically biased under these circumstances (e.g., 
I^ncho et al. 2006). However, these inhomogeneities 
themselves are thought to have formed through the in- 
fall and disruption of satellites and hence we actually 
have more information about the stars in these lumps 
than in a truly random sample. For example, we know 
that stars that are clearly part of a stream of debris were 
once all part of the same satellite. We can use this knowl- 
edge to map the mass distribution in the Galaxy. If we 
could measure the distances, angular positions, line-of- 
sight velocities and proper motions of debris stars, we 
could integrate their orbits backwards in some assumed 
Galactic potential. Only in the correct potential will the 
path of the stream stars ever coincide in time, position 
and velocity with that of the satellite (see Fig. [TT|) . 



SIM measurements combined with ground-based line- 
of-sight velocities should provide everything needed to 
undertake the experiment; however, obtaining precision 
trigonometric parallaxes of numerous distant debris stars 
would involve a significant investment of SIM observing 
time. On the other hand, distances could be estimated ei- 
ther by using accurate photometric parallaxes (calibrated 
with SIM) for red giant/horizontal branch stars or by 
exploiting our expectations for the orbital energy distri- 
bution in the debris. In the latter case, we know the 
mean offset of the leading a nd trailing debr is from the 
satellite's own orbital energy (|Johnstonlll998f) and hence 
can solve the energy equation for the distances to stars 
in each of these groups by as suming each has this mean 
energy (| Johnston et al.lll999f ). This should be more ac- 
curate than using a photometric parallax so long as the 
offset in orbital energy is less than 5d d^/dr, where 5d 
is the uncertainty in distances and d^ j dr is the gradient 
in the assumed Galactic potential. 

If we find a coherent stream but not the associated 
satellite, the same techniques apply, but with the parent 
satellite's position and velocity as additional free parame- 
ters. The usefulness of Galactic tidal streams for probing 
the Milky Way potential has long been recognized, and 
the results of astrometric space missions to provide the 
last two dimensions of phase space inform ation for stream 
stars has been eagerly anticipated (e.g., iJohnston et al.l 
a99_9; Pcnarrubia ct al. 2006). 

Applying this method to simulated data observed with 
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the nas yr^^ precision proper-motions possible with SIM 
and kms~^ radial velocities suggests that 1% accuracies 
on Galactic parameters (such as the flattening of the po- 
tential and circular speed at the Solar Circle) can be 
achieved with tidal tail samples as small as 100 stars 
([Johnston et al.lll999l : lMaiewski et al. 200 6). Dynamical 
friction is not an important additional consideration if 
the change in the energy of the satellite's orbit in iVorb 
orbits is less than the range in the energies of debris par- 
ticles. For iVorb = 3, this condition is met for all satellites 
except the LMC/SMC and Sagittarius. Evolution of the 
Galactic potential does not affect the current positions 
of tidal debris, which respond adiabatically to changes 
in the potential and therefore yield direct information 
on the prese nt Galactic mass distribution independent of 
how it grew (|Pefiarrubia et al . 2006). 

Ideally, stars from several different tails at a variety of 
distances from the Galactic center and orientations with 
respect to the Galactic disk would be probed. It is also 
important to sample the tidal tails out to points where 
stars were torn from the satellite at least one radial orbit 
ago and hence, have experienc ed the full range of Galac- 
tic potential along the orbit (|JohnstoiJl2001h . Finally, 
the stars also need to have proper motions measured suf- 
ficiently accurately that the difference between their own 
and their parent satellite's orbits are detectable — this 
translates to requiring proper motions of order 100 /ias 
yr^^ for Sgr but 10 /xas yr~^ for satellites that are farther 
away. 

The possible existence of a significant fraction of the 
halo in the for m of dark satellites has been debate d 
in recent years (jMoore et al.lll999t iKlvpin et al.lll999f ). 
These putative dark subhalos could scatter stars in tidal 
tails, possibly compromising their use as large-scale po- 
tential probes, but astrometric measurements of stars 
in these tails could, on the other hand, be used to as- 
sess the importance of substruct ure flbata et al.l 120021 : 
iJohnston. Spergel fc Havdnl |2002[ ). Early tests of such 
scattering using only radial velocities of the Sgr stream 
suggest a Milky Wa y halo smoother than predicted 
(jMaiewski et al.l[2004l ). but this represents debris from a 
satellite with an already sizable intrinsic velocity disper- 
sion. Because scattering from subhalos should be most 
obvious on the narrowest, coldest tails (e.g., from glob- 
ular clusters) these could be used to probe the DM sub- 
structures, whereas the stars in tails of satellite galaxies 
such as Sgr, with larger dispersions initially and so less 
obviously affected, can still be used as global probes of 
the Galactic potential. 

In the last few years a number of well-defined Galac- 
tic tidal tails have been discovered at a wide range of 
radii and, with SIM, can be used to trace the Galac- 
tic mass distribution as far out as the virial radius with 
an unprecedented level of detail and accuracy. For 
example, M-giant stars associated wit h Sgr have now 
been traced entirely around the Galaxy (jMaiewski et al.l 
12003( 1 ■ the globular clusters Pal 5 and NGC 5466 both 
have tid al tails traced to ove r 20 d e grees from their 
centers (iGrillmair fc JohnsonI l2006t iBelokurov et al.l 
l2006cl : IGrillmair fc Dionatosll2006al) . two apparently cold 
streams over sixty degrees long have been discovered 
in the Sloa n Digit al Sky Survev (IGr illmair fc Dionato^ 
I2006bt IGrillmair 2006al : IBelokurov eral. .2006af ), and the 
nearby "Anticenter", or "Monoceros" stream has been 



shown t o be broken up into dynamically colder "tribu- 
taries" ()Grillmaiill2006bl ) . In the outer Galaxy, there are 
suggestions of debris associated wit h the Ursa Minor , 
Carina, Sculptor and Leo I dSphs, (iPalma et al.ll2003l: 
Maicwski et al. 2005; Westfall et all l2006l : iMufioz et all 
i200 6a: Sohn_et al. 2007,) , and evidence for o uter halo de- 
bris from the Sgr dSph (|Pakzad et al.l[200l ). With SIM, 
for the first time it will be possible to probe the full 
three-dimensional shape, density profile, and extent of 
(and substructure within) an individual DM halo. 

10.2. Hypervelocity Stars 

Hypervelocity stars (HVS) were postulated by iHillsl 
(1988), who showed that the disruption of a close bi- 
nary star system deep in the potential well of a massive 
black hole could eject one member of the binary at speeds 
exceeding 1000 km s~^. HVSs can also be produced by 
the interaction of a sin gle star with a binary black hole 
(|Yu fc Tremaind 12003). These remarkab le objects have 
now been discovered: iBrown et al.l (|2006D report on five 
stars with Galactocentric velocities between 550 and 720 
km s^^, and argue persuasively that these are HVSs in 
the sense that they are "unbound stars with an extreme 
velocity that can be explained only by dynamical ejec- 
tion associated with a massive black hole". It is likely 
that many more HVSs will be discovered in the next few 
years, both by ground-based surveys and by the Gaia 
mission. 

The acid test of whether these remarkable objects are 
HVSs is whether their proper motions are consistent with 
trajectories that lead back to the Galactic center. The 
magnitudes of the known HVSs range from 16 to 20, so 
their proper motions should be measurable by SIM with 
an accuracy of a few /j,as yr^^ . At the estimated distances 
of these stars (20 to 100 kpc) a velocity of 500 km s~^ 
corresponds to a proper motion of 1000 to 5000 yuas yr~^, 
so SIM should be able to determine the orientation of 
their velocity vectors to better than 1%. 

If these measurements confirm that the HVSs come 
from the Galactic center, then we can do more. For many 
stars, SIM can measure all six phase-space coordinates, 
but for HVSs, the orbits are far mor e tightly constrained 
because we know the point of origin. iGnedin et all (l2005h 
have pointed out that the non-spherical shape of the 
Galactic potential — due in part to the flattened disk and 
in part to the triaxial dark halo — will induce non-radial 
velocities in the HVSs of 5-10 km s~^, corresponding to 
10-100 ^as yr~^. Each HVS thus provides an indepen- 
dent constraint on the potential, as well as on the solar 
circular speed and distance to the Galactic center. 

10.3. Dark Matter Within Dwarf Galaxies 

Dwarf galaxies, and particularly dSph galaxies, are the 
most DM-dominated systems known to exist. Due to the 
small scale sizes (about 1 kpc) and large total mass-to- 
light ratios (approaching IOOMq/Lq), the dSphs provide 
the opportunity to study the structure of DM halos on 
the smallest scales. The internal structure of the dSphs, 
as well as their commonality within the Local Group, 
make possible a new approach to determining the phys- 
ical nature of DM with SIM. 

CDM particles have negligible velocity dispersion and 
very large central phase-space density, res ulting in cuspy 
density profiles over observable scales (|Navarro et al.l 
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Fig. 11. — A demonstration of the sensitivity of SIM to the Galactic potential using stellar streams. A Sagittarius-like tidal stream 
was created by the disruption of a dwarf satellite in a time-independent Galactic potential through a semi-analytical N-body simulation 
IILaw. Johnston &: Maiewskill2005l) . This resulted in the tidal stream demonstrated in panel (a). Complete 6-D phase space information on 
stars in the stream will be derived using the astrometric data from SIM. With such data in hand, guesses may be made on the strength and 
shape of the Galactic potential, and the orbits of the individual stars in the tidal streams run backwards under these assumed potentials. 
Panels (b) and (c) demonstrate what happens when the strength of the Galactic potential is underestimated by varying degrees: When 
the orbits are run backwards, the tidal stream stars orbit at too large a radius and do not converge to a common phase space position. 
In panels (e) and (f) the strength of the Galactic potential has been overestimated, and when the clock is run backwards the tidal stream 
stars assume orbits that are too small and once again do not converge on a common phase space position. In panel (d) a Galactic potential 
of the correct strength was guessed, and when the stream star orbits are run backwards, the tidal stream stars collect back into the core 
of the parent satellite. 

are well-fitted by King profiles, modulo sliglit variations. 
In tlie context of equilibrium models, the measured 
velocity profiles typically imply at least an order of 
magnitude more mass in DM than in stars, and imply 
mass-luminosity ratios that increas e with radius — in 
some cases quite substantially (e.g.. iKlevna et aLll2002D 
— though at large radii t idal e ffects may com plicate 
this picture (iKuhnl [l993l : iKroupa 1997: Munoz et all 
120051: IMufioz et al.ll2006al: ISohn et al.ll200l . Knowing 
whether mass follows light in dSphs or the luminous 
components lie within large extended halos is critical to 
establishing the regulatory mechanisms that inhibit the 
formation of galaxies in all subhalos (t ^lO.ip . 

Unfortunately, for equilibrium models the solutions 
to the Jeans equation are degenerate in that the dark 
matter density profiles are equally well-fitted by both 
cores or cusps. In particular, there is a strong degen- 
eracy between the inner slope of the DM density pro- 
file and the velocity anisotropy, /3, of the stellar or- 
bits; this leads to a strong dependency of the derived 
masses on (3. Radial velocities alone cannot break this 
degeneracy (Fig. [T^ . even if the present samples of ra- 
dial velocities are increased to sever al thousand stars 
(jStrigari. Bullock fc Kaplingha"tll2007f) . The problem is 
further compounded if we add triaxiality, Galactic sub- 
structure, and dSph orbital shapes to the allowable range 
of parameters. 



llQQTtlMoore et al.lll998f l. Warm Dark Matter (WDM), in 
contrast, has smaller central phase-space density, so that 
density profiles saturate to form constant central cores. 
Due to the small scale sizes of dSphs, if a core is a result 
of DM physics then the cores occupy a large fraction of 
the virial radii, which makes the cores in dSphs more ob- 
servationally accessible than those in any other galaxy 
type. Using dSph central velocity dispersions, earlier 
constraints on dSph cores have excluded ext remely warm 
DM, such as standard massive neutrinos (|Lin fc Fabeil 
119831: iGerhard fc Spergell [l99l ). More recent studies 
of the Fornax dSph provide strong co nstraints on the 
properties of sterile neutrino DM ( Goerdt et al.l l2006l : 
IStrigari et al][200611 . 

The past decade has seen substantial progress in 
measuring radial velocit ies for large numbers of stars in 
nearb y dSph galaxie s (lArrnand roff. Olszewski fc Prvo: 
19951: iTblstov et al I 120041: I Wilkinson et al. 200 ' 



Mufioz et all 120051 l2006bnWalker et all [200l 7 In aU 



dSphs, the projected radial velocity dispersion profiles 
are roughly flat as far out as they can be followed, with 
mean values between 7 and 12 km s~^. From these 
measurements, the DM density profiles are obtained 
by assuming d ynamical equilibrium and so lving the 
Jeans equation (jRichstone fc Tremaind [l986l ) . In this 
analysis, the surface density of the stellar distribution 
is required, and in all dSphs these stellar distributions 
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The only way to break the mass-anisotropy degener- 
acy is to measure more phase space coordinates per star. 
The Jeans equation written for the transverse velocity 
dispersion probes the anisotropy parameter differently 
than that for the radial velocity dispersion. Thus, com- 
bining proper motions with the present samples of radial 
velocities will provide orthogonal constraints and has the 
prospect to break the anisotropy-inner slope degeneracy. 

The most promising dSphs for this experiment will be 
the nearby (60-90 kpc away) systems Sculptor, Draco, 
Ursa Minor, Sextans and Bootes for which the upper gi- 
ant branches require proper motions of stars with V ~ 
19. The latter four of these d Sphs include t he most DM- 
dominated systems known (|Mateol Il998t iMunoz et al.1 
120061/) as well as a system with a more modest M/L 
(Sculptor). While Sagittarius is closer still, its strong 
interaction with the Galaxy and obvious tidal distor- 
tions indicate a system clearly not in dynamical equi- 
librium; thus it offers an interesting, possibly alternative 
case study for establishing the internal dynamical effects 
of tidal interaction. To sample the velocity dispersions 
properly will require proper motions of > 100 stars per 
galaxy with accuracies of 7 km s~^ or better (less than 
15 fxas yr~^). Detailed analysis shows that with about 
200 radial velocities and 200 transverse velocities of this 
precision, it will be possible to reduce the error on the 
log-slope of the dark matter density p rofile to about 0.1 
(|Strigari. Bullock fc KaDhnghatj|2007[ ). This is an order 
of magnitude smaller than the errors attainable from a 
sample of 1000 radial velocities, and sensitive enough to 
rule out nearly all WDM models (see FigurefT^. Obtain- 
ing the required transverse velocities, while well-beyond 
the capabilities of Gaia, is well-matched to the projected 
performance of SIM. 

10.4. Galactic Satellite Proper Motions 

Knowing the bulk proper motions of the dSphs is 
key to modeling not only their structure and evolu- 
tion, but whether they correspond to the predicted dark 
subhalos and how they relate to the two primary solu- 
tions to the missing satellites problem, i.e., the "very 
massive dwarf" versus the "very old dwarf" pictures 
(|Mashchenko. Sills fc Couchmanll2006l ). The derived or- 
bits of the dSphs can be compared directly to predictions 
for the orbits of infallin g DM sub structure (jGhigna et al.l 
[1998; Benson 2005; Ze ntner et al . 2005). These orbits 
can also be used to test for the reality of the purported 
dyna r nical families (e.g.. iLvnden-Bell fc Lvnden-Belll 
119951: iPalma. Maiewski fc JohnstonI I2002D hvpothesized 
to be the product of the break up of larger systems 
(jKunk el 1979; Lynden-Bell 1982), to deter mine the pos- 
sible connecti on of their appa rent alignment to local fila- 
ments (|Libeskind et al.ll2005), and to verify w hether dis- 
tant systems like Leo I fe.g. JSohn et al.ll2007l ) are bound 
to the Milky Way and can be used as test particles for 
measuring the mass profile of the Galaxy. 

The proper motions of the more distant satellites of 
the Milky Way are expected to be on the order of 100s 
of /xas yr~^, and to derive transverse velocities good to 
~ 10 km s~^ (~ 10%) requires a bulk proper motion 
accuracy of ~ 10/zas yr~^ for the most distant satellites 
(Leo I, Leo II, Canes Venatici), in which the brightest 
giant stars have V ~ 19.5. The constraints are slightly 
more relaxed for satellites at roughly 100 kpc distances. 
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Fig. 12. — A demonstration of the ability to recover infor- 
mation on the nature of dark matter using observations of dSph 
stars, from analytical modeling by Strigari, Bullock & Kaplingh^ 
II2007I ). Each panel shows the the 68% and 95% confidence regions 
for the errors in the measured dark halo density profile (log) slope 
(measured at twice the King core radius) and velocity anisotropy 
parameter /3 for a particular dSph; small crosses indicate the fidu- 
cial input model values. In the upper plots, derived from line of 
sight (RV) velocities for 1000 stars, the two parameters are highly 
degenerate, for both a cusp model (left panel), and core model 
(right panel). The addition of 200 proper motions from SIM pro- 
viding 5kms~^ precision transverse velocities (lower panels) dra- 
matically reduces the uncertainty in both parameters. 



~ 20 fi'AS yr~^ for systems with brightest stars at V ^ 
17.5. Such measurements for individual stars are well 
within the capabilities for SIM, which could derive the 
desired bulk motions for Galactic satellites with only a 
handful of stars per system (suitably placed to account 
for rotation and other internal motions, ^ 310.3^ . 

With much less per star precision, but many more 
stars, attempts to measure dSph galaxy proper 
motion s from the gr o und wi th long time b as elines 
(e.g.. iScholz fc IrwinI 11994 ISchweitzer et al] 119951 : 
ISchweitzer. Cudworth fc Maiewski Il997f ). with the 
Hubble Space Tel escope and short time baselines 
(jPiatek et al.l l2002l et seq.), an d with combinations 
of both ground and HST data (iDinescu et al.l I2004D 
still lead to bulk motions measurements with uncer- 
tainties of order the size of the motion of the dSph, 
and with significantly different results (e.g., two re- 
cent measures of the proper motion of the Fo rnax 
dS ph differ at the 2a level — IDinescu et al.l 120041 ver- 
sus iPiatek et al.l l2002f ) . Even for the closer Magellanic 
Clouds, there still remain significant a nd distur bing vari- 
ations in the derived proper niotions (iJoneset al. 1994; 
iKrou pa. Roser fc Bastianlll994l:lKroup a fc Bastian 1997' 



Ang uita. Lovola fc Pedreroi~ 20(301: ;M omanv fc Zaggi 
120051: iKallivavalil et all 120061 : iPedreros et all 12006^ 
These inconsistencies not only reflect the difficulty of 
beating down random uncertainties with ^/N statistics 
in systems with limited numbers of sufficiently bright 
stars — a problem that will be severe in the most re- 
cently discovered dSphs, where the popu lations of giant 
stars number at most in the t e ns (e.g.. IWillman et all 
I2005bl: iBelokurov etlHI l2006bl: iZucker et all I2006D — 
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but in dealing with numerous systematic problems, 
such as establishi ng absolute proper motion zero-point s 
(see discussion in iMaiewskilfigoa iDinescu et al.ll2Q04[ ). 
With its /^as-level extragalactic astrometric reference 
tie-in, and intrinsic per-star proper motion precision 
SIM will make definitive measures of the proper motions 
of the Galactic satellites that will overcome the previous 
complications faced by satellite proper motion studies. 

The same stellar precisions conferred by SIM to the 
study of distant dwarf satellites and star clusters fij llO.Sp 
may also be applied to proper motions of individual dis- 
tant field stars, whose full space motions can be derived 
and used as additional point mass probes of the Galac- 
tic potential, whether they are part of tidal streams or a 
well-mixed halo population. 

10.5. Globular Cluster Proper Motions 

Globular clusters, which have a spatial distribution 
that spans the dimensions of the Milky Way and which 
readily lend themselves to abundance and age assess- 
ments, have long served as a cornerstone stellar popu- 
lation for understanding galaxy evolution. Several dis- 
tinct populat ions of globular clust ers are known : a disk 
(jArmandrofl fl989l and/or bulge (jMinnitil [l996l ) popu- 
lation and at least two kind s of h alo clusters (jZinn I 
|1996[) . Since ISearle fc Zi^ (| 19781 ). the notion that 
halo globular clusters were formed in separate environ- 
ments — "protogalactic fragments" — later accreted 
by the Milky Way over an extended period of time 
has been a central thesis of Galactic structure stud- 
ies. More recently, direct substantiation of the hy- 
pothesis for at least some clusters has come from the 
identificatio n of clusters that are parts of the Sagittar- 
ius stream (llbata. Gilmore fc Irwinlll995t IDin escu et all 
120001: iMaiewski et al.ll2004l: iBellazzin i et al."2002. 2003) 
and the Monoceros structure (ICrane et al. 2003; 
iFrinchabov et al.l 120041 : IBellazzini et all l2004f ). These 
systems represent more obvious "dynamical families" 
associated with recent mergers having readily identifi- 
able debris streams, but presumably such mergers were 
even more common in the early Galaxy. Ancient merg- 
ers may have tenuous stellar streams today, and their 
identification will require 6-D phase space information 
to locate objects of common energy and angular mo- 
mentum. Tracing ancient mergers via their identified 
progeny will provide key insights into the evolution of 
substructure and star cl usters in hierarchical cosmologies 
(jPrieto fc Gnedinll200l . 

Of course, halo globular clusters will serve as valu- 
able test particles for determining the halo potential, but 
these cluster data will also play an essential role in un- 
derstanding clusters as stellar systems. The dynamical 
evolution of small stellar systems is largely determined 
by e xternal influences s uch as disk and bulge shocks 
(e.g.. iGnedin et al]ll999( l. so determining cluster orbits 
by measuring their proper motions will dramatically im- 
prove our understanding of their evolution and address 
the long-standing issue of whether the present popula- 
tion of these systems is the surviving remnant of a much 
larger initial population. The formation of globular clus- 
ters remains a mystery that can be better constrained by 
understanding cluster orbits. By obtaining definitive or- 
bital data for the entire Galactic globular cluster system, 
SIM will clarify the range of extant cluster/satellite dy- 



namical histories and how the Galactic ensemble evolved 
and depopulated. 

At present, < 25% of Galactic globular clusters have 
had any attempt at a measured proper motion, and reli- 
able data generally exist o nly for those clusters closest to 
the Sun (see summaries iniDin escu. Girar d fc van Altenal 
119991: iPalma. Maiewski fcJo hns ton 2 00^ As in the 
case of the Galactic satellites fii ll0.4p . even in the rare 
cases when appropriate data for proper motion mea- 
surements exist, analyses are hampered by critical sys- 
tematic errors, notably the tie-in to an inertial refer- 
ence frame. Galaxies yield unreliable centroids and 
QSOs have too low of a sky density at typical magni- 
tudes probed. Moreover, most of the outer halo glob- 
ular clusters, including some newly discovered examples 
(jCarraro2005; Will man et al 20Q5a) have very sparse gi- 
ant branches, reducing the effectiveness of averaging the 
motions of numerous members to obtain a precision bulk 
motion. SIM will immediately resolve these problems. 

11. ASTROMETRIC MICROLENSING 

What would an unbiased census of Galactic objects, 
dark and luminous, reveal? At a minimum, it would 
yield the frequency of black holes, neutron stars, and old 
brown dwarfs, which are either completely dark or so dim 
that they defy detection by normal methods. It might 
also find a significant component of the dark matter, al- 
though the majority of dark matter cannot be in the form 
of co mpact objects (lAlcock et al.l l2000t iTisserand et al.l 
|2007| ). The only known way to conduct such a census 
is to put a high-precision astrometry telescope in solar 
orbit. 

Masses of astronomical bodies can be measured only 
by the deflections they induce on other objects, typically 
planets and moons that orbit Solar System bodies and 
binary companions that orbit other stars. Masses of lu- 
minous isolated fleld stars can be estimated from their 
photometric and spectroscopic properties by calibrating 
these against similar objects in bound systems. Hence, 
photometric surveys yield a reasonably good mass census 
of luminous objects in the Galaxy. 

Dark objects like black holes are another matter. Mass 
measurements of isolated fleld black holes can be ob- 
tained only by their deflection of light from more distant 
luminous objects. Indeed, it is difficult to even detect 
isolated black holes by any other effect. However, to go 
from detection to mass measurement (and therefore pos- 
itive identification) of a black hole is quite challenging. 

Gravitational microlensing experiments currently de- 
tect about 500 microlensing events per year. The vast 
majority of the lenses are ordinary stars, whose gravity 
deflects (and so magnifles) the light of a more distant 
'source star'. As the source gets closer to and farther 
from the projected position of the le ns, its magnifica tion, 
A, waxes and wanes according to the lEinsteinI ()1936[ ) for- 
mula 

where u is the source-lens angular separation (normalized 
to the so-called Einstein radius 9e), to is the time of 
maximum magnification (when the separation is uq) and 
is the Einstein radius crossing time, i.e., is — ^e/M: 
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where is the lens-source relative proper motion. The 
mass M cannot be directly inferred from most events 
because the only measurable parameter that it enters is 
^E, and this is a degenerate combination of M, /z and the 
source-lens relative TTrci: 



rel; 



(5) 



where n = 4G/{c'^A\J) 8 mas Mq\ 

It follows immediately that to determine M, one must 
measure three parameters, of which only one (^e) is rou- 
tinely derived from microlensing events. Another such 
parameter is 9-^, which could be routinely measured from 
the image positions, if it were possible to resolve their 
O(mas) separation. A third is the "microlens parallax" 
tte — TTrci/^E- Hence the l ens mass can be extracted from 
9e and tte alone (see e.g.. lGouIdll2000[ ): 



M = 



^E 

ktte 



(6) 



Just as ^E is the Einstein radius projected onto the 
plane of the sky, tte is related to te = AU/tte, the Ein- 
stein radius projected onto the observer plane. And just 
as 9-E could, in principle, be measured by resolving the 
two images on the sky, tte could be routinely measured 
by simultaneously ob serving the ev e nt from two l ocations 
separated by O(fE) (lRefsda]|ll966l : lGouldlll995l) . "Rou- 
tine" measurement of both tte and 9^ is essential. As 
of today, there have been a few dozen measuremen ts of 
these parameters separately (e.g.. lPoindexteiil2005 f). but 
only one very exceptional microlensing event for which 
both were measured tog ether with sufficient precision t o 
obtain an accurate mass (IGould. Bennett fc Alvesll2005l) . 

In fact, such routine measurements are possible by 
placing an accurate astrometric and photometric tele- 
scope in solar orbit. For current microlensing experi- 
ments carried out against the dense star fields of the 
Galactic bulge, tTi-oI ~40 /xas, so for stellar masses, 9e ~ 
500 fias and te ^ 10 AU. Hence, a satellite in solar orbit 
would be an appreciable fraction of an Einstein radius 
from the Earth. As a result, the photometric event de- 
scribed by equation ([4]) would look substantially different 
than it would from the ground. From this difference, one 
could infer te (and so tte). 

Determining 9^ is more difficult. As mentioned above, 
this would be straightforward if one could resolve the 
separate images, but to carry this out routinely (i.e., 
for small as well as large values of 9^) would require 
larger baselines than are likely to be available in next- 
generation instruments. Rather, one must appeal to 
a more subtle effect, the deflection of the centroid of 
the two lensed images. This deflection is given by 



tne two lensea images. inis denection is given by 
(IMivamoto fc Yoshiilfl995l: IHog et al.|[l995l: IWalked[l995h 



A6' = 



(7) 



Simple differentiation shows that this achieves a maxi- 
mum at u = V2, for which M = 9^/^/^, roughly 1/3 
of an Einstein radius. Hence, if the interferometer can 
achieve an accuracy of 0(10 /ias) at the time when this 
deflection is the greatest, then 9e can be measured to a 
few percent. 



There are some subtleties as well as some challenges. 
Satellite measurements of "te are subject to a four-fold 
discrete degeneracy, which can only be reso lved by ap- 
pealing to higher-order effects (jGouldl [l995D . It is not 
enough to measure the centroid location to determine 
the astrometric deflection: one must also know the un- 
deflected position to which the measured position is to 
be compared, and this can only be found by extrapolat- 
ing back from late-time astrometry. And the precision 
of the mass measurement depends directly on the signal- 
to-noise ratio of the underlying photometric and astro- 
metric measurements. This is important because space- 
based astrometric telescopes are likely to be photon chal- 
lenged and so to require relatively bright (and hence rare) 
microle nsing events to provide accurate mass measure- 
ments. IGould fc SalimI (|1999( ) estimated that ~ 1200 
hours of SIM time would yield 5% mass measurements 
for about 200 microlenses. Most of these lenses will be 
stars, but at least few percent are likely to be black holes, 
and several times more are likely to be other dark or 
dim objects like neutron stars, old white dwarfs, and old 
brown dwarfs. Since such a census is completely new, it 
may also turn up unexpected objects. 



12. DYNAMICS OF GALAXY MOTIONS: NUMERICAL 
ACTION AND SIM 

If one could measure the proper motions of galaxies 
with global accuracies of a few /xasyr^^, one could ob- 
tain another two components of phase space with which 
to construct flow models and determine histories and 
masses for galaxies and galactic groups. For a galaxy 
1 Mpc away, 4/xasyr~^ corresponds to 19kms~^ trans- 
verse motion, which is small compared to the expected 
transverse motions in the field, ~ 100 kms""'^. With an 
instrument such as SIM, one could measure accurate po- 
sitions of a few dozen stars, with V < 20, in each galaxy 
and after ^--^5 years obtain proper motions with adequate 
accuracy. For a typical dwarf galaxy, after averaging ran- 
domly located stars, the contribution to the error from 
the internal motions would be only a few kms^^. For 
larger galaxies, simple rotation models, adjusted to the 
observed velocity profiles, can be removed from the mo- 
tions for < 20kms~"'^ accuracy. There are 27 galaxies 
known (all within 5 Mpc) that have stars sufficiently 
bright. 

Just beyond the outermost accessible 21-cm isophotes 
of galaxies, the dark matter distribution becomes un- 
known. Within the light-emitting parts of galaxies, rota- 
tion curves are flat and not falling according to Kepler's 
law, which implies that mass grows roughly linearly with 
radius. The total mass-to-light ratio depends critically 
on where this mass growth ends, but this is generally 
not observed. As a result, critical questions about dark 
matter on scales of galaxies to groups remain: Do dwarf 
galaxies have lower or higher mass-to-light ratios than 
regular galaxies? Do the dark matter halos of galaxies 
in groups merge into a common envelope? How do these 
mass components compare with the warmer dark mat- 
ter particles smoothly distributed across superclusters or 
larger scales. At present, we can only detect dark mat- 
ter through its gravitational effects; therefore, a careful 
study of the dynamics of nearby galaxies is one of the 
few ways to resolve these issues. 
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SIM measurements of the deviation from Hubble flow 
wiU be of lasting importance in the modeling of the for- 
mation of the Local Group, several nearby groups and of 
the plane of the Local Supercluster. Note that with Gaia 
only M31 and M33 have sufficient numbers of stars that 
are bright enough to attempt a proper motion determi- 
nation. If Gaia achieves 25 /xasyr"^ (the presently stated 
goal) or ~ lOOkms"^ at 1 Mpc, it will obtain only ~ 1-a 
detections for these two galaxies. 

12.1. Peculiar Velocities 

Most analyses of peculiar velocity flows have ap- 
plied linear perturbation theory (appropriate for scales 
large enough that overden sities are <C 1) to spher- 
ical in-fall (Peebles 1980') . Pec u har velocity anal- 
ysis (iShava . T uUv fc Pierce] UMl : iDekel et all fl99l 
iPike fc HudsonI 120051) have proven the general concept 
that the observed velocity fields of galaxies result from 
the summed gravitational accelerations of overdensities 
over the age of the universe. They also agree with virial 
analyses of clusters and WMAP observations that in- 
dicate the existence of a substantial dark matter com- 
ponent strewn roughly where the galaxies are. But 
these studies apply only to large-scales of >10 Mpc, 
the scales of superclusters and large voids. Spherical 
infall (including "timing analysis" and "turnaround ra- 
dius" ) studies do not presume low overdensities and have 
been applied to the sma ller scale of the Local Group 
()Lvnden-Bell fc Linlll977[ ). These studies indicate mass- 
to-light ratios for the Local Group of roughly M/L ~ 
IOOM0/L0, but the model is crude; non-radial motions 
and additional accelerations from tidal fields and sub- 
clumping are expected to be non-negligible and would 
substantially alter the deduced mass. A more complete 
treatment of solving for self-consistent complex orbits is 
required. 

12.2. The Numerical Action Method 

The application o f the Numerical Action Method 
(NAM, iPeeblesI Il989f ) allows one to solve for the tra- 
jectories that result in the present distribution of galax- 
ies (or more correctly, the centers of mass of the ma- 
terial that is presently in galaxies). By making use of 
the constraint that early time peculiar velocities were 
small, the problem becomes a boundary value differential 
equation with constraints at both early and late times. 
For each galaxy, a position on the sky, either a redshift 
or a distance and an assumed mass (usually taken to 
be the luminosity times an assumed mass-to-light ra- 
tio) are required inputs. In addition, one must presume 
an age of the universe (which is known from WMAP). 
Several studies comparing NAM with N-body solutions 
have shown that NAM can recover accurate orbits and 
position s or v elocities (iBranchini^ ^^'^^l fc Nusseil 120021: 
iPhelpsI [200I iRomano-Di'az et al.l l2005f ). iSharpe et all 
( 2001h have used NAM to predict distances from red- 
shifts and then compared these to Cepheid distance mea- 
sures. 

Figure [T3] (upper) shows the output of a recent NAM 
calculation for the orbits of nearby galaxies and groups 
going out to the distance of the Virgo Cluster. The orbits 
are in comoving coordinates. This is just a single solution 
of a set of several solutions using present 3-d positions as 
inputs. The four massive objects (Virgo Cluster, Coma 



Group, CenA Group, and M31) have been adjusted to 
provide best fit to observed redshifts. Figure [T3] (lower) 
is the same calculation, but now the coordinates are real- 
space rather than comoving, and the supergalactic plane 
is viewed edge-on. It appears that the plane of the Lo- 
cal Supercluster was not created by matter raining in 
from large distances but rather that material has never 
achieved great heights away from the plane. 

Assuming a constant mass-to-light ratio provides only 
a rough guess for the galaxy masses. Adding proper mo- 
tions measured by SIM will allow us to solve for the ac- 
tual masses of the dominant galaxies. The dwarf galax- 
ies are essentially massless test particles of the potential. 
Also, the two components of proper motion would be 
extremely useful for determining several supplementary 
but crucial parameters: the halo sizes and density falloff 
rates at large radii, the mass associated with groups aside 
from that in the individual galaxies, and the amount of 
matter distributed on scales larger than 5 Mpc. 

12.3. Preliminary Science Data 

Ongoing ground- and space-based observations should 
help in the analyses of nearby galaxy dynamics. Ac- 
curate relative distances of nearby galaxies, often with 
accuracy as good as 5%, are being made at a rapid rate 
using methods such as Cepheids, tip of the red giant 
branch, eclipsing binaries, maser distances, etc. SIM 
should help reduce distance errors by providing better 
and therefore more consistent calibrations of these tech- 
niques. In addition, measurements of the proper motions 
of a few galaxies (only M33 and IC 10 thus far) using 
masers supplement our knowledge. For the fortuitous 
cases in which both a maser and SIM measurements are 
made for the same galaxy, the maser information can 
be averaged to beat down errors from internal motions. 
However, it is unlikely that future maser proper motion 
measurements, if any, will be in the same galaxies as 
those to be measured by SIM and thus they may provide 
completely complementary information. 

12.4. Future NAM Studies using SIM Data 

It appears that the mass-to-light ratio for systems 
with early type galaxies is higher than systems with late 
type galaxies. Cluster mass-to-light ratios and ellipti- 
cal galaxy X-ray data tend to M/L ~ 300 Mq/Lq while 
group virial masses, turnaround radii, and the timing 
analysis lead to ratios of 50-100 in the field. In NAM 
calculations, we are beginning to see this trend as well, 
with Virgo and Cen A requiring ratios of ~ 450 and ordi- 
nary groups at ~ 120. A detailed study of the flow on the 
5 Mpc scale with all three components of velocity that 
SIM would allow, will be able to determine if the normal 
groups have additional dark matter at larger radii that 
could bring the baryon-to-dark matter ratios of the two 
species into better alignment. 

The total masses of galaxies and the sizes of their ha- 
los are essential parameters for an understanding of dark 
matter and large-scale structure formation. It may be 
that the mass of the dark matter particle is in the range 
in which it is neither completely cold nor hot. If this is 
so, the clumping scale of the dark matter will provide 
a mass estimate of the dark matter particle. Or it may 
be that the dark matter particle is mildly dissipative. 
In which case, the halos will be flattened and the orbits 
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Fig. 13. — Upper - The trajectories of nearby galaxies and groups 
going out to the distance of the Virgo Cluster from a Numerical 
Action Method calculation with parameters M/L = 90 for spirals 
and 155 for ellipticals, Qm = 0.24, Q/^ = 0.76. The axes are in the 
supergalactic plane (SGX-SGY) in comoving coordinates. There 
are 21 time steps going from z = 40 to now. The large circle is 
placed at the present position and the radius is proportional to the 
square root of the mass. The present estimated distances were fixed 
and the present redshifts were unconstrained. Lower - A view of 
the galaxy trajectories, except that the coordinates here are real- 
space (proper) rather than comoving, and the axes are SGY-SGZ 
allowing a view of the collapse into a plane. 



of dwarf galaxies in different planes about some massive 
galaxy would be subjected to different effective masses. 
These would be important clues to revealing the physics 
and identity of the dark matter particle. If we combine 



the constraints on the 1 — 5 Mpc scale with flow studies 
on larger scales, it will be possible to strongly constrain 
a constant density component, resulting in a new astro- 
physical limit to the mass of the neutrino. 

13. QUASAR JETS AND ACCRETION DISKS 

Galaxies that possess spheroidal bulges also appear to 
possess central supermassive black holes whose mass is of 
order 0.1% of the bulge. Episodic accretion of gas, dust 
and stars onto the central black hole gives rise to the 
active galactic nucleus phenomenon. The potential SIM 
discovery space for AGN observations is broad because 
direct measurement on scales less than ~ 100 /xas has 
previously been done only at radio wavelengths. Most 
of what is known about the optically-emitting inner re- 
gions of AGN comes from optical variability studies of 
the Broad Line Region: the intrinsic size-scales are very 
small (light days to a light year, or ~ 0.1 — 10 /xas) in ex- 
tent. The ability of SIM to measure motions and position 
differences on microarcsecond scales means that we can 
study AGN on scales of tens to hundreds of Schwarzschild 
radii, perfect for studying accretion disks, jet coUimation 
and possible orbital motion of binary black holes. 

Some key questions that SIM can answer are: (1) Does 
the most compact optical emission from an AGN come 
from an accretion disk or from a rclativistic jet? Does 
this depend on whether the AGN is radio loud or radio 
quiet? (2) Do the cores of galaxies harbor binary su- 
permassive black holes remaining from galaxy mergers? 
(3) Does the separation of the radio core and the optical 
photocenter of the quasars used for the reference frame 
tie change on the timescale of the photometric variabil- 
ity, or is the separation stable? The use of quasars in 
the astrometric grid to approximate a perfectly inertial 
frame is described in Appendix [Cl We first briefly review 
the basic properties of AGN and the physics of the cur- 
rent accretion disk-jet paradigm, and then describe how 
SIM will answer these questions, including examples of 
what SIM should see for specific targets. 

13.1. Basic Quasar Properties: Radio Quiet and Radio 
Loud Sources 

The observational signature of a quasar is an optically- 
emitting source that is very bright intrinsically (up to 
10^'^ Lq) but very small in physical size (~ 10^^ cm). 
At a 1 Gpc distance, for example, the size of the cen- 
tral quasar engine is of order 0.1 /las and that of the 
broad line region of order 1 ^as. About 90% of all quasars 
are 'radio quiet'; their emission is dominated by optical 
and X-ray emission. The optical continuum spectrum 
consists of a fairly steep power-law, sometimes with a 
"Big Blue Bump" in the blue or near ult raviolet region 
(jBand fc Malkanlll989l: [Zheng et al.lll995f) . In addition, 
there are broad emission lines from highly-ionized ele- 
ments, which are produced rather close to the central 
source 10^^^^^ cm) and ionized by it. 

The remaining 10% of quasars are 'radio loud'; they 
have strong diffuse radio emission in addition to all the 
properties of radio quiet quasars, with radio jets extend- 
ing on both sides of the optical source out to the external 
radio lobes. In many cases these jets flow at relativistic 
speeds (Lorentz factors of 10 or more). About 10% of 
these (i.e., about 1% of all quasars) are blazars — strong 
and variable compact radio sources, which also emit in 
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the optical (mainly red and near infrared), and in X-rays 
and 7-rays. They are believed to be normal radio-loud 
quasars viewed by us nearly end-on to the jet. Rela- 
tivistic beaming toward the observer can produce an en- 
hancement of an order of magnitude or more in apparent 
radio luminosity, as well as apparent proper motions of 
jet components of up to 1000 //asyr~^. If they display 
similar internal proper motions in their optical jets, these 
would be readily detectable with SIM. 

13.2. The Location of the Most Compact Optical 
Emission in the Accretion Disk-Jet Paradigm 

In order to understand how SIM observations can pro- 
vide insights into the physical processes in AGN, we need 
to briefly review the major components and current con- 
straints on their parameters. A sketch of the canonical 
AGN model is presented in Figure [HI Observations on 
this size scale are indirect, so this picture is based largely 
on theoretical models. There are three possible sources 
of compact optical emission: the accretion disk, the disk 
corona or wind, and the relativistic jet. For a number 
of nearby AGN, SIM will be able to distinguish which 
component dominates and to study AGN properties on 
scales of 1 pc or less. 

13.2.1. The Big Blue Bump and Hot Corona 

In radio quiet quasars, there should be two sources of 
optical-UV emission: 1) thermal accretion disk emission, 
which produces the 'big blue bump', and which dom- 
inates in the blue; and 2) a nonthermal corona whose 
origin is a steep power-law ionizing source, which domi- 
nates in the red. Both of these emission regions should 
be physically centered on the black hole within ~ 1/ias, 
and both should have a similar ~1 /xas size. 

In the high-accretion case, the disk produces a ther- 
mal peak in the near ultraviolet region. For a typical 
10^ Mq black hole system, accreting at 10% of Msdci, 
the diameter of the portion of the disk that is radiating 
at a temper ature of K or above is ~ 3.6 x 10^^ cm, 
or 0.012 pc (jShakura fc Sunvaevl[l973l ). At the distance 
of M87 (16 Mpc) this region would subtend an angular 
size of ~ 160 fia.s, while at moderate redshift {z = 0.5) 
the angular size would be only 2 fias. 

The hot corona is believed to be responsible for 
exciting the observed emission lines of the high- 
ionization spe cies of car bon, silicon, nitrogen and oxy- 
gen, ( Ostcrbr ock &: Mat thews 1986). Coronal emission, 
therefore, is likely non-thermal — either optical syn- 
chrotron or inverse-Compton-scattered emission from a 
radio synchrotron source. M odels of this ionizing source 
(e.g.. 'Band fc Mal kanlll989l ) indicate a size of only ~ 70 
Schwarzschild radii. At moderate redshift (e.g., 3C345, 
z = 0.6) this subtends an angular size of only ~ 1 ^as, 
centered on the black hole and comparable in size to the 
Big Blue Bump. 

13.2.2. The Beamed Relativistic Jet in Compact Radio-Loud 

Quasars 

Radio-loud quasars possess powerful radio-emitting 
jets, in addition to the central optical-UV source seen 
in radio-quiet quasars. These jets are often relativistic 
flows, ejected in a direction approximately perpendicu- 
lar to the accretion disk (see Figure fT4|) by the effects of 



magn etic fields ()Meier. Koide fc Uchi"dall2001l : [McKinnevl 
l2006f ). If an AGN is radio- loud, and viewed at a small 
angle to the jet axis (a 'blazar'), then a third source of 
optical continuum emission may be seen, in addition to 
the big blue bump and hot disk corona: emission from 
the relativistic jet beamed toward the observer. This 
emission would be readily detectable by SIM, because it 
would be offset with respect to the other sour ces of emis- 
sion, and highly variable. The iKonigll (jlQBlf ) model for 
relativistic jets predicts that the majority of the optical 
emission comes not from synchrotron emission from the 
base of the jet but from synchrotron-self-Compton emis- 
sion in the region of the jet where the synchrotron emis- 
sion peaks in the radio or millimeter (iHutter fc MufsonI 
Il986f ). A detailed application of this model to 3C 345, 
for example, predicts the optical emission to be offset 
~ 80 /ias from the center of mass of the system and 
nearly coincident with the 22 GHz radio e mission which 
lies ^ 70/xas or 0.4 pc from the black hole (jUnwin et al.l 
I1994D . Not only is this readily detectable with SIM, but 
the vector position of the variability on the sky can be 
compared with the jet orientation seen on larger scales. 

13.3. Distinguishing AGN Models using 
Color- Dependent Astrometry 

SIM can directly test the above modes by measuring 
a color-dependent shift in the astrometric position. A 
displacement in optical photocenter between the red and 
blue ends of the passband can be readily measured, and 
is very insensitive to systematic errors; measurements 
will be primarily limited by photon noise not instrument 
errors. This color-dependent displacement, and its time- 
derivative in variable sources, is a vector quantity on the 
sky whose alignment can be compared to, say, the orien- 
tation of a radio jet imaged by VLBI or the VLA. 

In radio-quiet quasars we do not expect to see a color 
shift, because of the absence of any contribution of a 
relativistic jet whose optical emission might introduce 
an astrometric asymmetry. The red emission from the 
corona and the blue emission from the disk both should 
be coincident with the central black hole within ^ 1 /ias. 
Any color-dependent astrometric shift seen in radio-quiet 
quasars would challenge the current models of accreting 
systems in AGN. 

By contrast, the astrometric position may be strongly 
color-dependent in any object with strong optical jet 
emission. So while the blue end of the spectrum should 
be dominated by the thermal disk, the red region may be 
dominated by the beamed relativistic jet. The relative 
contributio ns may vary as the activity level changes (see 
iFerrero et al. 2006). Furthermore, we would expect any 
variability to be aligned on the sky with the direction of 
the color shift itself. 

An example of a moderate-redshift jet-dominated 
quasar is 3C 345 (z = 0.6) for which the red optical jet 
emission should be roughly coincident with the 22 GHz 
radio emission, at about 80 /ias from the black hole. For 
3C273 {z = 0.16), the separation may be as large as 
300 /ias. Not only is such a large shift readily detectable, 
but SIM could also detect variations in the offset with 
time. 

In the nearby radio galaxy M 87 we expect the red 
optical emission should be dominated by the accretion 
disk corona because its jets are not pointing within a few 
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Fig. 14. — Schematic diagram of the structure of a typical quasar on scales from 0.002 — 100 pc, after [Elvli II2000I) and Niall Smith 
(private communication, Cork Institute of Technology). Upper — On scales from ~ 0.1 — 100 pc, the main components are the broad- 
and narrow-line emission clouds and a dust torus that obscures the black hole in certain orientations. In radio-loud quasars, powerful (and 
often relativistic) jets are ejected from the central engine (Konigl 1981). These jets emit over the entire EM spectrum, with optical/IR, 
dominating on scales < 0.1 pc, and radio emission dominating on larger scales. Lower — On scales < 1 pc, a geometrically thin disk 
forms around the mass ive accreting black hole. In some quasars, the innermost portion of the accretion disk is seen as the 'big blue bump' 
flBand fc Malkarilll989l) . All quasars also have a source of non-thermal emission, a corona, which dominates in the red, and produces 
enough UV radiation to ionize the broad and even narrow line regions. In radio-loud quasars, the base of jet also will produce non-thermal, 
typically highly-variable, red emission that can outshine the central corona. SIM can probe the sub-parsec structure using time-dependent 
and color-dependent astrometry (see ii ll3.3l l. The blue disk and red corona are very compact and spatially coincident, so SIM would not 
expect to detect time- or color-dependent astrometric shifts in these components. However, if the jet dominates, then a color shift of 
~ 0.5 pc (for a 2 = 0.6 quasar; see text) would be seen, and any variability would be along the direction of the jet. SIM will be able to test 
these predictions for radio-quiet and radio-loud quasars. 



degrees of our line of sight. Therefore, SIM should not 
see a significant color shift in this source. However, M 87 
is so close that we might see an absolute position shift 



between the measured radio photocenter and the overall 
optical photocenter. This shift should be even larger for 
this low-rcdshift radio galaxy than for 3C 345 — perhaps 
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in the 1000 /^as range. 

13.4. Finding Binary Black Holes 

Do the cores of galaxies harbor binary supermassive 
black holes remaining from galaxy mergers? This is a 
question of central importance to understanding the on- 
set and evolution of non-thermal activity in galactic nu- 
clei. SIM can detect binary black holes in a manner 
analogous to planet detection: by measuring positional 
changes in the quasar optical photocenters due to orbital 
motion. If a quasar photocenter traces an elliptical path 
on the sky, then it harbors a binary black hole; if the mo- 
tion is random, or not detectable, then the quasar shows 
no evidence of binarity. If massive binary black holes are 
found, we will have a new means of directly measuring 
their masses and estimating the coalescence lifetimes of 
the binaries. 

An AGN black hole system (see Fig. [T5|) can occur near 
the end of a galactic merger, when the two galactic nuclei 
themselves merge. Time scales for the nuclei themselves 
to merge, and the black holes to form a binary of ~ 1 pc 
in size, are fairly short (on the order of several million 
years) and significantly shorter than the galaxy merger 
time (a few hundred million years). Furthermore, once 
the separation of the binary becomes smaller than 0.01 
pc, gravitational radiation also will cause the bi nary to 
coalesce in only a few million years (|Kroliklll999| ). How- 
ever, the duration of the 'hard' binary phase (separation 
of 0.01 - 1 pc) is largely unknown, and depends critically 
on how much mass the binary can eject from the nucleus 
as it i nteracts with ambient gas, stars, and o ther black 
holes (| Yul 120021 : iMerritt fc Milosavlievidl2005f l. Depend- 
ing on what processes are at work, the lifetime in this 
stage can be longer than the age of the universe, implying 
that binary black holes are numerous - or as short as the 
"Salpeter" time scale ~ {Mi+M2)/MEdd ^ 5 x lO^yr, 
implying that binaries might be rare. Therefore, the 
search for binary black holes in the nuclei of galaxies will 
yield important information on their overall lifetime and 
on the processes occurring in galaxies that affect black 
holes and quasars. 

A promising candidate would be an object that looks 
like a quasar (with broad lines and optical-UV contin- 
uum) but with absolute luminosity somewhat less than 
10 % of the Eddington limit expected from the central 
black hole. These might be objects with a large, but 
dark, primary central black hole that is orbited by a sec- 
ondary black hole of smaller mass; the secondary would 
have cleared out the larger hole's accretion disk inte- 
rior, but still will be accreting prodigiou sly from the in- 
ner e dge of the primary's disk (Milosavl ievic fc PhinnevI 
[20051) . In this case, the astrometric motion would indi- 
cate the full orbit of the secondary about the primary, 
which could be a few to a few hundred /ias, depending 
on the source distance. Figure [15] is a simplified sketch 
of this situation. OJ 287 {z = 0.3) is a candidate , based 
on brightness variations wi th 12-year periodicity (jKidgeil 
120001 : [Valtonen et al.ll2006l ). and we estimate that about 
14 yitas of orbital motion may be expected during a five- 
year span. 

14. COSMOLOGY WITH SIM 

In this section, we describe some contributions that 
SIM can make to the fields of cosmology and dark en- 
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Fig. 15. — Schematic diagram of a binary black hole system. 
Binary separations of order 0.01 — 1 pc may be expected. While 
the larger object will have a large disk, the smaller will clear out 
the center of the disk, rendering the larger hole quiet, and leaving 
only the lower-mass black hole as a source of ionizing radiation and 
jets, as in Fig. 1141 The orbital shift of the secondary hole will be 
readily detectable by SIM especially for nearby quasars {z < 0.2). 
Expected orbital periods for a 10^ Mq primary range from ^ 3 yr 
for a 0.01 pc orbit to ~ 90 yr for a 0.1 pc orbit. 

ergy (DE). There are at least three distinct types of con- 
tributions: 1) the measurement of accurate stellar ages 
which provide a robust lower bound for the age of the 
Universe (to), 2) the estimation of the primordial Helium 
abundance, and 3) the determination of percent-level dis- 
tances to several galaxies in the Local Group. Stellar 
ages and helium abundances are best determined from 
spectroph otometry and accurate parallaxes of eclipsing 
binaries (■Lestradell2000l ). For the present purpose, we 
will assume that the Local-Group distances are equiv- 
alent to determining the Hubble constant (Hq) itself. 
In fact, a robust maximum value for the equation of 
state of dark energy derives fr om the product to and Ho 
(jBothun. Hsu fc Murravll2006) . 

An error level for Hq of one to two percent is opti- 
mal for determining the EOS of Dark Energy and the 
total energy density of the Universe. Errors of the or- 
der of 1% and 0.1% for w and iltot, respectively, are 
achievable with Planck-like CMB data and a 2% error 
on Hq. Assuming ^/N statistics, 2% distance errors may 
be achieved with SIM observations of 43 stars in M31, 
200 stars in M33 and 2,025 stars in the LMC. Local- 
Group distances do not determine Ho, rather they allow 
for a re-calibration of secondary distance indicators that 
extend beyond the Local Group such as Cepheids and 
the TRGB (see ©. The disks of these galaxies encom- 
pass a wide range in metallicity, extinction and number 
density, and thus provide ample opportunity for an accu- 
rate calibration of methods such as the TRGB, Cepheids, 
RR Lyrae, surface-brightness fiuctuations and so forth. 

14.1. The Extragalactic Distance Scale 

When assuming a fiat ACDM (cold Dark Matter) 
model, the fluctuations of the cosmic microwave back- 
ground (CMB) as observed by WMAP imply a value of 
Ho = 73 ± 3kms"^Mpc-^ (jSDergel et al.l [2001 [200l : 
hereafter referred to as WMAP07). This value is very 
close to the HST-derived value of Ho = 74 ± 2 (ran- 
dom) ± 7 (systematic) kms~"'^Mpc~^ (Freed man et al.l 
l2001h . However, if the flatness assumption is abandoned, 
WMAP by itself hardly constrains Ho because WMAP 
measures the product of the normalized matter density 
(fi™) and h'^ (WMAP07, their Figure 20), as well as the 
baryon density (Slfcft.^). Here, h = lOOkms""^ Mpc~^/i7o 
and flmh^ = ojm — 0.126 ± 0.009. Thus, an independent 
and accurate determination of Ho would determine flm 
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and fib, and, if we assume a flat universe, the dark-energy 
density {flnE)- 

While an accurate determination of Hq has many ad- 
vantages, here we will concentrate on those with direct 
cosmological implications: 1) determining the equation 
of state (EOS) defined as the ratio of pressure to density 
{w = p/p) of dark energy (WMAP07), and 2) determin- 
ing the total density (fltot) of the universe. The equation 
of state of dark energy tells us something about its na- 
ture: the cosmological constant, strings, domain walls 
and so fort h predict different valu es for the equation of 
state fe.g-. lPeebles fc Ratra|[200l . 

According to Hu (2005, hereafter H05), a change in 
the EOS of Dark Energy by 30% leads to a change in Hq 
by 15% if a constant w is assumed, and by about one- 
half that much if w changes with redshift (H05, Figures 
3a and 3b). Thus, an accuracy of 1% in Hq leads to a 
percent-level determination of the EOS of dark energy. 
However, the absolute value of w and its rate of change 
can not be independently determined by fixing Hq. In 
fact, once Hq is fixed locally, the determination of w 
and/or its slope and/or the curvature term requires sub- 
percent level determination of Hq at z ~ 0.1 — 1.5 (H05, 
Figures 3c and 3d). 

As discussed above, if the CMB parameters are in- 
finitely well known, we expect that uncertainties in the 
EOS of Dark Energy remain at the level of the uncer- 
tainty in Hq. Below, we estimate analytically the effects 
of increasing both the accuracy of the CMB data and of 
Hq. The results of this exercise are consistent with the 
II05 findings. Because SIM is a targeted mission, a SIM- 
based determination of Hq can obtain, more or less, any 
required accuracy if an appropriate amount of observing 
time is allocated. Thus, our analytical results are useful 
in quantifying more precisely the resulting accuracy in 
the EOS of dark energy given a certain expenditure of 
SIM time, and vice versa. 

WMAP and other data currently constrain the EOS: 
w - (-0.826±0.109) - (0.557±0.058) w„ h-^ - -0.95± 
0.11, where we follow WMAP07 and assume a con- 
stant w , but allow for a non-zero curvature term 
[20073). This follows from various relations between the 
vacuum energy (JIa), the matter density, the spatial 
curvature of the universe {Hk) and the Hubble con- 
stant. Likewise, current data yield: Qtot = + ^ 
(0.9438 ± 0.0114) + 0.225 ujm h'^ - 0.996 ± 0.016. Thus, 
current data allow for the determination of the total den- 
sity of the universe to plus or minus one percent, while 
the EOS of dark energy is known to about 10%. 

To eliminate the uncertainty associated with h, one 
would want to determine the Hubble constant via 
trigonometric parallaxes of nearby galaxies. Unfortu- 
nately, this is not possible with foreseeable/planned 
technology. However, the 'Rotational Parallax' method 
(see S lTl?^ should be almost a s good (jPeterson fc Shaol 
Il997t lolling fc PetersonI l2000l hereafter referred to as 
OP2000).^ 

OUina (|2007af ) estimates the effects of more accurate 
CMB data and a better Hubble constant on the EOS of 
Dark Energy and finds that, even at Planck accuracy, 
the errors on w and ^Itot are only slightly smaller than 
the current values. The accuracy of the EOS of Dark 
Energy only improves significantly when the accuracy of 
Hq is improved. 



The results are summarized in Figure [TBI which shows 
the attainable accuracy on w as a function of the im- 
provement of our knowledge of the CMB, with respect 
to the WMAP 3-yr data, for four values of Hq accuracy. 
As expected, the error on w (e^) decreases as knowledge 
of the CMB parameters improves. However, these curves 
show that e„ approaches a limiting value. The curve 
that corresponds approximately to the current error on 
Hq {cho ~11%, top curve) indicates that an error of only 
about 9% can be obtained for w, even with Planck-like 
CMB accuracies (jEfstathiou et al.ll2005[ ). Note that this 
accuracy is only slightly better than our current knowl- 
edge as set by the WMAP 3-yr data. In order to signifi- 
cantly improve our knowledge of w, we need Planck data 
and also need Hq to much better accuracy. Achieving 
€w = 2.3%, requires an accuracy of 1.1 % on iJo. Errors 
on fltot behave much like those on w, but with roughly 
ten times smaller amplitude. 
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Fig. 16. — The accuracy of the equation of state (EOS) of Dark 
Energy ui as a function of the uncertainty in the CMB data, essen- 
tially the factor by which the error on u!m decreases, for four values 
of accuracy of the Hubble Constant (curves labeled in % accuracy 
on Ho). Vertical lines (from left to right) indicate the accuracies in 
the CMB parameters using WMAP 1-yr, 3-yr, and 8-yr data, and 
expected Planck data (Efstathiou et al. 2005) normalized to the 
3-year WMAP data. The resulting w errors are 8.9%, 4.8%, 3.0% 
and 2.3% for CMB accuracies at the level expected from Planck. 
A significant improvement in w requires both Planck CMB data 
and a more accurate Hubble Constant. 

Further improvements in our knowledge of the Dark 
Energy EOS may come from improvements in observa- 
tions other than Hq, though reducing the uncertainty 
in Hq makes the biggest difference. The other data 
sets used by WMAP07 are: large-scale structure ob- 
servations, galaxy redshift surveys, distant type-la Su- 
pernovae, Big-Bang nucleosynthesis, Sunyaev-Zel'dovich 
fluctuations, Lyman-a forest, and gravitational len sing. 
The Dark Energy Task Force (jAlbrecht et al.ll2006f ) con- 
siders how future instruments and space missions might 
significantly reduce the error in e^, through improvements 
in these datasets in four stages. It defines Stage I to be 
the current state of the art, and Stage IV to comprise a 
Large Survey Telescope (LST), and/or the Square Kilo- 
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meter Array (SKA), and/or a Joint Dark Energy Mission 
(JDEM). Each higher stage represents an improvement 
of by a factor of about three with respect to the previ- 
ous stage i Albrecht et al.|[200^ . It concludes that when 
these data are of Stage IV quahty, a reduction in the 
error on Hq by a factor of two matters at most 5 0% for 
the acc uracy of the Dark Energy EOS. However. lOUingj 
(l2007aD argues that the effects of smaller errors on Hq 
are especially important at the early stages, estimating 
that a decrease of by a factor of ten yields an im- 
provement of by factors of: 3.9, 3.0, 2.4 and 1.6 for 
Stages I, II, III and IV, respectively. 

14.2. Rotational Parallax Distances 

The method of Rotational Parallaxes (RP) combines 
proper motions (fi) and radial velocities (Vr) of stars in 
external galaxies to yield bias-free single-step distances 
with attainable accuracies of the order of one percent. 
This r nethod is analogous to the orbital parallax tech- 
nique (|Armstrong et al.lll992D . 

For a nearby spiral galaxy at distance D (in Mpc) 
which is inclined by i degrees and with a rotation 
speed of Vckms^^, the proper motion due to rotation 
is Vc/{kD) ^aayT"^, with k ~ 4.74 kms^V(AU yr^^). 
For M33, M31, and the LMC we find: ^i -24, -74 and 
~192 /iasyr"^, respectively, which are easily measured 
by SIM. A simplified RP method uses stars on the major 
and minor axes at similar galactocentric distances of a 
galaxy with a flat rotation curve. The minor-axis proper 
motion measures the circular velocity divided by I?, while 
the ratio of minor to major-axis proper motions is sim- 
ply cosi (jPeterson fc Shaolll997f l. However, stars need 
to be close to the principal axes, making it difficult to 
find enough targets. Generalizing this method we find: 



K y X^x + V'l^y' ' 

with X and y' the major-axis and minor-axis p osition 
of indi vidual stars in the target galaxy (OP2000, OUind 
[2007^. The attainable distance errors per star are esti- 
mated to be 13%, 28% and 90% for M31, M33 and the 
LMC, respectively. 

14.3. Realistic Rotational Parallaxes 

The rotational parallax method allows distances to be 
made to M 31, M 33 and the LMC in a single step, to ac- 
curacies of about 0.92, 2.0 and 6.4% respectively. Achiev- 
ing these bias-free single-step distances requires care- 
ful modeling of non-circular motions (due to spiral-arm 
streaming motions, perturbations from nearby galaxies, 
a bar, and warps) which could otherwise bias the dis- 
tance determination (OP2000). These effects are likely 
to be significant for the LMC. In order to achieve er- 
rors of several percent, it will be necessary to correct 
for any sizable deviations from circular motion. OP2000 
show that a correction can be achieved with SIM-quality 
proper motions. For a disk geometry, the combination 
of SIM proper motions and ground-based radial veloci- 
ties will yield four of the six phase space coordinates for 
individual stars, where the missing components can be 
chosen to be, for example, the vertical displacement (z) 



of the star with respect to the galaxy plane and the av- 
erage z velocity (|Onindl2007a[) . For M31 and M 33, this 
assumption is likely to be reasonable. Furthermore, we 
can make different assumptions regarding the fifth and 
sixth phase-space parameters and require consistency be- 
tween the results. Thus, the rotational parallax method 
will yield very reliable distances. This contrasts other 
techniques that claim to yield extra-galactic distances at 
the few-percent level such as eclipsing binaries, Cepheids, 
nuclear water masers and so forth. In fact, these other 
techniques rely on additional assumption and/or inaccu- 
rate slopes and/or zero points, albeit that many of those 
problems are likely to be c alibrated by SIM and Gaia 
(see review m lOlling|[2007ah . 

Because non-circular motions can be correlated on 
large scales, a substantial number of stars need to be 
observed (spread out over a large area of the galaxy) to 
be able to apply reliable corrections. OP2000 estimate 
that at least 200 stars are required to achieve the 1% 
distance error for M31 noted above. 



14.4. Other Local Group Distances 

SIM can be used to apply the "orbital parallax" tech- 
nique to binary stars in other nearby Local Group dwarf 
galaxies. This would increase the number and the range 
of types of galaxies that can be used to calibrate other 
rungs of the distance ladder. Let us estimate the required 
accuracies by neglecting inclination effects and eccentric- 
ity. Then, the orbital velocity, w, the semi-major axis, a, 
and the period, P, for each of the components yield the 
distance: 

2?ioofepc = 10 « .^^^I^ ^ 298 -^1^ (9) 

where Vkm/s = K2TTaAu/P and uau = a^i,s/DMpc- 
Because the distance error (Ad) scales according to 
(Ad/D)^ = {AJaf + {Ap/Pf + {A^/vf, a 1% dis- 
tance error requires that the errors on semi-major axis 
(Aa/a), period (Ap/P), and radial velocity (A„) are all 
at the sub-percent level. 

Given sufficient observing time, the period and the or- 
bital velocity and their errors can be determined from the 
ground with the required accuracy. Also, short period bi- 
naries are unlikely to have survived the expansion of the 
primary, while long-period orbits will be poorly sampled 
during the SIM mission and ground-based observations, 
so we will assume periods between 2 and 10 years. We 
simulate a population of binaries at 100 kpc with a 1 Mq 
primary. Secondaries are drawn from the stellar initial 
mass function and are in circular orbits with a period dis- 
tribution according to iDuquennov fc Mavod (|T991). We 
assume a 100% binarity rate. The result is that 1% of 
stars lie in the required period range (median is 7 yr) if 
a^ias > 12.5. For these systems, the median projected 
orbital velocity is Vo;km/s ~ 1-9. To achieve the 1% dis- 
tance accuracy goal, SIM would need (100/12.5)^ — 64 
observations per star, while the radial-velocity program 
would need to reach an accuracy of 19 m s~^. At 138 
(200) kpc, only 0.4% (0.042%) of binaries satisfy these 
criteria and have: Py^. — 7.9 (9.3), a^as ^ 14.9 (13.2) 
and Vo.km/s ^ 1-9 (2.5), where 138 kpc corresponds to 
the Fornax distance. 
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Excluding the Sagittarius dwarf and the Magellanic 
Clouds , we id entify five galaxies within 100 kpc in the 
iMateol ()1998[ ) compilation of local group dwarfs, and 
seven within 200 kpc. If we make the optimistic assump- 
tion that the total luminosity of these systems comes only 
from stars 1 mag below the Tip of the Red Giant Branch 
(TRGB), then these systems have about 10 TRGB bina- 
ries with the right properties per galaxy. These are all 
fainter than ~ 18 and would be fairly expensive in SIM 
observing time, but binaries in the Sagittarius dwarf, the 
Magellanic Clouds and possibly the Sculptor and Fornax 
dwarfs are observable. 

15. IMAGING WITH SIM 

SIM will demonstrate synthesis imaging at optical 
wavelengths in space, thereby showing the viability of 
this approach for the design of the next generation of 
UV/Optical/IR imaging telescopes for space astronomy. 
The telescopes we now have in space at these wavelengths 
(GALEX, HST, & Spitzer) are aU of the "filled aper- 
ture" type, which ultimately limits their angular reso- 
lution to that dictated by the size of the largest rock- 
ets available for launching them. High sensitivity is au- 
tomatically provided with such high-resolution instru- 
ments, but precision collecting area is expensive and may 
not always be required. Synthesis imaging using phase- 
stable UV/Optical/IR interferometers provides, for the 
first time, the possibility of separately choosing the reso- 
lution and the collecting area of space telescopes in order 
to provide a more cost-effective match to the specific as- 
trophysical problems to be addressed. Such flexibility is 
essential for the future if instruments are to provide ever- 
increasing angular resolution and still be affordable. SIM 
will break new ground by demonstrating these imaging 
techniques at optical wavelengths in space. 

Since SIM must have very high fringe phase stabil- 
ity over long periods of time in order to accomplish its 
astrometric goals, it behaves at optical wavelengths in 
much the same way as a ground-based radio interferom- 
eter. The results of the observations can be put into the 
same general form of complex fringe visibilities. Data 
can be obtained over a range of baseline orientations, 
mimicking the rotation of the Earth for ground-based 
synthesis instruments. SIM will also have two baselines 
available, one at 9.0-m and the other at 7.2-m, and data 
for imaging observations can be taken at both baselines 
(although not simultaneously). Figure [T7b shows a typ- 
ical coverage of the aperture (m, v) plane which can be 
achieved. A software simulator for this mode of SIM is 
available, based on t he initial work on this subject by 
iBoker fc AlTenl ()1999l ). The current version of the SIM 
imaging simulator, imSIM, incorporates the latest infor- 
mation on the expected performance parameters of SIM, 
and can generate simulated observations for a variety of 
source models. 

Figure [T7] shows the image which could be produced 
by SIM on a (model) source field consisting of 5 stars 
ranging in brightness from V = 15 — 17, and the {u,v) 
coverage of the simulated data. The stars were observed 
for a total of 300 seconds of on-source integration time 
(30 points, 10 sec per point). The point spread function 
of HST (convolved with the pixels of the HRC) is shown 
for comparison. 

For the imaging demonstration, targets of modest com- 
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Fig. 17. — Upper panel - Example of the (u, »;)-coverage for 
a single channel at mean wavelength A = 500 nm. Both the long 
(~ 9.0m) and short {~ 7.2m) baselines are included, and data taken 
at increments of baseline orientations of 10° and 15°, respectively. 
Only half of the («, v)-plane needs to be observed; the rest can 
be computed using the fact that the field to be imaged is purely 
a real function. Note that SIM can provide its own zero-spacing 
data as well. The (u, i;)-coverage for longer wavelength channels 
would lie inside these circles, and the synthesized PSFs would be 
proportionally wider. Lower panel - A model cluster of 5 stars 
observed for ten seconds at each of the 30 {u, v) points shown in 
the left panel. The image has been GLEANed and restored with a 
gaussian beam of FWHM 0.010 arcsec, and the residuals added in. 
Contours are logarithmic at intervals of 1.25 mag. The brightest 
object (at 0,0) is a 15th mag AO star, the faintest member of the 
cluster (at 0, 0.025) is a 17th mag GO star. Other objects on this 
image are noise and artifacts of the restoration. The FWHM ~ 0.06 
arcsec of the HST/AGS/HRG camera PSF is drawn in (the circle 
at the lower left) for comparison. 
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plexity but showing a wide range of surface brightness 
will be chosen, consistent with the limited number of 
physical baselines available in the SIM instrument. The 
observations for this demonstration will be carried out 
soon after the launch of SIM during the in-orbit check- 
out of the instrument. Subsequently, this capability will 
be available for general use. 

15.1. Performance features of SIM as a synthesis 
imaging instrument 

In addition to its high resolution, SIM as an imager 
has several other novel features. First, the instrument is 
a simple adding interferometer instead of a correlation in- 
terferometer (common in radio astronomy), so that the 
total flux of photons received in the FOV is retained 
in the measurement. If that value is not significantly 
contaminated by photons from extraneous objects in the 
FOV, it can provide an estimate of the "zero-spacing" 
data, as indicated by the dot at the origin in Figure flTk. 
Second, and perhaps most significant, images made syn- 
thetically with SIM will have exceptionally high dynamic 
range. This is a result of the extraordinarily high phase 
stability required of SIM in order to achieve its goal of re- 
peatable microarcsecond astrometric precision. Since this 
stability means that the theoretically-computed point- 
spread function (PSF) is an excellent approximation to 
the actual PSF, the former can be used to remove bright 
stars from the image, leaving only photon noise and low- 
level artifacts from residual phase instability. In marked 
contrast to the situation with filled apertures, these lat- 
ter instabilities will be entirely uncorrelated from one 
measurement location in the (m, v) plane to the next. 
The result is that the remaining low-level artifacts make 
equal contributions everywhere over the synthesized im- 
age rather than piling up around the location of the 
bright stars. After PSF subtraction, the remaining faint 
sources can be detected anywhere in the difference images 
with equal sensitivity, even as close as one interferome- 
ter fringe (« 11 mas at A = 500 nm) from a bright star. 
Preliminary results from our simulations indicate that 
we can expect the dynamic range (defined as the mag- 
nitude difference between the peak value on the bright- 
est star and the faintest detectable companion) on SIM 
images to routinely exceed 6 mag for target stars with 
V < 13, which is comparable to that achieved at radio 
wavelengths with the VLA without any additional data 
processing. However, in contrast to the case with the 
VLA where values of dynamic range as high as 10 mag 
(factor 10'') have been achieved by modeling the residual 
instrumental effects, the dynamic range on SIM's syn- 
thetic images will ultimately be limited by the photon 
noise from the stars in the FOV. With SIM's modest col- 
lecting area, it will be difficult to significantly improve 
the dynamic range on SIM images much further without 
large amounts of observing time. 

15.2. Synthetic imaging science with SIM 

Although intended primarily as a technology demon- 
stration, the synthesis imaging mode of SIM will also 
open new possibilities for science applications on tar- 
gets of high surface brightness, and where resolving 
structure beyond the hmits of the HST/ACS/HRC PSF 
would be significant to an understanding of their nature. 



iBoker fc AIIctI (fl999l ) were the first to consider this ques- 
tion in some detail, and concluded that targets as varied 
as the cores of dense Galactic globular clusters, AGNs 
in galaxies, and dust disks around nearby stars were 
feasible. Based partly on their analysis, a proposal to 
map the stellar distributions and kinematics of the nu- 
clear regions of M 31 was included in the in itial science 
program for SIM (|Unwin fc Turvshevl[200l . Since that 
time, instrument design changes have reduced the num- 
ber of baselines available, and the nulling capability has 
been removed. These two de-scopes of the SIM design 
have rendered the study of dust disks around nearby stars 
significantly more difficult. However, the M31 project 
appears to be still feasible, as are observations of com- 
pact high-surface-brightness targets in general. The ex- 
ample of observations of a dense stellar cluster described 
in the previous section makes it clear that SIM will have 
significant things to say about such targets. 

The successful demonstration of synthesis imaging at 
optical wavelengths with SIM is expected to have a sig- 
nificant impact on the design of future space-based in- 
struments for astronomical imaging. SIM imaging will 
also offer new capabilities for the study of the structure 
of targets presently only barely resolved by HST, or oth- 
erwise confused by the inability to adequately remove 
diffraction spikes and scattered light from nearby bright 
stars. The design of SIM is now quite stable, and many 
details of just how the instrument will operate are now 
known. A new review of potential targets for imaging 
science with SIM is under way. 

16. MEASURING THE PPN PARAMETER 7 

Our current lack of understanding of quantum grav- 
ity and dark energy demands new physics. Searches 
for gravitational waves and attempts to understand the 
nature of dark matter a re further motivated by the 
disco very of dark energy (|Turvshev. Shao fc Nordtvedti 
l2006h . Promising theoretical models involve new 
gravitational effects that differ from General Rela- 
tivity (GR), some of which could manifest them- 
selves as violations of the equivalence principle, varia- 
tion of fundamental constants, modification of the in- 
verse square law of gravity at short distances, Lorenz- 
symmetry breaking, and large-scale gravitational phe- 
nomena. These effects are amenable to study with 
space-based experim ents. For example, sc alar-tensor ex- 

brane- world 
and modi- 
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tensions of gravity (iDamour et al 
gravitational models^ (jPvali et al 
fied gravity on large scales, motivate new searches 
for deviations from a level of 1 0^~^ below the level 
currently tested by experiment (Turvshev et al.l 120041 : 
iTurvshev. Shao k Nordtvedtll2006i r 

The param eterized-post-Newtonian (PPN) parameter 
7 (| Willi l200"6l ) is currently measured to differ from unity 
by no more than (2.1 ± 2.3) x 10~^, as obtained us- 
ing radio-m etric tracking data received fro m the Cassini 
spacecraft (jBertotti. less fc Tortoral [2003D during a so- 
lar conjunction experiment. This accuracy approaches 
the region where multiple tensor-scalar gravity mod- 
els, consi stent with the recent cosmological observa- 
tions (Spe rgel et al.|[2007l ). predict a lower bound for the 
present value of this parameter at the level of (1 — 7) '--^ 
10-6 _ 10-7 (jPamour et al.ll2002 f). motivating the pro- 
posal of space-based experiments to improve the mea- 
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surement of 7. 

SIM will routinely operate at about this level of accu- 
racy - the Sun produces an astrometric bend of 4 mas 
at 90°, more than 1000 times the accuracy of individual 
SIM grid stars. Gravitational effects in the Solar System 
must therefore be included into the global astrometric 
model and the corresponding data analysis, formulated 
such that 7 is a measured parameter. By performing 
differential astrometric measurements with an accuracy 
of a few /ias over the instrument's FOV of 15° SIM will 
provide this precision as a by-product of its astrometric 
program and could measure the parameter 7 with accu- 
racy of a few parts in 10^ (|Turyshevl l2002f ) . a factor of 
10 better than the Cassini result. A precision SIM mea- 
surement would aid the search for cosmologically relevant 
scalar-tensor theories of gravity by looking for a remnant 
scalar field in today's Solar System. 

17. CONCLUSIONS 

The Space Interferometry Mission (SIM) represents a 
revolutionary step forward in the development of astron- 
omy's most ancient measurement technique, determining 
the positions of celestial objects. With an absolute pre- 
cision of 3 /ias (microarcseconds) covering a wide range 
of magnitudes (from V < mag to F > 19 mag), SIM 
will measure positions, distances, and proper motions for 
almost every type of astronomical object with a parallax 
accuracy of 10% or better across the entire Milky Way. 
Within a few parsecs, the parallax precision will be better 
than 1%. These precision measurements will challenge 
our theories of stellar (§[ni[S]) and galactic structure (§[7]), 
while adding to our understanding of dark matter and 
cosmology through refined knowledge of the distance to 
Cepheids (§[5]), the value of the Hubble Constant f^ fli)) . 
the ages and distances of Globular Clusters (§ [5]), and 
the motions of dwarf spheroidal and Local Group galax- 
ies (§ [T^ . Precision SIM astrometry even promises to 
probe small scale phenomena within the cores of quasars 
and active galactic nuclei, including the possible presence 
of binary supermassive black holes (§ [T3| . 

With a differential astrometric precision of 1 microarc- 
second, SIM will measure the orbital motions of many 
classes of multiple objects, including every sort of normal 
and exotic star from low-mass stars (§ [5]) to X-ray bina- 
ries, evolved AGB stars, neutron star and black hole bi- 
naries (§111), and microlensing systems f^i fTU)) . From these 
dynamical measurements, we will be able to determine 
masses for astronomical objects, perhaps the most funda- 
mental single parameter in understanding their physical 
nature. Coupling new physical information from SIM 
with previously known characteristics of such systems, 
such as ages and metallicity, will allow a new level of 
understanding of the evolution of astronomical objects. 

From a cultural standpoint and in the eyes of the pub- 
lic, as well as in the research community, SIM's most dra- 
matic contribution may come from opening up the next 
level of performance in the reconnaissance of our nearest 
stellar neighbors by looking for planets with masses that 
are equal to or just a few times above the mass of our 
own Earth and located in the range of orbits conducive 
to development of habitable environments (§ As de- 
scribed in ^ \2.4l within a deep 5 year survey, SIM could 
find planets with a mass of 1 Af^ orbiting within the 
habitable zones of over 100 of the most favorable stars. 



SIM is the first step in a long term strategy for search- 
ing for other habitable worlds suitable for subsequent 
follow-up with direct detections of planet-light at visible 
or infrared wavelengths. SIM also offers dramatic new 
capabilities for finding planets orbiting stars not acces- 
sible to radial velocity or transit studies, such as young 
stars very massive stars, or stars with highly variable 
photospheres or weak spectral lines. SIM will investigate 
the formation and migration of planets in young plane- 
tary systems, finding out whether Jupiter-mass planets 
are common or rare when orbiting stars younger than a 
few million years old (§ E]). 

As we noted in the §[1] this paper is intended to high- 
light some of the many areas of astrophysics that preci- 
sion astrometry will address, specifically those problems 
amenable to flexibly targeted observations at very high 
precision. As a guide to possible opportunities, we show 
an overview of the assignment of SIM observing time in 
Table m Key Project teams have the task of optimizing 
their science within their allocations, which requires a 
careful trade of the numbers of targets, their magnitudes, 
and the astrometric accuracy needed for the science ob- 
jective. We note that for faint targets, say V > 18, 
accuracies of < 10/xas require hours of mission time, and 
such targets will be selected with substantial care; many 
of the Sections in this paper discuss the considerations 
in target selection. As this process of target definition 
for the Science Team is still underway, we show only the 
formal allocation to the Team. The remaining time will 
be allocated via a future peer-reviewed Guest Observer 
(GO) program, and will be completely open with respect 
to science topic. The number of selected programs will 
be set by NASA; it is expected that extended mission 
would be entirely open to competition. We show a set of 
hypothetical, but representative, programs that together 
would complete the SIM time assignment. The first GO 
program will likely be assigned to further searches for 
terrestrial planets to supplement the NASA Key Projects 
(see §[2]) . The number of targets, magnitudes, and accu- 
racies in Table [H] should be viewed as illustrative, rather 
than the final word on target selection. Observing time 
calculations depend on the science objective (e.g., paral- 
lax and proper motion require different optimizations). 

The technology pioneered by the SIM mission, the first 
long baseline Michelson Interferometer in space, repre- 
sents an important investment for the future of space as- 
tronomy. Visible- and infrared- imaging on the milliarc- 
second scale demands widely separated apertures. SIM's 
long-baseline interferometric capabilities (§ [T5)) will make 
simple images with 10 milliarcsec resolution on y = 15 
objects. Future telescopes, potentially separated by hun- 
dreds of meters, will build upon the technology demon- 
strated by SIM. 

The scientific measurements that SIM will provide can- 
not be duplicated by other means. As a pointed observa- 
tory, SIM will complement the all-sky survey planned for 
the Gaia mission by enabling order-of-magnitude more 
precise observations of both very faint and very bright 
objects. The combination of SIM and Gaia will move as- 
tronomy and astrophysics into a new era of precision dy- 
namical and kinematic knowledge unprecedented in the 
2,000 year-old history of astrometric measurements. 

ACKNOWLEDGMENTS 



40 Unwin et al. 



TABLE 6 

Strawman Key Project and Guest Observer Time Assignment 



ScIghcg Progr3,rn 


Number 


Target V 


Observing 


Acci_ir3,cy 


Mission 




of Targets 


Magnitude 


Mode 


(/ias)^ 


Fraction 


Prime Mission 










1.00 


ociGnce _LeciiTi IvGy Projects 


~ lUUUU 




W ide/ IN arrow 




0.36 


Ciuest Observer Oall 1 


~ 75 


6-9 


Narrow 


0.6-1.2 


0.27'' 


(Terrestrial planets) 












Guest Observer Call 2 


~ 10000 


12 


Wiae 


8 


0.05 


(Open) 












Astrometric grid^ 


1302 


9-10.5 


Wiae 


3 


0.24 


Engineering 










0.08 


Extended Mission 










1.00 


Guest Observer Call 3 










0.7 


(Sample programs) 










0.2^ 


GO Program 1 


~ 180 


6-9 


Narrow 


1.5 


GO Program 2 


- 5000 


12 


Wide 


5 


0.1 


GO Program 3 


~ 20000 


12 


Wide 


8 


0.1 


GO Program 4 


~ 12000 


14 


Wide 


10 


0.1 


GO Program 5 


~ 750 


18 


Wide 


12 


0.1 


GO Program 6 


~ 1500 


18 


Wide 


20 


0.1 


Astrometric grid 


1302 


9-10.5 


Wide 


2.5"= 


0.2 


Engineering"^ 










0.1 



Note. — The SIM instrument and operations are designed for a 10-year total lifetime. Prime mission has a duration of 5 years; 
extended mission, a further 5 years. Science Team Key Project time was assigned by NASA Announcement of Opportunity in 
2000. All Guest Observer (GO) Programs will be competed. The study of terrestrial planets described in §[2]assumes that the first 
GO Call will be devoted to extending the target list for planet searches. GO Call 3 is represented here by programs which span a 
range of magnitudes and accuracies, and which are purely intended to be illustrative. 

In wide-angle mode, accuracy is mid-epoch position, one axis, at end of prime mission. In narrow-angle mode, accuracy is for 
a single measurement in a local reference frame, one axis'' Mission fraction is for 100 visits per target in each of two orthogonal 
axes"^ Grid allocation also includes ~ 50 quasars (~ 1.5% mission fraction) to define the absolute reference frame (see Appendix 
[C)l '' Includes instrument calibration and scheduling margin"^ Estimated accuracy of mid-epoch position, one axis, at end of extended 
mission 
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APPENDIX 
NARROW ANGLE ASTROMETRY 

SIM is designed to perform narrow angle (relative) astrometry in three different modes where it is capable of 
measuring periodic motions with a precision of 0.6 //as per axis in a narrow-angle frame. The modes include the 
standard Narrow Angle (NA) mode, Gridless Narrow Angle Astrometry (GNAA) mode (jShaklan. Milman Sz PanI 
l200 3f). and Grid-based Differential Astrometry (GBDA) mode (jPan et al. 2005). Because of its high efficiency, NA 
mode will be used for the vast majority of narrow angle observations. Early mission results can be obtained with 
GNAA, while GBDA is useful once grid stars are known to the level of 0.1 arcsec. 

NA mode is linked to the wide-angle grid campaign. It uses the end-of-mission grid star positions to determine the 
instrument baseline orientation and length at the time that a target and reference stars were observed. Two roughly 
orthogonal baseline orientations are needed to make a 2-dimensional measurement of the target star motion relative 
to the reference frame. Because all measurements are tied to the absolute grid, the target parallax and proper motion 
are determined with high precision and their errors do not contribute significantly to the single-measurement error. A 
single visit typically requires < 30 min to achieve 0.6 /xas precision. 

GNAA is akin to traditional relative astror netry with a sing le telescope. Relative astrometry typically employs a 
photographic plate, a CCD, or Ronchi Ruling ()Gatewoodlll987l ) to record the positions of a target star and a reference 
stars over several months to years. The reference stars are used to anchor a least-squares conformal transformation that 
matches the scale, rotation, translation, and potentially higher order field-dependent and time-dependent terms into 
a common reference frame, which is then applied to the target. GNAA mode combines measurements of a target, at 
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least four reference stars, and at least three baseline orientations (jShaklan. Milman &: Panll2003[ ). and uses a conformal 
model to solve for instrument parameters (baseline length and orientation, and a phase-delay constant term). Since 
the parallaxes and proper motions of the target and reference stars are not known to microarcsecond precision, this 
approach is not useful for long-period (>0.5 yr) signals. 

GBDA is a hybrid mode that uses coarse (early mission) grid star positions to solve for baseline parameters. GBDA 
requires the measurement of several grid stars using two nearly-orthogonal baseline orientations in addition to the 
target and reference star measurements. It is more efficient than GNAA (two vs. three baselines) but is still limited 
to the detection of short periods. In GBDA, grid catalog errors combine with stochastic baseline repointing errors to 
cause baseline orientation, length, and constant-term errors. These in turn mix with reference frame catalog errors 
resulting in random noise on the target star estimate. To keep the errors between the target and reference stars below 
a fraction of a microarcsecond, the baseline re-orientation requirements for subsequent visits are 0.1 degrees (around 
the line of sight) and 0.01 degrees (tilt toward the target). These are easily met using star trackers. 

To demonstrate the GBDA approach, we modeled the detection of the known planet orbiting Tau Boo (|Butler et al.l 
|1997() . The planet has a measured period of 3.3 days and (minimum) mass 3.9Mjup, from which we derive the 
astrometric amplitude is 9.0 /ias. In the model we assumed that the stellar reflex motion was on one axis in order to 
demonstrate that the magnitude of the astrometric errors is independent of the signal amplitude. Tau Boo's proper 
motion was assumed to be known to 1 masyr"^. We modeled 20 GBDA measurements over a 10 day period. We 
used realistic catalog errors for the reference and target star positions and proper motions. The NA reference frame 
consisted of the 10 brightest reference stars within a 75 arcmin radius of Tau Boo (Fig. [18k). We modeled instrument 
systematic errors based on the flight instrument design, and added signal-dependent noise added in quadrature to the 
target and reference star observations. For our adopted observing scheme, this yields an astrometric error of 0.65 fias 
per axis per visit. Figure 118b shows the simulated signal and measurements, and that the Tau Boo planet is easily 
detected by the GBDA technique. 
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Fig. 18. — Simulation of Grid-based Differential Astrometry (GBDA) with SIM of the 3.9A/jup planet around Tau Boo. (a) 10 brightest 
reference stars within 75 arcmin of Tau Boo. (b) Simulated measurement of the 9 fias amplitude signal (assumed to be along one axis). 
Residual proper motion and parallax are removed from the plot but were included as error sources in the astrometric model. 



ALL-SKY ASTROMETRIC GRID 

SIM observes stars sequentially within a 15° patch of sky termed a 'tile', while the instrument is held inertially stable 
by observing a pair of guide stars with guide interferometers. One guide star lies close in angle to the target star; the 
other is roughly 90° away. Details of the operatio n of SIM, and how the guide stars are used to derive microarcsecond 
precisions on science targets, is discussed in, e.g.. iLaskinl (|2006f ). Here we describe the development of the astrometric 
grid of stars to which the science measurements are referenced, the astrophysical selection of those stars, and the 
'frame tie' of the grid to an inertial frame defined by distant quasars. 

SIM will construct an all-sky astrometric grid including 1302 pre-selected stars and a smaller number of quasars. 
The grid accuracy is expected to be about 3 /j,as in mean-epoch position, based on the current best estimate of 
instrument performance. Each grid object will be observed about 200 times during the nominal 5-year mission. The 
limited number of grid objects and the moderate density of the grid is explained by its special role that is more 
utilitarian than research-driven. SIM will largely rely on self-calibration of the instrument and on the determination 
of the baseline orientation from its own interferometric measurements of the grid objects interspersed with routine 
measurements of science targets. A number of key instrument parameters, for example, the baseline length, can not 
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be determined to the required picometer-level accuracy by external metrology techniques and, therefore, should be 
derived from observations of stars. 

Since the number of instrument and attitude parameters is large, it is necessary to observe a global grid of stars 
multiple times during the mission to construct an overdetermined system of linear equations. The grid solution is a 
one-step direct Least Squares adjustment of a system of ~ 300, 000 linear equations of ~ 160, 000 unknowns. Of the 
latter set, only some 6,500 unknown parameters are related to the grid objects (i.e., the mean positions, parallaxes and 
proper motions that constitute the actual astrometric grid) . An efficient and fast algorithm has be en developed to solve 
the gr id equations eliminating of the numerous attitude terms by QR factorization, described in iMakarov fc MilmanI 
(|2005l ). 

The astrometric accuracy of the grid can be directly evaluated from the global covariance, as far as random errors 
of measured positions are concerned. Systematic errors are not always possible to predict or to model, and are less 
amenable to straightforward covariance analysis. However, the manageable size of the reduced design matrix allows us 
to employ rigorous mathematical analysis of possible systematic errors and compute the com plete space of those pertur- 
bations that have a strong adverse effect on the grid accuracy. As explained in more detail in M akarov &: MilmanI ((2005), 
measurement errors (accidental and systematic alike) propagate non-uniformly in different singular vectors spanning 
the parameter space, which is likely to result in large-scale (so-called 'zonal') correlated spatial errors in the grid. 
They may be represented by spherical orthogonal functions for the parallax error distribution, a nd by vector spherical 
harm onics for the proper motion error distribution, which is a vector field on the unit sphere (|Kopeikin fc Makarovl 
l2006l ). The spatial power spectrum of grid errors is 'red', meaning the signal is larger at low frequencies. 

We explored the grid performance for a wide range of possible instrument systematic errors developed as a by-product 
of the detailed flight instrument design. This allows us not only to evaluate the expected RMS grid performance and 
the properties of the zonal errors, but also to estimate confidence levels, for instance, how likely the actual SIM grid 
will be better than a given value. 

An effective way to mitigate these zonal errors is to use a number of quasars as grid objects; quasar parallaxes are 
vanishingly small and can be constrained to zero in the global solution. Numerical simulations and covariance analysis 
have shown that with only 25 additional optically bright quasars in the grid program, the grid parallax astrometric 
accuracy improves by ^ 28 %, meeting the mission goal (4 /ias) with considerable margin. Using the current best 
estimate for the instrument performance, the grid should achieve 3 /uas. Quasars also dramatically improve the grid 
confidence: the 99% confidence limit on grid parallax error drops from 6.2 /Lias to 3.1 ^as (with the same simulated 
mission but with quasars included in the grid). 

Analytical considerations showed that the astrometric accuracy on bright, frequently observed science program stars 
should be equal to or slightly surpassing that of grid stars. Spreading observation times evenly across the mission time 
is not the best strategy for parallax. Significantly better parallax can be obtained for a given star by simply scheduling 
the chosen number of observation at the most favorable times, determined by a semi-empirical optimization algorithm 
implemented in the wide angle processing code. In a similar way, the observation schedule can be optimized for any 
object on proper motion performance, or a combination of all three types of astrometric parameters. A realistic SIM 
schedule takes into account the desired science objectives of the measurements, and optimizes across the ensemble of all 
science targets both the integration time invested and the observation epochs. While this ideal is likely impossible to 
realize in practice, observation planning will definitely play a major role in extracting the best astrometric performance 
from the instrument. 

ASTROMETRIC REFERENCE FRAME 

SIM will define a reference frame accurate to 3 fia.s using 50-100 quasars to 'tie' the stellar grid to a presumed inertial 
frame. The SIM frame will complement the Gaia catalog, as the SIM frame will have much more precise positions 
than Gala's, albeit for a smaller number of sources. These distant sources are assumed to have negligible parallaxes 
and proper motions, and as noted above, quasar observations will also reduce the amplitude of zonal errors in the 
astrometric grid. Because of the physical effects discussed in §[131 radio-quiet quasars are to be preferred, as there 
are less likely to be asymmetries, or worse, time-dependent position shifts, in the optical photocenters relative to the 
massive black holes that serve as inertial reference points fii ll3.2p . Most SIM science programs are not affected by 
frame rotation. For instance, parallaxes are quite insensitive to a rotation rate. Some topics require the frame to be 
inertial to ~ 2/iasyr~^. Examples include galactic structure, including tidal tails (§[9|) and extragalactic rotational 
parallaxes (§[13]). This requirement should be easily met using 50-100 radio-quiet quasars, selected from bright quasar 
catalogs. Note that a tie to the radio (ICRF) frame plays no direct role in these investigations. 

Tying the SIM inertial frame (defined by radio-quiet quasars and grid stars) to the ICRF serves two fundamental 
purposes. First, future astrometric catalogs will require a registration with the ICRF since that frame is the basis 
for all current astrometric data. Second, accurate registration between reference frames is important for science that 
combines astrometric data in two frames, for instance SIM and VLBI data. This requires a tie to the ICRF, and hence 
SIM will include (radio-loud) quasars in the astrometric grid. The study of non-thermal radio emission associated 
with stars will lead to a better understanding of the mechanisms giving rise to an individual source's spectral energy 
distribution (see ii l6.6p . We showed in §[T3|that astrometry is a powerful tool for probing the structure of blazars. A 
tie to radio (VLBI) images is important for a full understanding of the physical processes in quasars that manifest 
themselves as astrometric motion. 

The tie must of course be done using radio loud quasars, since the ICRF is fundamentally a radio reference frame. 
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Fig. 19. — Histogram of the simulated period distribution of stellar companions to candidate SIM grid stars, showing the effectiveness 
of radial velocity (RV) screening in eliminatin g binary stars with significan t astrometric signatures from the grid. The binary mass and 
orbit distributions are assumed to be those of IDuguennov fc Mavon 1 119911 ). Gray symbols represent the input period distribution. RV 
filters virtually all binaries with periods less than 104 days, leaving only the long-period companions (white symbols). Very few of these 
remaining (i.e., missed) long-period companions produce astrometric motion detectable by SIM (jitter > 4/xas RMS, shown in black). 



The accuracy of the radio-optical frame tie wiU therefore depend on the quasars in common between the two frames. 
Since the SIM reference frame will be much more accurate than the ICRF, errors in tying the frames are likely to 
be dominated by the VLBI observations that define the ICRF. Quasars will be chosen from the 212 ICRF sources; 
there is approximately one ICRF quasar per 15° diameter SIM tile, with fewer in the southern hemisphere and in the 
plane of the Milky Way Galaxy. By establishing an accurate link between the optical SIM frame and the radio ICRF, 
high resolution imaging data at these disparate wavelengths can be accurately aligned to allow absolute positional 
correlation. 

Quasars with radio emission may not make ideal astrometric reference targets, due to variability. It is not yet clear 
how the known structure changes seen at milliarcsecond levels (from VLBI imaging) in the most variable radio-loud 
objects translate into astrometric shifts in the optical. Section [T3] discusses the physical models that allow us to make 
estimates of the effects, prior to SIM launch. Although there is a lot of experience with the ICRF defined by radio 
observations, the physics is not understood well enough to make a simple prediction. Indeed, we fully expect that SIM 
will provide key insights into this problem. 

The quality of the radio-optical frame tie can be improved in a number of ways. First, the quasars in common 
should be selected to avoid the most highly variable targets. Second, repeated astrometric measurements taken during 
the SIM 5-year mission can be checked for consistency and obvious outliers rejected. Third, the number of quasars 
in common should be large to allow us to average out the offsets, which will be uncorrelated from target to target. 
Since the investment of SIM observing time may be traded between number of quasars and accuracy of the astrometric 
measurements, intrinsic astrometric variability would lead one to prefer to increase the number of quasars in common; 
this would also make the identification of outliers easier. Those objects which show astrometric signatures as large 
as perhaps 15 — 20 /xas will be easily eliminated from the reference frame program but will make for very interesting 
astrophysical studies. Finally, we note that detailed studies of the physical processes in quasars, as probed by the 
astrometric program described in §[131 can provide insights into the most effective selection of quasars to use in the 
end-of-mission solution for the reference frame, allowing some quasars to be selectively omitted from the solution. 

SELECTION OF SIM GRID STARS 

The SIM grid comprises 1302 K giant stars spread quasi-uniformly over the sky. K giant stars were chosen for 
several reasons. K giant stars are numerous and located at all galactic latitudes. K giant stars are intrinsically much 
brighter than dwarf stars. Compared to similarly numerous F and G dwarfs (the most common dwarfs in a y « 9 — 12 
magnitude limited sample), K giants are intrinsically brighter by 4 to 5 magnitudes at V. Thus, for a given brightness, 
K giant stars are 5—10 times more distant than F and G dwarfs. Astrometric motions induced by unseen stellar 
and planetary companions are of cou rse minimized by the increased distance (e.g., fCpuld 2001: IfVink et al] l2001l : 
[Peterson. Liu &: Portegies Zwartl[200l) . 

Candidate K giant stars were selected from the Tycho-2 and 2MASS catalogs based upon a compiled catalog of 
photometric colors {BVJHKg), proper motions, and distances (if known) of the stars. A near-infrar ed reduced proper 
motio n diagram was constructed to separate candidate K giant stars from dwarf stars of similar color (| Gould fc MorganI 
120031 ) . The visual extinction for each candidate K giant star was est imated to produce extinction corrected magnitudes 
and colors (i.e., spectral types) and distance estimates ()Ciardil l2004). The Tycho-2 catalog giant star extraction was 
supplemented with a ground-based s urvey, sparsely coverin g the whole sky with a roughly uniform pattern of 'bricks', 
each covering w 0.5 sq-deg per brick ([Patterson et al.ll2001f ). The survey utilized the Washington M, T2, and DD051 
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filter system to identify candidate low metallicity K giant stars (|Maiewski et al.ll2000f) . Together the two methods 
provided a total of ^ 170, 000 K giant stars from which the grid star candidates could be selected. For each 'brick', six 
candidate grid stars were selected for a total of 7812 candidate grid stars. A minimum distance of 500 pc was required 
for candidate selection. The candidate grid stars have a median visual magnitude of = 9.9 mag, with 90% of the 
stars between V — 9.0 — 10.5 mag. The median distance of the grid stars is d = 600 pc, with 90% of the stars between 
d = 500 - 1000 pc. 



— 1 1 I I I r 

La Silla Observatory 



Gaussian lu widLh — 51 m/s 
Binsize - 15 m/s 





-1000 



in-IlJ QlLO I In n I i-inl-i 




-1 I 1 I I I I 

KV Difference beLweeii 2 
observations separated 
by 300+ days 

X = 847 Pairs 

|AHV|". 1000 m/s: 9% 



JZl m n I n n 



-500 500 

Radial VclociLy Difference | m/s 



Fig. 20. — Frequency distribution of the difference between consecutive radial velocity measurements separated by more than the 300 
days. Data taken at La Silla Observatory as part of the SIM grid star RV monitoring program. 



Detailed simulation s of the effectiveness of radial velocity (RV) 'screening' of grid star candidates has been performed 
bv lCatanzaritd (|2004f ). using reasonable assumptions about the population of binary companions that induce astromet- 
ric perturbations in the primary star. The results show that without RV screening, the grid is likely to be significantly 
contaminated, but that screening to a fairly modest precision of 50 m s~^, with 3 observations spread over at least 
3 years, are sufficient to reduce the contamination to a very low level (Fig. [T9|) . RV screening mostly fails to reject 
brown dwarfs and giant planets, but most of the companions remaining produce small astromet ric signatures, and 
have a negligible effect on the grid as a whole. A pilot study of the RV stability of K giant stars bv lFrink et al.l ()2001l ) 
showed that early K giant stars {B — V 0.8 — 1.1 mag) were intrinsically stable at a level of 20 — 30 m s~^ . Analytic 
studies and Monte Carlo simulations ljGouldll2001l : iFrink et al.ll200ll ; iPeterson. Liu &: Porte gies Zwart"2002') indicated 
that a modest RV screening program (30-50 m s~^) would reduce the fraction of grid stars with umnodclcd astrometric 
motions detectable by SIM to less than 1 — 10%. In September 2004, a high precision (< 50 m s~^) RV monitoring 
campaign was begun to remove those stars from the candidate list with unseen stellar companions. Observations are 
currently being performed at the 1.2m Euler Telescope at La Silla. The RV program spans four years and is structured 
such that half of the candidates (3 per brick) are observed twice in years 1 and 2 with a minimum separation between 
observations of 9 months. Those stars which display large RV excursions {xt > 4) are removed from the candidate 
list. The removed stars are replaced and are observed during years 3 and 4, along with the surviving candidates from 
years 1 and 2. The program has finished its second year. 

The effective sensitivity of the radial velocity program to detect companions is a convolution of the stellar at- 
mospheric jitter and the instrumental jitter. Thus, the effective precision of the RV measurements is given by 
<^totai = V i'^inst + '''*) where (Jinst is the RV precision of the instrument/data reduction, and (t^, is the intrinsic 
stellar atmospheric jitter. In Figure 1201 we present the frequency distribution of the difference between two consecu- 
tive RV measurements separated by 300-1- days for 847 candidate grid stars. The precision for each RV measurem ent 
is approximately ainst = 20 — 30 m s~^. The intrinsic stellar jitter of early K giants, as found bv lFrink et al.l ()2001[ ). is 
a^, « 25 m s^^, yielding an expected distribution width (in the absence of companions) of Utotai ~ 30 — 35 m s~^. The 
central portion of the distribution is well modeled with a gaussian of width Utotai = 51 m s""'^, slightly larger than the 
expected width, indicating the presence (not unexpectedly) of binaries in the sample or a slightly higher atmospheric 
jitter (40 vs. 25 m s""'^). Additionally, the high velocity wings (|ARV| > 100 m s~^) are likely the result of the 
presence of unseen companions to the grid stars. Using a reduced chi-square of < 4 as the cut-off for acceptance, 
approximately 60% of the stars in this sub-sample would pass through for further radial velocity vetting in years 3 
and 4 of the program. 
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